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Abstract 

We present a collection of signatures for physics beyond the standard 
model that need to be explored at the LHC. The signatures are orga- 
nized according to the experimental objects that appear in the final state, 
and in particular the number of high px leptons. Our report, which in- 
cludes brief experimental and theoretical reviews as well as original re- 
sults, summarizes the activities of the "New Physics" working group for 
the "Physics at TeV Colliders" workshop (Les Houches, France, 11-29 
June, 2007). 
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Parti 

Introduction 

G. Brooijmans, A. Delgado, B.A. Dobrescu, C. Grojean and M. Narain 

The exploration of the energy frontier will soon enter a dramatic new phase. With the 
startup of the LHC, planned for later this year, collisions at partonic center-of-mass energies 
above the TeV scale will for the first time be observed in large numbers. 

The Standard Model currently provides an impressively accurate description of a wide 
range of experimental data. Nevertheless, the seven-fold increase in the center-of mass energy 
compared to the current highest-energy collider, the Tevatron, implies that the LHC will probe 
short distances where physics may be fundamentally different from the Standard Model. As a 
result, there is great potential for paradigm-changing discoveries, but at the same time the lack 
of reliable predictions for physics at the TeV scale makes it difficult to optimize the discovery 
potential of the LHC. 

The TeV scale has been known for more than 30 years to be the energy of collisions 
required for revealing the origin of electroweak symmetry breaking. The computation of the 
amplitude for longitudinal WW scattering [1] shows that perturbative unitarity is violated un- 
less certain new particles exist at the TeV scale. More precisely, either a Higgs boson or some 
spin-1 particles that couple to WW (as in the case of Technicolor or Higgsless models) are 
within the reach of the LHC, or else quantum field theory is no longer a good description of 
nature at that scale. 

Given that ATLAS and CMS are multi-purpose detectors, it is commonly believed that 
they will provide such an in-depth exploration of the TeV scale that the nature of new physics 
will be revealed. Although this is likely to be true, one should recognize that the backgrounds 
will be large, and an effective search for the manifestations of new physics would require a large 
number of analyses dedicated to particular final states. Other than the unitarity of longitudinal 
WW scattering, there are no clear-cut indications of what the ATLAS and CMS experiments 
might observe. Furthermore, recent theoretical developments have shown that the range of 
possibilities for physics at the TeV scale is very broad. Many well-motivated models predict 
various new particles which may be tested at the LHC. Hence, it would be useful to analyze 
as many of them as possible in order to ensure that the triggers are well-chosen and that the 
physics analyses have sufficient coverage. 

The purpose of this report is to provide the LHC experimentalists with a collection of 
signatures for physics beyond the Standard Model organized according to the experimental 
objects that appear in the final state. The next four sections are focused on final states that 
include, in turn, three or more leptons, two leptons, a single lepton, and no leptons. Section 6 
then describes an interface for event generators used in searches for physics beyond the Standard 
Model. 

Whatever the nature of TeV scale physics is, the LHC will advance the understanding of 
the basic laws of physics. We hope that this report will help the effort of the particle physics 
community of pinning down the correct description of physics at the TeV scale. 



Multi-Lepton Final States 
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Part 2 

Four leptons + missing energy from one 
UED 

M. Gigg and P. Ribeiro 
Abstract 

Minimal Universal Extra Dimensions (MUED) models predict the pres- 
ence of massive Kaluza-Klein particles decaying to final states con- 
taining Standard Model leptons and jets. The multi-lepton final states 
provide the cleanest signature. The ability of the CMS detector to 
find MUED final state signals with four electrons, four muons or two 
electrons and two muons was studied. The prospect of distinguish- 
ing between MUED and the Minimal Supersymmetric Standard Model 
(MSSM) is then discussed using simulations from the event generator 
Herwig-i-i-. 

1. DISCOVERY POTENTIAL FOR THE FOUR LEPTON FINAL STATE 
1.1 Introduction 

The Universal Extra Dimensions (UED) model [2] is an extension of the sub-milUmeter extra 
dimensions model (ADD) [3, 4] in which all Standard Model (SM) fields, fermions as well 
as bosons, propagate in the bulk. In the minimal UED (MUED) scenario [5] only one Extra 
Dimension (ED) compactified on an orbifold is needed to create an infinite number of excitation 
modes of Kaluza-Klein (KK) particles with the same spin and couplings as the corresponding 
SM particles. The mass spectrum of the KK particles is defined by three free parameters: R~^, 
the size of the ED, given in terms of the compactification radius; AR, the number of excitation 
modes (KK levels) allowed in the effective theory; and mn, the SM Higgs boson mass. KK 
partners are indicated with the subscript related to the n-th mode of excitations (e.g. at the first 
level they are gi, Zi, uli, cri, 71). A direct search for MUED in the multi-lepton channel at 
Tevatron energy of 1.8 TeV [6] set a lower bound on the size of ED of > 280 GeV. Also, 
constraints from dark matter infer 600 < < 1050 GeV [7]. 

In this section a summary report on the discovery potential of the CMS experiment [8] for 
MUED is presented. The complete analysis is described in [9]. The experimental signatures for 
production of first level KK states at hadron colliders are isolated leptons and/or jets radiated 
in the cascade decay process, in addition to the transverse missing energy carried away by the 
lightest KK particle (LKP). These characteristics were exploited to discriminate the signal from 
the background. The four lepton final state constitutes the cleanest channel. The KK mass 
spectrum, however, is highly degenerate since the masses of the KK particles with respect to 
the corresponding SM particles at tree level are m\ = n^/ i?^ + nig]^, where n is the excitation 
mode. Furthermore, radiative corrections do not introduce an additional large splitting and 
typically, within the same excitation mode, there is a difference of about ICQ GeV between the 
heaviest and the lightest KK particle. Therefore, the average values of the lepton momentum 
and the missing transverse energy are typically smaller than average values which characterise 
searches for supersymmetric events. 
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1.2 Signal and background processes 

The MUED signal is produced in a pp collision as a pair of two KK strongly interacting parti- 
cles, gluons (gi) or quarks (qi). Three significant subprocesses were considered: 

pp^gm, PP ^ Qi/qiQi/qi, pp^giQi/Qi, 

Singlet and doublet KK quarks of the first generation were taken into account. Four points of 
the MUED parameter space have been chosen for the study: nin = 120 GeV, AR = 20 and 
e {300, 500, 700, 900} GeV. The total cross section strongly depends on the compactifi- 
cation radius being equal to 2190, 165, 26 and 5.86 pb for R-^ = 300. 500. 700 and 900 GeV 
respectively. The four lepton final state signature can provide a discrimination against the SM 
background and is considered as in the following: 

gi^QiQ, Qi^Z^Q, Z.^L^l^, L,^LKP{-f,)e^. (1) 

The KK gluon (gi) decays into a KK quark (qi) and a SM anti-quark; then, the qi decays 
into the KK boson (Zi) and a SM quark. Subsequently, Zi decays into a pair of leptons, one 
being a KK lepton (li = singlet 1ri or, mainly, doublet 1li). Finally, li can decay only into 
the LKP photon (71) and a SM lepton. Instead, if qi is produced initially, then the decay 
cascade is shorter. The B.R. of the four lepton final state is about 10""^ — 10~^. Within a 
decay branch the pair of SM leptons (i'^P) has the same flavour and opposite sign. Three 
possible combinations of four leptons arise, namely 4e, 4// and lel/j,, studied in three separated 
channels. Signal events were generated with CompHEP with particle definitions and Feynman 
rules taken from [10] at the LO approximation. The background to MUED signals results 
from SM processes with four leptons in the final state. The dominant sources are from the 
continuum production of (Z*/7*) (Z*/7*) and real ZZ production, from processes involving pair 
production of heavy quark flavours such as tt and bbbb, and the associated production of Zbb. 
Background events were generated with PYTHIA and ALPGEN. Signal and background events 
were processed with full detector simulation using official CMS software (OSCAR version 
3.6.5). Underlying events from minimum bias interactions were superimposed to generated 
events, assuming an average number of 5 inelastic colUsions, including diffractive interactions, 
at each beam crossing, simulating the effect of pile-up at C — 2 ■ 10^^cm~^s~^ (LHC low 
luminosity scenario). The reconstruction of physics objects was based on the dedicated CMS 
software ORCA (version 8.7.3/4). 

1.3 Event Selection 

First, Level 1 (LI) and High Level trigger (HLT) requirements for £ = 2 ■ 10^"^cm^^s~^ are 
applied to the simulated events. We then require the presence of at least two pairs of OSSF 
leptons. The leptons should be isolated (4 iso) and are required to be within the following 
kinematical boundaries (52): 

• electrons with px >7.0 GeV and \r]\ < 2.5, 

• muons with pt >5.0 GeV and \r]\ < 2.4. 

Because a substantial fraction of the background leptons results from b-quark leptonic de- 
cays, we reject events where one or more b-jets are identified (Bveto). Due to the soft KK mass 
spectrum, leptons from the MUED cascade (eq. 1) have on average lower transverse momentum 
than some of the background channels, like for example the background from top quark decays. 
For this reason we apply upper bound cuts on the lepton transverse momentum (lept pt) of 70, 
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60, 40, 30 GeV for the P*, 2""^, 3'^'*, 4*^* lepton sorted in pt, respectively. A missing transverse 
energy cut of '^t > 60 GeV proves to be important especially for high values where the 
is higher due to the massive LKPs, as the background is significantly rejected with respect 
to a small reduction of signal events. Finally, we apply a selection on the invariant mass of the 
lepton pairs, according to which an event is rejected if it has one or more OSSF lepton pair with 
Minv < 5 GeV or Minv > 80 GeV, aimed at rejecting the ZZ background. The selection cuts 
were chosen so that the signal efficiency is maximum for (R~^ = 900 GeV), where the signal 
cross section is lowest. 

The summary of all selection cuts is presented in figure 1 in two ways: as an efficiency 
of each cut after the previous one (left), and as a cross section after each cut (right). After all 
selection cuts, a S/B greater than five is achieved for all studied points of the parameter space. 

1.4 Results 

The CMS discovery potential of MUED in the four lepton channel, defined as the integrated 
luminosity needed to measure a signal with a significance (Scp) of five standard deviations is 
shown in figure 2. The significance estimator ScP gives the probability to observe a number 
of events equal or greater than Nobs = Nsignai + Nekg, assuming a background-only hypothe- 
sis, converted to the equivalent number of standard deviations of a Gaussian distribution. The 
dashed (solid) lines show results including (not including) systematical uncertainties. The sys- 
tematic uncertainties include a 20% uncertainty on the background cross section, the effect of 
jet energy scale on the missing energy distribution (3-10%, pr dependent) and a 5% uncertainty 
in the b-tagging algorithm efficiency. For the three four-lepton channels ( 4e, 4/i and 2e2fx) 
and for the integrated luminosity of 30 fb^^ , the signal significance is above the background 
by a few standard deviations and therefore the MUED signal could be detected at the CMS 
experiment during the first few years of data taking. In the 4/x and the 2e2// channels alone, a 
significance of five standard deviations for — 500 GeV could be reached with less than one 
fb-^ of data. 

2. MODEL DISCRIMINATION 

The four lepton channel described above can also produced in supersymmetry (SUSY). The 
focus of this section will be in comparison of different signatures from SUSY and UED using 
the monte carlo event generator Herwig++. 

As of version 2.1 of Herwig++ [1 1] BSM physics was included for the first time with both 
the minimal supersymmetric standard model (MSSM) and minimal universal extra dimensions 
models implemented including spin correlations in production and decay [12].^ This allows 
the comparative study of both models within the same general purpose event generator. To 
compare the two models in the most sensible manner the mass spectra should be the same. 
The simplest way to achieve this is to chose a scale for the MUED model and then adjust the 
parameters in the MSSM so that the relevant masses are matched. Two scales were chosen for 
this work, R""^ = 500 GcV and R"^ = 900 GcV, both with AR = 20 which produced the mass 
spectra shown in tables 1 and 2. The MSSM spectrum file and decay tables were produced 
using SDECAY version 1.3 [13]. 

The previous section tells us that the four lepton channel could give a sizeable signal 
compared to the standard model background. Plotting the invariant mass distribution of the 



'All plots are made with version 2.1.1 ofHerwig++ which included some minor bugfixes. 
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qLi 


uri 


dm 


hi 


71 


626.31 


588.27 


576.31 


574.90 


514.78 


500.98 



Table 1: The UED mass spectrum for R ^ = 500 GeV and KR = 20. The SUSY counterparts are matched to 
this. All values are in GeV. 
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9li 


um 


dm 


Ili 


7i 


1114.25 


1050.50 


1028.84 


1025.28 


926.79 


900.00 



Table 2: The UED mass spectrum for R ^ = 900 GeV and AR = 20. The SUSY counterparts are matched to 
this. All values are in GeV. 

di-muon pairs in a four muon final state event for MUED and SUSY gives the results shown in 
figure 3. It is apparent that larger values of the masses within the spectrum make it more difficult 
to distinguish between the two models. Moreover, in the = 500 GeV case the shapes of the 
distributions are similar it is just the overall number of events in the SUSY case that is larger 
due to the size of the relative branching ratios. Given the possibility that the distributions could 
be so similar it will be necessary to make use of other combinations of invariant mass plots. 
The most logical is the invariant mass of a quark plus one of the lepton or antileptons.^ Since it 
is possible to distinguish between leptons and antileptons in a detector one can make separate 
distributions for the jet^ plus lepton and the jet + antilepton cases. Since these now take into 
account the helicity of the quark these distributions will be more sensitive to spin effects. 

Figure 4 shows the invariant mass of a jet plus a lepton while figure 5 shows the distribu- 
tion for a jet plus an antilepton for the two scales under consideration. Again there is a greater 
difference at a lower value of the compactification radius where it would seem that the shapes of 
the distributions differ more at higher invariant mass values. The main reason, however, for the 
similarity in the distributions under study is that they have combined effects from opposite sets 
of spin correlations. This result can be attributed to firstly the lack of distinction between quark 
and antiquark in the jet/lepton distributions which means that two sets of data with opposite 
spin correlations appear on the same plot thereby cancelling the effect out. Also the run was 
set up so that both left and right-handed partners to the quarks were produced in the initial hard 
collision and when these decay they will have, again, opposite correlations. It is these kinds of 
effect that will cause the most trouble in trying to distinguish between the two models. 

A useful quantity in trying to achieve this at the LHC will be the asymmetry, defined as 




where dP/dmji+ and dP/drriji- are the antilepton and lepton distributions respectively. Its 
usefulness stems from the fact that the LHC is a proton-proton collider and will produce an 
excess of quarks over antiquarks. The result will be a slight favour in one helicity mode over 
the other meaning the asymmetry should be the most sensitive to the underlying physics model. 
The distributions from Herwig-i-i- are shown in figure 6. 

^The theoretical distributions for the case where one distinguishes between quark and antiquark are given in 
[14] and will not be reproduced here. 

^We are working at the parton level so we define a jet as simply a quark or an antiquark. 
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3. CONCLUSIONS 

The CMS experiment will be able to detect evidence of MUED model in the four lepton final 
state up to = 900 GeV with an integrated luminosity of 30 fb~^. For the purpose of dis- 
crimination between the MUED and SUSY scenarios it is apparent that the analysis of the four 
lepton signature alone is insufficient. The best hope is using asymmetry distribution as this is 
most sensitive to the spin differences in the underlying physics model. 
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Figure 1: Cumulative efficiency and cross section after each selection cut for the MUED signal (R ^ G 
{300,500,700,900} GeV, nin = 120 GeV, AR = 20) and the background for all channels. Only the upper 
statistical uncertainties are shown. 
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Figure 2: The discovery potential of MUED signals (R-^ e {300, 500, 700, 900} GeV, nin 120 GeV, AR ^ 20) 
in the four-lepton channels is defined as the integrated luminosity needed to measure a signal with a significance 
(Scp) of five standard deviations. The dashed (solid) lines show results including (not including) systematical 
uncertainties. The uncertainties due to the limited understanding of the detector performance and characteristic of 
the early phase of the LHC data taking are not considered and may limit the sensitivity below one /6^^ (horizontal 
'First data uncertainty' line). 
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Figure 5: The invariant mass of a quark or an antiquark with an antUepton for (a) R ^ = 500 GeV and (b) 
R-i = 900 GeV. Solid: UED, dot-dash: SUSY. 
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Figure 6: The asymmetry, as defined in equation 2, for (a) R-^ = 500 GeV and (b) R-^ = 900 GeV. SoUd: UED, 
dot-dash: SUSY. 
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Part 3 

LHC Events with Three or More Leptons 
Can Reveal Fermiophobic Bosons 

R.S. Chivukula and E.H. Simmons 

1. INTRODUCTION 

Events with three or more leptons plus either jets or missing energy can lead to the discovery 
of fermiophobic W bosons associated with the origin of electroweak symmetry breaking. One 
possibility is the process pp {W*) W'^*^Z WZZ jji+£~e+e- where the W 
is assumed to decay hadronically and i can be an electron or muon. Another is the process 
pp — > WZjj where the W and Z re-scatter through the W resonance; the final state of interest 
here includes three leptons, two jets, and missing energy. This section describes a general class 
of "Higgsless" models that include fermiophobic W bosons, specify the particular model used 
in our phenomenological studies, and then describe the calculations and results for each multi- 
lepton channel in turn. 

The pp — > WZjj signal is the classic l^VF-scattering process studied for a strongly inter- 
acting symmetry breaking sector, with the W playing an analogous role to the technirho boson. 
In the three-site higgsless model considered below it is possible to calculate this process in a 
fully gauge-invariant manner, rather than using the traditional method that involves separately 
calculating the signal (by using a model of tttt scattering in conjunction with the effective W 
approximation) and background (usually done by considering the standard model with a light 
Higgs boson). 

2. HIGGSLESS MODELS IN GENERAL 

Higgsless models [15] provide electroweak symmetry breaking, including unitarization of the 
scattering of longitudinal W and Z bosons, without employing a scalar Higgs boson. The most 
extensively studied models [16, 17] are based on a five-dimensional SU{2) x SU{2) x U{1) 
gauge theory in a slice of Anti-deSitter space, and electroweak symmetry breaking is encoded 
in the boundary conditions of the gauge fields. Using the AdS/CFT correspondence [18], these 
theories may be viewed as "dual" descriptions of walking technicolor theories [19-24]. In addi- 
tion to a massless photon and near-standard W and Z bosons, the spectrum includes an infinite 
tower of additional massive vector bosons (the higher Kaluza-Klein or KK excitations), whose 
exchange is responsible for unitarizing longitudinal W and Z boson scattering [25]. To provide 
the necessary unitarization, the masses of the lightest KK bosons must be less than about 1 
TeV. Using deconstruction, it has been shown [26] that a Higgsless model whose fermions are 
localized {i.e., derive their electroweak properties from a single site on the deconstructed lattice) 
cannot simultaneously satisfy unitarity bounds and precision electroweak constraints. 

The size of corrections to electroweak processes in Higgsless models may be reduced by 
considering delocalized fermions [27-29], i.e., considering the effect of the distribution of the 
wavefunctions of ordinary fermions in the fifth dimension (corresponding, in the deconstruction 
language, to allowing the fermions to derive their electroweak properties from several sites 
on the lattice). Higgsless models with delocalized fermions provide an example of a viable 
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effective theory of a strongly interacting symmetry breaking sector consistent with precision 
electroweak tests. 

It has been shown [30] that, in an arbitrary Higgsless model, if the probability distribution 
of the delocalized fermions is related to the W wavefunction (a condition called "ideal" delocal- 
ization), then deviations in precision electroweak parameters are minimized. Ideal delocaliza- 
tion results in the W resonances being fermiophobic. Phenomenological limits on delocalized 
Higgsless models may be derived [31] from limits on the deviation of the triple-gauge boson 
(WWZ) vertices from their standard model value; current constraints allow for the lightest KK 
resonances to have masses as low as 400 GeV. 

3. THREE-SITE MODEM IN PARTICULAR 

Many issues of interest, such as ideal fermion delocalization and the generation of fermion 
masses (including the top quark mass) can be illustrated in a Higgsless model deconstructed 
to just three sites [32]. The electroweak sector of the three-site Higgsless model incorporates 
an SU{2)o x SU{2)i x U{1)2 gauge group, and 2 nonlinear {SU{2) x SU{2))/SU{2) sigma 
models responsible for breaking this symmetry down to U{l)em- The extended electroweak 
gauge sector of the three-site model is that of the Breaking Electroweak Symmetry Strongly 
(BESS) model [33]. The mass-eigenstate vector bosons are admixtures of the seven gauge- 
bosons in SU{2Y X U (1), with one massless photon, three corresponding to the standard model 
and Z, and three nearly-degenerate W'"^ and Z'. For the reasons described above, the 
masses of the W'^ and Z' in the three-site model (and indeed for the lightest bosons in any 
Higgsless model with ideal fermion delocaUzation) must be between roughly 400 GeV and 1 



The left-handed fermions are doublets coupling to the two SU (2) groups, which may be 
correspondingly labeled ipio and ipLi- The right-handed fermions are a doublet coupling to 
SU (2)i, and two singlet fermions coupled to U{1)2, denoted ur2 and dR2 in the case of 
quarks. The fermions V'lo> '4>li, and V'ri have U (1) charges typical of the left-handed doublets 
in the standard model, +1/6 for quarks and —1/2 for leptons. Similarly, the fermion ur2 has 
U (1) charges typical for the right-handed up-quarks (-1-2/3), and dR2 has the U (1) charge asso- 
ciated with the right-handed down-quarks (—1/3) or the leptons (—1). With these assignments, 
one may write the Yukawa couplings and fermion mass term 



Here the Dirac mass M is typically large (of order 2 or more TeV), is flavor-universal and 
chosen to be of order Mw/Mwn to satisfy the constraint of ideal delocalization, and the Er 
are proportional to the light-fermion masses (and are therefore small except in the case of the 
top-quark) [32]. These couplings yield a seesaw-like mass matrix, resulting in light standard- 
model-Uke fermion eigenstates, along with a set of degenerate vectorial doublets (one of each 
standard model weak-doublet). 

The three-site model is sufficiently simple that it is possible to implement the model 
in CompHEP [34] or MADGRAPH [35], and - by including both fermion and gauge-boson 
couplings - to do so in a way that is fully gauge-invariant. 



TeV. 




(1) 
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4. 4 LEPTONS PLUS 2 JETS 

LHC events with 4 charged leptons and 2 jets can reveal [36] the presence of a fermiophobic W 
boson that is produced through pp — > WZZ with the Z bosons decaying leptonically and the 
W decaying hadronically. The signal for the W boson comes from associated production pp — > 
(W*) W' Z followed by the decay W WZ. The backgrounds include the irreducible 
SM background pp WZZ jji^i^i^i^; a related, but reducible, SM background pp 
ZZZ jji^i^i^i^ in which the hadronically-decaying Z is mis-identified as a W; and 
all other SM processes leading to the same final state pp jji^i^i'^i^ through different 
intermediate steps, including processes in which one or more of the jets is gluonic. 

Ref. [36] has calculated the full signal and background in the context of the three-site 
higgsless model [32] using the cuts described here. One set of cuts is used to suppress the SM 
backgrounds: 

Mjj = 80 ±15 GeV, AR{jj) < 1.5, $^]?r(^) + 5ZPr(j) = ±15 GeV. (2) 

z j 

The first of these selects dijets arising from on-shell W decay (leaving a margin for the exper- 
imental resolution [37]); the second reflects the dijet separation of the signal events; and the 
third exploits the conservation of transverse momentum in the signal events. In addition, a set 
of transverse momentum and rapidity cuts are imposed on the jets and charged leptons 

p„ > 10 GeV, \7]^\ < 2.5 , Ptj > 15 GeV, \7]^\ < 4.5 . (3) 

for particle identification. 

These cuts essentially eliminate the first two sources of background and reduce the third to 
a manageable size, with the signal peak standing out cleanly. Ref. [36] concludes that fermio- 
phobic W bosons will be visible in this channel at the LHC throughout their entire allowed 
mass range from 400 — 1200 GeV. The integrated luminosity required for detecting the W in 
this channel is shown here as a function of W mass in Fig. 1. 

5. 3 LEPTONS PLUS MISSING ENERGY 

LHC events with 3 charged leptons, missing energy, and forward jets can reveal the presence 
of a fermiophobic W boson produced through the scattering process pp — > WZqq', where 
both vector bosons decay leptonically [38^0]. In this case, the signal arises from the W 
contribution to the vector boson fusion subprocess WZ — > WZ. 

An initial estimate of the W' signal and related backgrounds was presented in in ref. [41] 
for a 5d higgsless sum rule scenario. Ref. [36] has improved on this by performing the first 
tree-level calculation that includes both the signal and the full electroweak (EW) and QCD 
backgrounds for the 2 — > 4 scattering process pp WZjf in the context of a complete, gauge- 
invariant higgsless model (the three-site model [32]). A forward-jet tag is used to eliminate the 
reducible QCD background [38] from the annihilation process qq — > WZ. The irreducible 
QCD backgrounds pp — > WZjj with jj — qg, gg serving as forward jets are suppressed by 
the cuts 

E- > 300 GeV , pj,^ > 30 GeV , |r;,- 1 < 4.5 , | Arj-^ | > 4 , (4) 

where Ej and pj<j are transverse energy and momentum of each final-state jet, r]j is the forward 
jet rapidity, and \Ar]jj | is the difference between the rapidities of the two forward jets. The cut 
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on |A?7jj | is especially good at suppressing the QCD backgrounds pp — > WZgg, WZqg in the 
low Mwz region [42]. The following lepton identification cuts are also employed 

Pti > 10 GeV, \%\ < 2.5. (5) 

While one must specify a reference value of the SM Higgs boson mass in computing the SM 
EW backgrounds, the authors of [36] found that varying the Higgs mass over the range Mh = 
115 GeV - 1 TeV had httle effect. 

Ultimately, ref. [36] reports both the signal and backgrounds for the transverse mass dis- 
tribution of the vector boson pair, where M^{WZ) = KP{m) + p%{m) + \pf'-^\f - 
\prp{M) + p^^^^^\'^. Counting the signal and background events in the range Q.ShMw' < Mt < 
l.ObMw', yields the integrated LHC luminosities required for 3a and 5a detections of the W 
boson in this channel, as shown here in Fig. 1. 
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Figure 1 : Integrated luminosities required for Scr and 5a detection of W signals as a function of Miy> . The dotted 
and dashed-dotted curves are for the WZZ channel, while the dashed and solid curves are for the WZjj channel. 
From ref. [36]. 



6. CONCLUSIONS 

Both the pp W'Z WZZ jjii and pp W'jj WZjj i^Sijj channels are 
promising for revealing fermiophobic W' bosons, like those in Higgsless models, at the LHC 
[36]. The WZZ channel has a distinct signal with a clean resonance peak. The WZjj' channel 
has a larger cross section when Mw' is heavy, but the measurement is complicated by the 
missing Et of the final-state neutrino. Hence, confirming the existence of both signals for the 
W boson, as well as the absence of a Higgs-like signal in pp — *■ ZZqq — > 4^ qq, will be strong 
evidence for Higgsless electroweak symmetry breaking [36]. In particular, as shown here in 
Fig. 1, for Mw' = 500 (400) GeV, the 5cr discovery of W' requires an integrated luminosity of 
26(7.8)fb-^ for: pp WZZ jjA£, and 12(7)fb-i f or pp WZjj uSijj. These are 
within the reach of the first few years' run at the LHC. 
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Low-Scale Technicolor at the LHC 

G. Azuelos, K. Black, T. Bose, J. Ferland, Y. Gershtein, K. Lane and A. Martin 
Abstract 

If technicolor is responsible for electroweak symmetry breaking, there 
are strong phenomenological arguments that its energy scale is at most 

a few hundred GeV and that the lightest technihadrons are within reach 
of the ATLAS and CMS experiments at the LHC. Furthermore, the 
spin-one technihadrons pt, and ot are expected to be very nar- 
row, with striking experimental signatures involving decays to pairs of 
electroweak gauge bosons (7, W, Z) or an electroweak boson plus 
a spin- zero ht- Preliminary studies of signals and backgrounds for 
such modes are presented. With luminosities of a few to a few tens 
of femtobams, almost all the spin-one states may be discovered up to 
masses of about 600 GeV. With higher luminosities, one can observe 
decay angular distributions and technipions that establish the underly- 
ing technicolor origin of the signals. Preliminary ATLAS studies show 
that, with 50-100 fb~^ and assuming Maj, — l.lMp^, both processes 
a| ^ Z^W^ (with Mp^ ~ 500 GeV) may be seen in the E+^-E^Pt 
final state and p^, ^ Z°7rJ (up to Mp^ ~ 400 GeV) in £+£-6jet, 
where l — e, ji. 

1. INTRODUCTION 

Technicolor (TC) is a proposed strong gauge interaction responsible for the dynamical break- 
down of electroweak symmetry [43,44]. Modern technicolor has a slowly-running ("walking") 
gauge coupling [19,22,45,46]. This feature allows extended technicolor (ETC) [47] to gen- 
erate realistic masses for quarks, leptons and technipions (ttt) with the very large ETC boson 
masses (10^-10^ TeV) necessary to suppress flavor-changing neutral current interactions. (For 
reviews, see Refs. [48,49].) The important phenomenological consequence of walking is that 
the technicolor scale is likely to be much lower and the spectrum of this low- scale technicolor 
(LSTC) much richer and more experimentally accessible [50-52] than originally thought [53]. 
The basic argument is this: (1) The walking TC gauge coupling requires either a large number 
Nd of technifermion doublets so that Ktc ^ 250 GeV/y/No ^100 GeV, or two TC scales, 
one much lower than 250 GcV.' (2) Walking enhances tct masses much more than those of their 
vector partners, px and ujt- This effect probably closes the all-TTx decay channels of the lightest 
techni- vectors. In LSTC, then, we expect that the lightest pt and ljt lie below about 0.5 TeV 
and that they decay to an electroweak boson (7, W, Z) plus ttt; a pair of electroweak bosons; 
and //, especially t^l~. These channels have very distinctive signatures, made all the more so 
because pT and lot are very narrow, T{pt) ^ 1-5 GeV and T{ujt) — 0.1-0.5 GeV. Technipions 
are expected to decay via ETC interactions to the heaviest fermion-antifermion flavors allowed 
kinematically, providing the best chance of their being detected. 

^For an alternate view based on a small TC gauge group, SU{2), see Refs. [54, 55]. 
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Many higher-mass states are reasonably expected in addition to pt and ujt- In Refs. [56, 
57] it was argued that waUdng TC invalidates the standard QCD-based calculations of the 
precision-electroweak ^'-parameter [58-61]. In particular, the spectral functions appearing in S 
cannot be saturated by a single px and its axial-vector partner a^. Thus, walking TC produces 
something like a tower of vector and axial-vector isovector states above the lightest px and 
ot. All (or many) of them may contribute significantly to the S'-parameter.^ Most important 
phenomenologically, in models with small S, the lightest ar and px likely are nearly degen- 
erate and have similar couplings to their respective weak vector and axial-vector currents; see, 
e.g., Refs. [62-66]. The 3 vrr-decay channels of the ot are closed, so these states are also very 
narrow, T{aT) ^ 0.5 GeV. 

The Pt, cot, clt, and tt^ of low-scale technicolor that we consider are bound states of 
the lightest technifermion electroweak doublet, (T[/,Td). The phenomenology of these tech- 
nihadrons is set forth in the "Technicolor Straw-Man Model" (TCSM) [66-68]. The TCSM's 
most important assumptions are: (1) There are Nd isodoublets of technifermions transforming 
according to the fundamental representation of the TC gauge group. The lightest doublet is 
an ordinary-color singlet.^ We use = 9 in calculations; then, the technipion decay con- 
stant Ft ^ 246 GeV/\/ND = 82 GeV. The technipion isotriplet composed of the lightest 
technifermions is a simple two-state admixture, 

\Up') = smx\W^'')+cosx\4''), (1) 

where Wl is a longitudinally-polarized weak boson and ttt is a mass eigenstate, the lightest 
technipion referred to above, and sinx = Ft/246 GeV = l/\fW^. This is why the lightest 
spin-one technihadrons are so narrow: all their decay amplitudes are suppressed by a power of 
sin X for each Wl emitted and by a power of e = sin 9w for each transversely-poldonzQd l^_j_.In 
addition, decays to -kt are phase-space limited."^ The technihadrons' principal decay modes are 
listed in Table 1 . (2) The lightest bound-state technihadrons may be treated in isolation, without 
significant mixing or other interference from higher-mass states. (3) Techni-isospin is a good 
symmetry.^ These assumptions allow the TCSM to be described by a relatively small number 
of parameters. The ones used for the present study are, we believe, fairly generic; they are listed 
below. 

The main discovery channel for low-scale technicolor at the Tevatron is pt W^ttt — > 
i'^uebq, with at least one tagged 6-jet. At the LHC, this channel is swamped by a it background 
100 times larger than at the Tevatron. There the discovery channels will be WZ, 'yW and 7Z, 
with the weak bosons decaying into charged leptons.^ In the TCSM each of these modes is 
dominated (generally ^ 80%) by production of a single resonance: 

p±^iy±Z°, a^^-fW^, ujt^iZ^. (2) 

^These higher mass states are also important in unitarizing longitudinal gauge boson scattering at high energies. 

^Some of these doublets may be color nonsinglets with, e.g., three doublets for each color triplet. Tech- 
nifermions get "hard" masses from ETC and ordinary color interactions and will have some hierarchy of masses. 
We expect that the lightest will be color-singlets. 

^Because the interactions of the techni-vectors with electroweak gauge bosons (and fermions) are suppressed 
by sin x, they can be light, ^ 200 GeV, without conflicting with precision electroweak and Tevatron data. 

^Also, something hke topcolor-assisted technicolor [69J is needed to keep the top quark from decaying co- 
piously into TTy 6 when M^^ ^ 160 GeV. Thus, if tt^ is heavier than the top, it will not decay exclusively to 
th. 

^The channel W^Z'^ was studied by P. Kreuzer, Search for Technicolor at CMS in the pxc W + Z 

Channel, CMS Note 2006/135, and Ref. [70]. To consider the decay angular distributions, we employ somewhat 
different cuts than he did. 
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Table 1: Amplitude factors for the dominant decay modes of px Wlttt, WlWl and pt, (jJt, o.t Gttt, 
GWl [66,68]. Here, Wl is a longitudinally-polarized and G = 7, W±, Z± a transversely-polarized electroweak 
gauge boson. Technifermion charges are Qu — Qd + 1, and sinx = Ft/246 GeV — 1/^/Ne,. Amphtudes 
for Wlttt and WlWl are proportional to gp^ = ^/47ra^, where = 2. 16 {3 /Ntc)'-> those for emission of a 
transverse gauge boson are proportional to e = \/47ra. 



In Sects. 2-4, the PGS detector simulator [71] is used for our preliminary studies of these sig- 
nals and their backgrounds. None of these LHC discovery modes involve observation of an 
actual technipion (other than the ones already observed, and There are other strong- 
interaction scenarios of electroweak symmetry breaking (e.g., so-called Higgsless models in 
five dimensions [26,72] and deconstructed models [15-17,27]) which predict narrow vector 
and axial- vector resonances, but they do not decay to technipion-like objects. Therefore, ob- 
servation of technipions in the final state is important for confirming LSTC as the mechanism 
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underlying electroweak symmetry breaking. It is possible to do this at high luminosity with 
the decays — > Z^tt^ — > t^£~hq. This channel also provides the interesting possibility 
of observing both and in the same final state. This analysis, using ATLFAST [73], is 
summarized in Sect. 5 [74]. 

In addition to the discovery of narrow resonances in these channels, the angular distribu- 
tions of the two-body final states in the techni- vector rest frame provide compelling evidence of 
their underlying technicolor origin. Because all the modes involve at least one longitudinally- 
polarized weak boson, the distributions are 

da{qq ^ ^l^l) da{qq p$ t^t^l) ■ 2 



d cos 6 ' d cos 6 



oc sin^ 9 ; (3) 



da{qq ^a^ ^ Wt) ^ ^^M^i^^^^^ ocl + cos^ ^ . (4) 
d cos B d cos B 

It is fortunate that each of the two-electroweak-boson final states is dominated by a single 

technihadron resonance. Otherwise, because of the resonances' expected closeness, it would 
likely be impossible to disentangle the different forms. Our simulations include these angular 
distributions. 

For Les Houches, we concentrated on three TCSM mass points that cover most of the 
reasonable range of LSTC scales; they are listed in Table 2. In all cases, we assumed isospin 
symmetry, with Mp^ — M^^^ and M^^, = l.lMp^; also, the pt and ot constants describing 
coupling to their respective weak currents were taken equal; sinx = 1/3; Qu + Qd = 1; 
Ntc = 4 for the TC gauge group SU{Ntc); and My,,,., = Ma,,^-, = Mp^ for the LSTC mass 
parameters controlling the strength of pr, ujt, cit decays to a transverse electroweak boson plus 
ttt/Wl or pt/out [66-68]. Pythia [75] has been updated to include these and other LSTC 
processes, according to the rules of the TCSM. The new release and its description may be 
found at www . hepf orge . org. 

The simulations presented here, especially those using the PGS detector simulator [71], 
are preliminary and in many respects quite superficial.^ E.g., no attempt was made to optimize 
S/B in the PGS studies. Nor did we carry out a serious analysis of statistical, let alone sys- 
tematic, errors.^ Still, we believe the simulations establish the LHC's ability to discover, or 
rule out, important signatures of low-scale technicolor. Beyond that, it is our intent that the 
present studies will stimulate more thorough ones by ourselves and by the ATLAS and CMS 
collaborations. 



2. p$ ^w^zi^ £±z/^£+r 

The cross sections for — > W^Z^, including branching ratios to electrons and muons, are 
listed in Table 1? Signal events were generated with the updated Pythia [75]. The principal 
backgrounds are in Table 3, along with their generators and parton-level cuts. The ALPGEN 
backgrounds were passed through PYTHIA for showering and hadronization. 

'All PGS simulations were done using the ATLAS parameter set provided with the PGS extension of MAD- 
GRAPHv4.0 [35, 76]. The relevant parameters are: calorimeter segmentation Ar? x Ac/) = 0.1x0.1, jet resolution 
AE/E = 0.8/Ve, and electromagnetic resolution AE/E = 0.1/\/E+0.01. To model muons more realistically, 
we changed the sagitta resolution to 5{)firn. Vv^ith this set of parameters, the PGS lepton identification efficiency is 
« 90% in their kinematic region of interest, Pt > 10 GeV and < 2.5. 

^Potentially important sources of systematic error are the higher-order QCD corrections to signal and back- 
grounds. The ii' -factors can be quite large, ~ 1.5. 

^For the TCSM parameters we use, about 20% of these —^ W"^ Z° rates involve one transverse gauge boson. 
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Case 




















A 


300 


330 


200 


400 




110 


168 


19.2 


158 


B 


400 


440 


275 


500 




36.2 


64.7 


6.2 


88.6 


C 


500 


550 


350 


600 




16.0 


30.7 


2.8 


45.4 



Table 2: Masses (in GeV) and signal cross sections (in fb) for the lightest technihadrons for the three TCSM mass 
points in this study. Isospin symmetry is assumed. The tt^' is an isosinglet, color-singlet technipion expected in 
TC models; we have assumed it so heavy that the pT, uit and ot cannot decay to it. The cross sections combine 
contributions from px, lot and ar, but tend to be dominated by a single resonance. Branching ratios of the W and 
Z to electrons and muons are included in all cross sections. 



Background 


Cross section (fb) 


Comments 


WZ ^M + u 


430 




zz 


52 




z + bb^ e+e- bb 


7600 


Prib) > 15.0 GeV, 


Vb 


< 3.5 


tt 2e 2iy bb 


22,800 


Pythia generator 




2560 


Pril) > 40 GeV, < 3.5 


Wiet £u-f (fake) 


3180 


PT(jet) > 40 GeV, \r]^,t\ < 3.5 
Includes 0.1% fake rate 




700 


Pril) > 40 GeV, | 




< 3.5 


Zjet ^ £+£-7 (fake) 


315 


PT(jet) > 40 GeV, |r;jet| < 3.5 
Includes 0.1% fake rate 



Table 3: Backgrounds to the W^Z°, ^W^ and 7^" signals of low-scale technicolor. The generator is ALP- 
GENvl3 [77J unless indicated otherwise. Branching ratios of the W and Z to electrons and muons are included in 
the cross sections. 

Events were selected which have exactly three leptons, electrons and/or muons, with two 
having the same flavor and opposite sign, \rii\ < 2.5, at least one having pt > 30 GeV, and 
the others with pr > 10 GeV. No cut on was applied in this analysis, though it may 
improve S/B to do so. The Z was reconstructed from two same-flavor, opposite-sign leptons 

with the smallest \Mg+i- — Mz\ < 7.8 GeV. In reconstructing the W, Pt{i^) = — J2^t was 
assumed, and the quadratic ambiguity in Pz{i^) was resolved in favor of the solution minimizing 
the opening angle between the neutrino and the charged lepton assigned to the W, as would be 
expected for a boosted W}^ 

Figure 1 shows various distributions for case A with Pt{W), Pt{Z) > 50 GeV and 
i?T(jets) = X] £'r(jet) < 125 GeV. The Ht cut significantly reduces the tt background. 
The integrated luminosity is 10fb~^, and a strong signal peak is clearly visible above back- 
ground in the first panel. Fitting the peak to a Gaussian, its mass is 311 GeV. Counting signal 
and background within twice the fitted resolution of 25 GeV, only J Cdt = 2.4 fb~^ is required 
for a S/ y/ S + B = 5a discovery of this resonance. Table 4 contains the final-state mass res- 

'*^The efficacy of this procedure, which was adopted at Les Houches ("the LH algorithm") was compared to a 
"TeV algorithm" used at the Tevatron. The TeV algorithm chooses the Pz(j^) solution which gives the smaller W 
energy. In ATLAS [78] and CMS-based analyses [79], it was found that the TeV algorithm does sUghtly better at 
choosing the correct solution. 
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olutions and 5 a discovery luminosities for the two-electroweak-boson modes considered here. 
The poorer resolution in the WZ and channels is due to the resolution. 



wz 




a (GcV) 


•*^min J 


Pj^ cut 


A 


311 


25.6 


2.4 


v^(W Z) > 50 GeV 


B 


414 


34 5 


7 2 


•On^iW Z) > 75 GeV 


c 


515 


41.0 


14.7 


Pt(W^, ^) > 75 GeV 




Mpeak (GeV) 


a (GeV) 


>C^in(fb-') 


Pt cut 


A 


328 


31.2 


2.3 


Pr(7, W^) > 75 GeV 


B 


439 


39.1 


4.5 


Pt(7,W^) > 100 GcV 


C 


547 


39.3 


7.8 


Pr(7, > 125 GeV 


7Z 


Mpeak (GeV) 


a (GeV) 


C^in (fb-^) 


Pt cut 


A 


299 


7.3 


16.8 


Pt(7, ^) > 80 GeV 


B 


398 


9.4 


45.5 


PTil,Z) > 110 GcV 


C 


498 


12.0 


97.2 


Pt(7, ^) > 150 GeV 



Table 4: PGS simulation data for the spin-one technihadrons decaying to a pair of electroweak gauge bosons. A 
simple Gaussian fit is made to determine the mass and width of the resonance; signal and background events are 
counted within ±2 a of the peak value to determine the minimum luminosity needed for S/V S + B = 5cr. 

The second panel in Fig. 1 shows the total and signal WZ angular (| cos 6'|) distribution. 
The distribution is folded since the signal and WZ background are even functions of cos 9. The 
total distribution reflects the forward-backward peaking of the standard WZ production. The 
signal distribution (open black histogram) is much flatter than the expected sin^ 9, presumably 
because of poorly-fit W's and their effect on determining Mwz and the WZ rest frame. To 
remedy this, we take advantage of the LSTC technihadrons' very small widths and require 
280 < Mwz < 340 GcV; see Fig. 2. The signal distribution now has the expected sin^ 9 shape. 
The remaining large background at | cos6'| ^0.7 can be fit and subtracted by measuring the 
angular distribution in the sidebands 220 < Mwz ^ 280 GeV and 340 < Mwz ^ 400 GeV. 
We believe that 10 fb~^ is sufficient to distinguish this angular distribution from 1 + cos^ 9, but 
detailed fitting is required to confirm this; see Ref. [78]. 

The I cos 9\ distribution in the p^-resonance region is shown in Fig. 3 for cases B (Mp^ = 
400 GcV) and C (500 GcV). ThepriW, Z) cuts of 75 and 100 GcV were chosen to accept signal 
data over the same cos 9 range, 0.0-0.9, as in Case A. The luminosities of 40 and 80 fb~^ were 
chosen to give roughly the same statistics. For the higher luminosities, the effects of pile-up on 
calorimetry and tracking were not considered. 

Finally, for the 50 GeV splitting used in case C, it appears possible to see and as 
separate peaks in the Mwz distribution. This was studied in Refs. [78,79]. With cuts similar to 
those used above (except thatpT(H^, Z) > 50 GeV) a simulation [78] using ATLFAST was per- 

"The pT cuts listed in this table were not optimized for discovery; rather they were chosen partly to reveal 
as much of the angular distributions as possible consistent with background reduction. Presumably, in a real 
search, harder cuts would be employed to reveal the signal. Once it was found, the px cut could be loosened 
and the final-state mass cut tightened to focus on the angular distribution. The upward shift of the peak 
mass, evident in their non-Gaussian high-mass tails, may be due to —^ W^Z^ at about 20% the strength of 

. These issues are being considered with more sophistication using ATLFAST [73] and CMS Fast Simulation 
(https : //twiki . cern . ch/twiki/bin/view/CMS/WorkBookFastSimulation) in Refs. [78,79]. 
They find discovery luminosities about 15-30% lower than estimated here. 
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Figure 1: Signal and background distributions of a 300GeV W^Z^ i^vd+i- for lOfb"^ at the LHC; 

Pt{W, Z) > 50 GcV and iJ^ljcts) < 125 GeV. The open black histograms are the signal contributions. 




Figure 2: WZ angular distribution of the signal and backgrounds for a 300 GeV W^Z'^ —+ i^vii^l for 

10 fb"^ at the LHC; 280 < Mwz < 340 GeV; other cuts are listed in the text. The color code is given in Fig. 1. 



formed. The result is seen in Fig. 4 where a luminosity of 50 fb was assumed. The pr and ax 
were modeled as Gaussian distributions and the background above 300 GeV as a falling expo- 
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Figure 3: W Z angular distributions of the signal and backgrounds for W^Z'^ i^vgi^i Left: Mp^ — 

400 GeV with 380 < Mwz < 440 GeV for 40 fb"\ Right: Mp^ = 500 GcV with 480 < Mwz < 540 GcV for 
80 fb~^. Other cuts are listed in the text. The color code is given in Fig. 1. 



ATLAS PRELIMINARY, 5D1b^ 




Red curves represent the signal resonances, blue the backgrounds, and green the total [78]. 



nential. The ax appears as a high-mass shoulder. The CMS analysis finds a similar result [79]. 
Both pt and ax can be observed and a 5 a (combined) discovery achieved with luminosity 
^ 9.5 fb^^ provided that data is reasonably described by the simulated I^t resolution. 
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Figure 5: Signal and backgi-ound distributions of a 440 GeV —> ^W^ li^vi for 20 fb ^ at the LHC; 
Pt{7, W) > 100 GeV and I^t > 20 GeV. The open black histograms are the signal contributions. 

3. ^ -fWi 

The axial-vector isovector ay is a new addition to the TCSM framework, motivated by the 
arguments that the S parameter problem of technicolor is ameliorated if px and ax are nearly 
degenerate and have nearly the same couplings to the vector and axial-vector weak currents. 

On account of space limitation, we show here only the results of PGS simulation of 
Case B, for which Maj, = 440 GeV. For the decays to a pair of electroweak bosons con- 
sidered in this report, a{aT — > ''yW)B{W are the largest; it is 65 fb in case B. 
Signal and background events were generated with pri'j) > 40 GeV. As noted, the discovery 
search could impose a higher threshold. Events were selected with exactly one lepton, having 
Pt > 10 GeV and \r]\ < 2.5. Distributions are displayed in Fig. 5, in which I^t > 20 GeV 
and Pt{i), Pt{W) > 100 GeV. The principal backgrounds are in Table 3. A jet ^ 7 fake 
rate of 10~^ was assumed for W + jet [37]. Another possible background, 7 + jet where the 
jet fakes a lepton, is negligible after the I^t cut. The luminosity of 20 fb^^ was chosen to give 
reasonable statistics for the signal's angular distribution. Consequently, the resonant peak 
has a significance of 10 a. 

It is clear from these distributions that the backgrounds are a more severe impediment to 
observing the signal's angular distribution (1 + cos^ 6) than they were in the case of pr — > WZ. 
In particular, there is no obvious cut to remove them other than one on M^vk- The result of 



29 



S + B.L = 20fb- 



S + B.L = 40fb-^ 





"0 0.1 0.2 0.3 0.4^Q, 



,0.6 0.7 0.8 0.9 1 



Figure 6: Angular distributions for the signal and backgrounds for a 440 GeV jW^ with 420 < M^w < 

460 GeV (left) and a 300 GeV ivt 7^° with 290 < M^z < 310 GeV (right) at the LHC; other cuts are listed 
in the text. The luminosities are 20 and 40 fb^^, respectively. The color codes are given in Figs. 5 and 7. 



requiring 420 < M^w < 460 GeV is in Fig. 6. Even though the signal's expected forward- 
backward excesses are eliminated by the pril) cut, this clearly is a flatter distribution than 
the sin^ 9 ones above. As in that case, subtracting the background by measuring the sidebands 
should reveal the signal. Careful fitting to see that it is consistent with 1 + cos^ 9 after cuts is 
work for the future. 

4. UT^lZl 

The CUT is as important to find as the pr, with which it is expected to be nearly degenerate, and 
the ax- Yet it is the most challenging to see of the light techni- vectors. At the Tevatron, the 
primary discovery mode is ujt — > 7vr^ — > 766. Backgrounds to this may make this channel 
difficult at the LHC; studies need to be done! Two other channels have much lower branching 
ratios, but are much cleaner: cut — ^ 'jZl and cux — * £^£^ . We discuss the first of 

these in this section.'^ 

As with the search just described, the main backgrounds to cot 7^° are standard 
Z° + 7 and Z + jet(^ 7) production. The cot signal, however, is about 10 times smaller than 
the ttT one (see Table 2), so considerably higher luminosities are required to see a significant 
signal peak and the characteristic 1 +cos^ 9 distribution. Here we discuss the PGS simulation for 
case A, in which M^j, = 300 GeV. Events were selected with two same-flavor, opposite-sign 
leptons, each having px > 10 GeV and rapidity \t]\ < 2.5. The leptons were required to satisfy 
\M£+£- — Mz\ < 7.8 GeV. Distributions are shown in Fig. 7 for a luminosity of 40 fb^^ and for 
Pt{i), Pt{Z) > 80 GeV. Note the much better final-state mass resolution than for pr WZ 
and ax 'jW. The significance of the signal peak is about 8 a; the high luminosity is needed 

'^Preliminary studies of ujt M^M > including its angular distribution have been carried out by J. Butler and 
K. Black. The lut e+e^ mode is ripe for picking. 
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Figure 7: Signal and background distributions of a 300 GeV ojt 7^° ^ 7^+^" for 40 fb ^ at the LHC; 
\Mi+i- — Mz I < 7.8 GeV and pt(7, Z) > 80 GeV. The open black histograms are the signal contributions. 



to accumulate statistics for the angular distribution.^^ 

To expose the angular distribution, we take advantage of the superior 7Z mass resolution 
and impose a tight cut on M^^ of 300 ± 10 GeV. The result is in Fig. 6. Because of the more 
stringent px cuts, the data above | cos6'| > 0.8 are lost. While quite acceptable, the angular 
distribution's signal-to-background is not as favorable as it was for 'jW. As in that case, 

detailed fitting beyond our scope is needed to determine how well the measured distribution fits 
the expectation. And, as there, the backgrounds can be subtracted by measuring the angular 
distribution in sidebands. 



5. a$ Z^n^ e+£-bq 

Even if narrow resonances in the WZ, "fW and 7Z channels are found as described above at the 
LHC, all with nearly the same mass and with the expected angular distributions, it will remain 
essential to discover a technipion to cement the technicolor interpretation of these states. In this 
section we present an analysis of Z^n^ carried out for the ATLAS detector [74]. The 

large backgrounds to this signal require large luminosity. On the other hand, for the masses 
assumed here, this channel has the extra advantage that the rate for — > Z^tt^ is only 2-4 
times smaller than for ^l^T' creating another opportunity for observing both resonant 

'^For case B, the corresponding luminosity is 80 fb^^. 



31 




Figure 8: Left: Mp^ — A/^tt ^^d Ma^ ~ M-^^ for cases A, B, C inATLAS. Right: Mp^^ar ~ -^ttt vs. Af^^ signal 
and background events per 2x10 GeV^ per 100 fb^^ for case A [74]. 



peaks in the same (well, similar) final state. 

As noted, txt are expected to decay into the heaviest fermion-antifermion flavors kinemat- 
ically allowed. For the range of M^rj, considered here, this implies vr^t tb and 7r° bb or it 
are dominant. Actually, something like topcolor-assisted technicolor [69] is needed to produce 
rrit — 175 GeV, and this implies that the coupling of ttt to t-quarks is suppressed by a factor 
of about mh/mt from its naive value. Thus, while the tt^ considered here is massive enough to 
decay to tq, it should still have an appreciable branching fraction to cb and ub. The latter are the 
decay channels considered here. It will be interesting to consider the tq modes. However, they 
are not yet included in Pythia and, therefore, B{n^ — > bq) = 0.87 in the simulation reported 
here. This branching ratio may decrease substantially when the top modes are included. That 
would change the search strategy, but we expect that vrj can be seen in tq as well. 

The signal cross sections are a{pp Z°7r^)5(Z° = (99, 59) (A), 

(71, 17) (B), and (37, 9) fb (C), where the two numbers are approximately the and con- 
tributions. The principal backgrounds and their leading-order cross sections are: it (500 pb) 
and, including the branching ratio of to e+e^ and fi^fi', Z^jj (344 pb), Z^bb (56 pb) and 
Z^bj (11 pb).'"^ See Ref. [74] for generation details. The ATLAS detector simulation used ATL- 
FAST [73]. An additional factor of 90% was applied to the simulation for lepton identification 
efficiency. The 6-jet tag efficiency used was 50%; this corresponds to a light-jet mistag rate 
of 1% and a c-jet mistag rate of 10%. 

To satisfy ATLAS trigger and high-luminosity (100 fb^^ per year) running conditions, 
events were preselected with (1) two same-flavor opposite-sign leptons with pt > 20 GeV and 
(2) at least one 6-tagged jet and one non-6-tagged jet, both with pt > 20 GeV; the two highest- 

'"^Recall footnote 8 regarding systematic errors on such backgrounds. Background contributions from processes 
with even more jets are possible; they are partly accounted for by the leading-log parton showering approximation 
and initial and final state QCD radiation in Pythia. The Z''^bb rate here is much larger than in Table 3 because it 
was generated with greater 6-jet acceptance. 
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Pt jets satisfying these conditions are the TTT-candidate jets. For the p^, these selections resulted 
in 548 (A), 382 (B) and 184 (C) signal events per 100fb"\ For the a^, there were 297 (A), 
117 (B) and 34 (C) events. The total background event numbers, dominated by Zjj and Zhh, 
were: (6930, 10670) (A), (7505, 6285) (B) and (3015, 2550) (C). Here, the first number in each 
pair refers to pr and the second to ot; the background is the number of events in an elliptical 
region in MT^^-{Mp^^aT ~ -^ttt) space centered at the mean and with widths corresponding to 

1.5(7. 

The following cuts were then applied to optimize the signal significances: (1) -^t < 
35 GeV to suppress tt; (2) the highest-pT jet had pr > 80 (A), 115 (B), 150 GeV (C); (3) the 
second highest-pr jet had px > 65 (A), 80 (B), 100 GeV (C); (4) there is exactly one 6-tagged 
jet; and (5) M^+c- = 91 ± 5 GcV. After these cuts, the number of remaining signal events is 
(344, 215) (A), (242,75) (B) and (126,21) (C) for (p^, a^).The backgrounds under these signals 
are (403,900) (A), (346,242) (B) and (96,69) (C). For the parameters used in this simulation, 
then, only in case C is the ar not observable in the Zttt channel in 100 fb~^. 

The resolution in M^^ varies from 16 to 23 GeV and in Mp^^ar from 19 to 30 GeV. Most 
of the error comes from the ttt jets' energy measurements. Therefore, much of it cancels in 
Q = Mp^^ar ~ ^TTT- This is shown on the left in Fig. 8, where the resolution in this difference 
ranges from 3 to 5 GeV. This sharpness will facilitate the discovery of pr, Zt^t and 

other technivector-to-technipion decays. The signals and background for case A are on the 
right in Fig. 8. The twin peaks stand out dramatically (looking rather like Boston's Back Bay). 

In summary: For the TCSM parameters used here, there should be no difficulty seeing 

Pj, and in the Z^tt^ channel in case A, and the and a strong indication of the 0^ in 
case B. In case C, only — > Z'^Hj-. can be seen in 100 fb^^. The minimal cross sections 
(times B(Z ^ e+e" / ji^ jjT)) and luminosities required to see the p^ and signals at 5 cr 
significance are in Table 5. 
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Pt 


29 


28 


14 




8.3 


15 


15 




41 


18 


18 




48 


106 


390 



Table 5: Minimal cross-section times branching fractions (in fb, left) and minimal luminosities (in fb , right) 
required for 5 a significance in cases A, B, C. 



6. CONCLUSIONS AND OUTLOOK 

Low-scale technicolor (with No = 0(10) isodoublets transforming as SU{Ntc) fundamen- 
tals) is a well-motivated scenario for strong electroweak symmetry breaking with a walking TC 
gauge coupling. The Technicolor Straw-Man framework provides the simplest phenomenology 
of this scenario by assuming that the lightest technihadrons — p^, ujt, and ttt — and the 
electroweak gauge bosons can be treated in isolation. This framework is now implemented in 
Pythia. 

'^The Q-value was used to advantage in the most recent CDF search 
for Pt W^TTT ^ e^uebiet; see CDF/ANAL/EXOTIC/PUBLIC/8566, 

http://www-cdf.fnal.gov/physics/exotic/r2a/20061025.techcolor/, and its impor- 
tance was emphasized in Ref. [66]. 
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We used Pythia and (mainly) the generic detector simulator PGS to study the final- 
state mass peaks and angular distributions for the LSTC discovery channels at the LHC: 
W^Z^, — > ^W^ and lot — 7^", with leptonic decays of the weak bosons. We also carried 
out an ATLFAST simulation for p^, ^ Z^n^ i~^i~b jet. The results are very promising. 
For the fairly generic TCSM parameters chosen, the technivector mesons can be discovered up 
to about 500-600 GeV in the two-gauge boson modes, usually with a few to a few tens of fb~^. 
The angular distributions, dispositive of the underlying technicolor dynamics, can be discerned 
with a few tens to 100 fb~^ (except for a higher mass lot — 7-^°)- Taking advantage of the 
superb resolutions in Q = Mp^^ar ~ ^i^d M^+i- for px, aT — > Zttt £~^£~biet, both 
resonances and the technipion can be seen for Mp^, < 500 GeV and Maj, ^ 400 GeV. 

Still, these studies just scratch the surface of what can and needs to be done to gauge the 
potential of the ATLAS and CMS detectors for discovering and probing low-scale technicolor. 
Simulating detector response to the signals and backgrounds of the relatively simple processes 
we considered requires considerably more sophistication, in both depth and breadth, than we 
have been able to deploy. Issues such as the accuracy with which technivector masses and decay 
angular distributions can be determined as a function of luminosity are especially important. 
While we believe that the TCSM parameters — sinx, Qu + Qd, M^^, My. and M^., Ntc — 
we chose are reasonable, relative branching fractions can be fairly sensitive to them, as Table 1 
indicates [68]. It would be valuable to reconsider the processes examined here for a range of 
these parameters. Finally, there are other modes we have not been able to consider but which 
are nevertheless of considerable interest. Two outstanding examples are p° , cot, 
and LOT, Pt ~^ 7^t ~^ 7^^- Thus, the main goal of our Les Houches studies, as it is for the 
other "Beyond the Standard Model" ones started at Les Houches, is to motivate the ATLAS 
and CMS collaborations to broaden the scope of their searches for the origin and dynamics of 
electroweak symmetry breaking. 

"Faith" is a fine invention 
When Gentlemen can see — 
But Microscopes are prudent 
In an Emergency. 

— Emily Dickinson, 1860 
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Technivectors at the LHC 

J. Him, A. Martin and V. Sanz 
Abstract 

Assuming composite spin-1 states to be the most relevant particles pro- 
duced by EW scale strong interactions, we model them with a simple 
parametrization inspired by extra dimensions. Our flexible framework 
accommodates deviations from a QCD-like spectrum and interactions, 
as required by precision electroweak measurements. 

1. INTRODUCTION: 

As was emphasized in the Les Houches non-SUSY BSM working group, very few LHC sim- 
ulations of dynamical electroweak symmetry breaking (DEWSB) scenarios are available. In 
order to remedy this situation, Les Houches 2007 called for an effective description of strong 
interactions, flexible enough to interpolate between some known models of resonance inter- 
actions in 4D or 5D, yet economical enough to have a tractable parameter space. Ultimately, 
such a framework could play the same role for strong interactions as was played by minimal 
Supergravity (mSUGRA) for the case of SUSY: simplifying assumptions reduce the number of 
parameters from O(IOO) down to a few, enabling a slew of phenomenological studies. 

As a step towards an effective DEWSB description, we present a flexible yet manageable 
model of interactions between spin-1 resonances and the Standard Model (without Higgs) using 
the framework of Holographic Technicolor (HTC) [64]. In HTC the number of parameters in 
the effective lagrangian of resonance interactions is reduced by deriving the interactions from a 
precursor 5D lagrangian, as suggested by the AdS/CFT correspondence [80,81]. In practice, we 
see no reason to be constrained by a strict 5D formulation [82] : we simply model interactions 
between 4D resonances, but resort to 5D techniques to compute the parameters. 

HTC is similar to Higgsless [17] models, but contains deviations from pure AdS 5D ge- 
ometry in the form of effective warp factors that differ for the various fields. These warp factors 
are a departure from true 5D modelling [82] but they are motivated by the requirement of small 
deviations from the SM in the gauge sector (oblique corrections [64], and cubic couplings, 
see below). With nonstandard 5D geometry we achieve a different resonance spectrum from 
rescaled QCD, confirming and building upon previous 4D results [56, 62, 83]. 

We also refrain from modelling the fermions in the extra-dimension, as this would not 
reduce the number of parameters: in the present study, the couplings of fermions to resonances 
gffv free parameters, set to pass experimental constraints, while the couplings of fermions 
to Z are assumed to exactly obey SM relations. This can be relaxed in the future. 

In sections 2 and 3 we present the basics of HTC parameterization and the constraints 
on its parameter space. In sections 4 and 5 we describe some LHC signals using a Mad- 
Graph [35]/BRIDGE [84] implementation of HTC. In the future we hope to provide a more 
complete package to allow further study. 
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2. HOLOGRAPHIC TECHNICOLOR: 

In HTC, as in Higgsless models, the SM SU{2)^ and U{l)em gauge fields are the lightest 
Kaluza-Klein (KK) states of 5D SU{2)l ® SU{2)r gauge fields. The higher KK excitations of 
the same 5D fields are interpreted as new spin-1 resonances. Provided they are light enough, 
these spin-1 resonances assist in unitarizing WW scattering in the absence of a Higgs. 

The masses and interactions of the resonances are dictated by the geometry of the 5th 
dimension, which is set by a warp factor w{z). The warp factor appears in the 5D metric ds"^ = 
where the extra coordinate 2; is restricted to the interval/o ^ z ^ h. Instead of 
a single warp factor, in HTC we allow the axial (A) and vector (V) combinations of 5D gauge 

fields to feel different backgrounds. Specifically we define wx = {lo/z)exp 

X = A,V. The power 4 in (z/li)^ was based on walking technicolor arguments [22], but 
is irrelevant for LHC phenomenology: one can absorb the effect of a different power in the 
Ox value. Pure AdS geometry corresponds to ov,a — 0. Choosing boundary conditions that 
preserve only C/(l)em leads to a massless photon, and light W, Z compared to the resonances. 
Although 5D provides a tower of resonances, we restrict our study to the Ughtest two triplets of 
resonances {W^2: -^1,2)- These resonances are narrow, T ~ GeV ^ 

While we assume the strong interactions themselves to be parity symmetric ^, the coupling 
to the EW sector (set by boundary conditions) leads to physical mass eigenstates that are an 
admixture of axial and vector components. The mass splitting between resonances is directly 
affected by nonzero oy, oa, as are their couplings. Specifically, the permutation symmetry 
among triboson couplings does not hold: for B, C, D representing three different HTC spin-1 
particles (including 7, W^, Z) 

gBCD{d^,B-f+^D'^) + gBCD2{d^,C^^D'^B-^) + gBCD^{d^,Dl^B-^C+''), (1) 

we find Qbcd 7^ 93002 7^ Qbcdz- We modified both MadGraph and BRIDGE accordingly. 

In summary, the HTC description is very economical. The remaining free parameters are: 
the size of the ED (Zi), which sets the overall mass scale for the new resonances, M ~ 1/Zi, 
the amount of departure from AdS geometry {ov,a) and the coupling of the resonances to SM 
fermions (gf/v)- 

3. PARAMETER CONSTRAINTS: 

In a pure-AdS model (oy = oa = 0), consistency with precision electroweak measurements 
(especially the S parameter) requires fermiophobic resonance couplings, gjjy ~ [27]. This 
is not true in HTC, where we find regions of parameter space in which S is small due to can- 
cellations between resonance multiplets. In the lefthand side of Fig. 1, we show the line along 
which oblique corrections cancel. Along that line, the lightest two resonances are separated by 
only > 100 GcV, though the exact spacing depends on the full set of 5D parameters. The mass 
separation between the Wi and W2 greatly impacts the phenomenology, as we see in section 4. 

Having narrowed down the ov,o^ region of interest, for the present study we simply 
set the couplings of fermions to the W, Z to follow the SM relations, and take the couplings 
gffY of fermions to resonances as free parameters. The gffv are still constrained by direct 
Z' . W Tevatron cross section bounds [86,87] and by contact interaction limits [85,88]. In any 

'The tower of narrow resonances is also expected in 4D large N gauge theories. 
^Meaning 551, = gr,^, where g5L,g5R are the 5D gauge coupUngs of SU{2)l ® SU (2)ij. 
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Figure 1: Lefthand side: line of S* = in the oa,ov plane with the two sample points on the o^i = axis. 
Righthand side: setting = and Mwi — 500 GeV, value of the trilinear gauge coupling as a function of 
ov- The two horizontal Unes correspond to the Icr bounds [85]. 



particular HTC model, one must also check that the resonances do not disrupt the measured 
Tevatron diboson cross sections [89,90] and high pT,z,PT,'y distributions [89,91]. 

For a given resonance mass, the geometry parameters oy, oa are constrained by LEP 
hmits on anomalous triboson couplings [85], as depicted in the righthand side of Fig. 1. As a 
first application of the HTC framework we now summarize the LHC signals of the two HTC 
points indicated in Fig. 1. For the details of the analysis, see Ref. [92]. 

4. S-CHANNEL PRODUCTION: 

Because HTC resonances need not be fermiophobic they can be produced as s-channel reso- 
nances. Li our setup, WZ is the dominant decay mode for charged resonances ^, therefore in 
Ref. [92] we considered the mode: 

PP W^Z, W^Z ^U + v (2) 

In figure (2) we show the invariant mass distributions in the W^Z channel for the two 
sample HTC points in [92]. For these points we have set the values oi gjfv to be compatible 
with Tevatron-LEP limits and yet both resonances could be discovered within the first few fb^^ 
at the LHC. These points are just an example, chosen because they have large signals at the 
LHC. 

Qualitatively, the overall size of the signal is set by the fermion-resonance coupling a oc 
gj^y and the mass scale of the new resonances {Mw^ ~ l/^i), while the relative height of 
the peaks is determined by the relative strengths of the couplings gw^wz, gWiWZ and the mass 
separation — M^/^, both of which depend (primarily) on the geometry parameters oy, o^. 
Preliminary studies of the oy,o^ and li dependence of gwiwz^gw2wz are being carried out 
[93], however a more thorough analysis of the HTC parameter space remains to be done. 

The large s-channel signals to W Z look similar to signals of Low-Scale Technicolor 
(LSTC) [66-68]. However, in LSTC the interactions are carefully chosen to preserve an ap- 
proximate techni-parity symmetry. With this symmetry only the vector resonance couples to 
longitudinal W and Z polarizations. In HTC, we are not free to tune the interactions. All in- 
teractions, including vector-axial mixing, are determined by the 5D parameters and boundary 



•'Because our 5D setup mimics a strong sector with minimal chiral symmetry, SU (2) ^ ® SU{2) the spectrum 
doesn't contain any uneaten pseudo-Goldstone bosons (technipions). 
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Figure 2: Invariant mass distributions for Wi_2 signal and background in the WZ channel for two HTC points and 
assuming £ = 10 fb"\ In the left plot, oy = -10, = 0,^1= 6.3 TeV"\gffv = 0.1 52 (92 is the SM 
SU{2)y, gauge coupUng), while in the right plot oy — —22.5, oa — 0,li — S TeV~^,g//y — 0.05 52- The 
dominant background is SM WZ production. All backgrounds are greatly reduced by imposing cuts on the px of 
the W and Z. The tt background is suppressed further by cutting on the maximum px carried in jets. 



conditions. In the region of interest (viable with electroweak constraints) we find that techni- 
parity is not a good approximation for HTC, so both low-lying resonances couple to Wl, Zl. 
Since the resonance contribution to pp WZ is dominated by W^i 2 W^Zl , LSTC predicts 
only one peak in figure (2), while in HTC we see two. 

Another interesting s-channel production mode is 

PP Wt^2 ,W^iiy (3) 

Of the conventional three vector boson terms, the only permutation consistent with U{l)em 
gauge invariance is g'yWi,2wid[^A,y]{W[2[^Wu]) + h.c), i.e. where the derivative acts on the 
photon field. A nonzero value for only one triboson coupling permutation is not possible in tra- 
ditional, AdS-based Higgsless models. However, this final state as been considered recently [94] 
in the context of LSTC, exhibiting only one resonance. 

This channel was also investigated in [92], and in figure (3) we plot the invariant W + 'j 
mass for the same sample HTC points used in figure (2). 

As in the WZ case, the signals for both these HTC points are dramatic and could be seen 
within the first few fb \ However, the difference between HTC parameters sets is more evident 
here than in the WZ case. When oy, oa are such that the separation between the resonances 
is > 100 GeV, as in the righthand plot, only the lightest resonance is visible, whereas in the 
WZ case the second peak is still visible. The main reason for this is that the decay modes 
W2 WiZ, ZiW are open, thus suppressing the branching ratio to W'j. Since the BR to 
Wj is smaller than to WZ, only when the resonances are very degenerate, as in the lefthand 
example, are both resonances visible. 

Finally, neutral resonances Zi 2 can also be produced in the s-channel, but the most 
promising final state is into leptons rather than gauge bosons WW. Despite the smaller cross 
section, the cleaner dilepton channel may reveal both resonances within C ~ few fb^^ for the 
two points presented here. 
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Figure 3: Invariant mass for Wi^2 signal and background in the Wj channel assuming £ = 10 fb"\ The HTC 
points are the same as in Figure (2). The dominant backgrounds are SM W + j, and W + jet where the jet fakes 
a photon. The background is suppressed through hard cuts on pt^j,pt,w- 



5. VECTOR BOSON FUSION: 

Although s-channel processes will be the most important in the early years of the LHC, alter- 
native channels do exist and can expose different aspects of the resonance theory. One example 
process is vector boson fusion (VBF): pp Wi2jji ^1,2 jj- For resonances decaying to a pair 
of gauge bosons, VBF directly probes WlWl scattering and thus it is important to study regard- 
less of the fermion-resonance coupling One VBF channel, pp W^2j3 ~^ W^Zjj, was 
studied in [92]. For the two HTC points, the existence of two nearby resonances could be seen 
as two edges in the transverse mass ^ My, though only with luminosity C > O(IOO) fb^^ VBF 
signals which require the gwi,2W'y coupling, such as pp W^'^jj, may also be interesting. 

6. CONCLUSIONS: 

An effective lagrangian description of two new triplets of vector resonances would introduce 
(9(100) new parameters. In this paper we perform a first step towards an economical parametriza- 
tion of models of Dynamical EWSB: Holographic Technicolor (HTC). Although a departure 
from 5D modelling, HTC uses 5D techniques to parameterize a wide class of models in terms 
of 4 parameters: Zi, oy, o^, Qffv- After imposing current experimental constraints, we identi- 
fied the relevant region in the HTC parameter space. We have chosen two sample HTC points 
and discussed the early discovery (£ ~ 10 fb~^) of two nearby resonances in the s-channel 
pp Wi^2 WZ, W'j. The framework presented here can be extended to add new particles, 
e.g. techni-pions, techni-omegas and composite Higgs. 
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Also, in fermiophobic models where gffv ~ 0, channels such as VBF are the only way to discover the 
resonances. 

^The transverse mass is defined as: M| = {y^ M^{m) + p^m) + \E^^ss,T\f - \pt{M) + E^iss^r?- 



Dilepton Final States 



41 



Part 6 

Generic Searches for New Physics in the 
Dilepton Channel at the Large Hadron 
Collider 

G. Landsberg 

Abstract 

We propose a model-independent framework applicable to searches for 
new physics in the dilepton channel at the Large Hadron Collider. The 
feasibility of this framework has been demonstrated by the D0 searches 
for large extra dimensions. The proposed framework has a potential to 
distinguish between various types of models and determine most favor- 
able parameters within a particular model, or set limits on their values. 

1. INTRODUCTION 

This letter is devoted to searches for signals for new physics at the Large Hadron Collider (LHC) 
in the dilepton channel^ This channel has been historically fruitful for discoveries: J /psi and 
T mesons, as well as the Z boson were all discovered using dileptons. The LHC may not be an 
exception! 

The advantages of the dilepton channels for searches for new physics are numerous: 

• Easy triggering; 

• Relatively low instrumental and standard model (Drell-Yan) backgrounds; 

• Well known (NNLO) standard model (SM) cross section; 

• Number of theoretical models that predicts relatively narrow resonances with non- vanishing 
decay branching fraction to dileptons; a typical example is a generic type of models with 
an extra U{\) group, which leads to the existence of a Z' boson, often with non-zero 
couplings to dileptons. 

2. THE MODEL 

In the SM, the dilepton final state at hadron colliders is produced via the s and t-channel ex- 
change of virtual photons or Z bosons. While the s and i-channel diagrams interfere, their main 
contributions are well separated in the phase space of the dilepton system: high-p^- dileptons are 
dominantly produced via the s-channel exchange, while the t-channel process mainly results in 
very forward leptons. While this doesn't really matter for the generic formalism discussed be- 
low, we will focus on the high-p^ dileptons in the s-channel, as the most promising signature for 
new physics at the LHC, with the exception of one case - t-channel exchange of a leptoquark 
(LQ) or a SUSY particle in the models with i?-parity violation (RPV). 



In what follows by "dileptons" we will imply the dielectron and dimuon channels, and is focussed on the early 
discovery potential in this channel. The discussed formalism also applies to the ditau channel, but since this is a 
more challenging channel experimentally, we are not pursuing it here for the purpose of early searches at the LHC. 
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In the presence of additional diagrams contributing to the dilepton final state via exchange 
of new particles, the overall cross section for the dilepton production is given by the interference 
of the SM diagrams with the new ones coming from new physics. Consequently, it makes 
sense to parameterize the double-differential cross section, (P a /dMn/d cos 9*, where Mu is the 
dilepton invariant mass and cos 9* is the cosine of the scattering angle in the dilepton c.o.m. 
frame, in the following form: 

d^a _ / d^a \ / d^a \ / d^a 
dMiid cos 9* ~ \dMud cos 9* J g^^^ \dMud cos 9* ) .^^^^^ \dMud cos 9* 

Here the first term describes the SM contribution, the second term corresponds to the interfer- 
ence between the SM and new physics contributions, and the third one describes direct contri- 
bution from new physics. In terms of matrix elements, the first three terms are proportional to 
the appropriate derivatives of jTVlsMp, |-^sm-^np + -^sm-M^pI, and |A^npP> respectively. 

Eq. (1) describes well general case of new physics due to, e.g., compositeness-like oper- 
ator. However, in case new physics appears in a form of a narrow resonance, the corresponding 
matrix element has a Breit-Wigner pole, and therefore it makes sense to explicitly specify it. 
Moreover, in the case of relatively narrow resonance, the interference effect nearly cancels out 
when integrating over the width of the resonance, so it simply could be added to the above 
equation: 

_ / d^a \ / d^a \ / d^a \ 

dMiidcos9* ~ \dMiid cos 9* J gy^^ \dMiid cos 9* J .^^^ \dMiid cos 9* J 

BW(M,„M.,r.)-^, (2) 

where BW is the line-shape for a Breit-Wigner resonance with the mass Mq and width Fq, and 
Q is the angular distribution of its decay products. The remaining NP contribution described 
by the second two terms now does not include the resonance, which has been explicitly treated 
separately. In case of more than one resonance (e.g., Kaluza-Klein tower of resonances), addi- 
tional terms similar to the last one can be added. However, since e focus our attention on early 
searches for new physics, chances are that only the lowest mass resonance is going to be visible. 

Note that the advantage of using double-differential cross section for description of the 
process is that the two variables, Mu and cos 9* define the tree-level 2 — > 2 process completely. 
Thus, at leading order, they contain entire information about both the SM and new physics 
contributions, i.e. offer the most powerful separation between the NP signal and SM background 
possible in the entire phase space. 

An additional advantage of using full double-differential cross section is that the "stan- 
dard candle" - the Z peak - is contained in the data. That allows in situ calibration of the 
search sample: in particular the limits on new physics cross section, or the measurement of its 
cross section can be expressed in terms of the Z production cross section, which is well-known 
theoretically. Moreover, normalization to the Z-peak results in the reduced systematic error, as 
many important uncertainties, such as the uncertainty of the luminosity measurements, signal 
acceptance, and lepton identification efficiency would largely cancel out. 

3. THE METHOD 

The proposed method to look for generic deviations due to NP effects is based on the idea 
outlined in Ref. [95], generalized to a more complete case of non-resonant and resonant NP 
contributions. 
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The fist step is to generate templates corresponding to each of the four terms of Eq. (2) for 
a particular model. The simplest one would have a flat resonance decay angular distributions, 
i.e. dO,/d cos 9* =const, and a particular non-resonant NP model, e.g., compositeness. These 
templates are result of parton-level Monte-Carlo (MC), followed by the detector simulation, 
which includes acceptance requirements as well as smearing of particle momenta. Once the 
templates are generated, one would fit experimental cP a/ dMud cos 9* spectrum (typically in a 
form of a 2-dimensional histogram) to the sum of four terms. The results of the fit can then be 
used to evaluate the relative weights of the SM-only and SM -|- NP hypotheses and distinguish 
between various types of new physics. The MC templates also include next-to-leading order 
corrections in the form of non-zero transverse momentum of the dilepton system. Since the final 
state is colorless, the new particles contributing to the final state must be color-singlet as well. 
That implies that the only source of next-to-leading order corrections is initial-state radiation, 
which is expected to be the same for the Drell-Yan dilepton production and for contributions 
from new physics. Thus the pt spectrum of the dilepton system can be reliably modeled by 
using the well-known px spectrum of Drell-Yan pairs. 

Since for a heavy narrow resonance experimental mass resolution is typically worse than 
the internal width, F, after the detector effects the Breit-Wigner in Eq. (2) will be replaced 
with a Gaussian with the r.m.s. a determined by the experimental mass resolution. (In case the 
resonance is not sufficiently narrow, a convolution of the Breit-Wigner and a Gaussian should 
be used instead.) 

We can further quantify effects of non-resonant new physics via a parameter rj, which 
shows the relative strength of the NP contribution. For example, if the new physics has char- 
acteristics of compositeness with the scale A, a good choice for this parameter is 77 = 1/A^. 
Putting it all together, we obtain a modified Eq. (2): 



dMiid cos 9* 



fo{Mu, cos 9*)+rifi{Mu, cos 9*)Wf2(Mu, cos 9*)+N x G(Mo, a)fs(cos 9*), 

(3) 

where /o, /i, /2 and fs are the templates discussed above, G{Mo, a) is the Gaussian, and N is 
the normalization. An example of templates from Ref. [95] is shown in Fig. 1. 

The fit of the double-differential distribution in data with the function described in Eq. (3) 
yields the values and the uncertainties on the two free parameters of the fit: i] and N. Depending 
on the values of these parameters, the following broad classes of new physics can be identified 
and possibly distinguished from one the other by doing steps outlined in the last column of the 
Table: 

A simpler version of this method has been successfully used by the D0 experiment in 
searches for large extra dimensions [96, 97]. It is therefore expected that the above method 
would work well at the LHC. 

4. CONCLUSIONS 

A framework for generic searches for new physics in the dilepton channel at the LHC is dis- 
cussed. This framework would allow to statistically distinguish between various models of new 
physics and to determine parameters (or set limits on their values) within a particular model. 
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Table 1: Contributions of various types of new physics in the dilepton final state. 



New physics model 




N 


Next steps 


Standard Model 


= 


= 


Look elsewhere 


Compositeness 




= 


Look at diphotons, dijets; correlate dielectrons and dimuons. 


Large extra dimensions 


-I n 


= U 


Look at diphotons, search for black holes and monojets. 


LQ's, RPV SUSY 


7^0 


= 


Look for pair-produced leptoquarks in the dilepton + dijet 
channel; focus on a dedicated search for RPV SUSY. 


Z' 


= 


7^0 


Confirm absence in the 77 channel, determine couplings. 


Randall-Sundrum model 


= 




Confirm presence in the 77 channel, look for black holes. 


Technirho/Tecniomega 


= 




Confirm in the W + dijet channel. 




7^0 


7^0 


Correlate destructive interference with the peak height; 
look for the next excitation. 
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Figure 1: Examples of the templates for Drell-Yan production at the Tevatron: a) the SM (/□); b) /i for models 
with large extra dimensions; c) /2 for models with large extra dimensions. Here t] is related to the ultraviolet cutoff 
for the Kaluza-Klein graviton tower, Ms as r/ = 1/M|. From Ref. [95]. 
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Part? 

Rapidity and Couplings to Quarks 

TM.P. Tait 

Abstract 

The rapidity distribution of a Z' produced at the LHC encodes infor- 
mation about the relative sizes of the couplings to up quarks and down 
quarks which is different from the inclusive cross section. Thus, by 
measuring the Z' production rate at different rapidities, we can help 
pin down the coupling to up quarks independently from the coupUng to 
down quarks. 

1. INTRODUCTION 

A massive neutral vector Z' decaying into a pair of leptons would be a fascinating discovery of 
physics beyond the Standard Model, one the LHC could hope to identify with a relatively small 
sample of data provided it couples to up- or down-type quarks and decays into charged leptons, 
i'^£~ — e+e" or fx'^fx~. After the initial excitement of discovery, it will be important to measure 
as many of the Z' couplings independently as possible, to unravel the underlying model from 
which it arose. 

The primary initial observables will be the mass, given by the position of the excess in the 
lepton invariant mass distribution, and the cross section, given by the magnitude of the excess. 
Assuming the Z' has flavor universal couplings (motivated so that it does not introduce large 
FCNCs), the cross section can be parameterized at a hadron colUder by [98], 

a{pp Z' ^ £^£~X) = CuWu + CdWd (1) 

where the and Wd factors contain the parton distribution functions and higher order QCD 
corrections (exactly at NLO and to good approximation to NNLO) and thus contain only SM 
inputs, and c„ and q contain the Z'-dependent quantities, 

Cd = {gl + g^^BR{Z' ^£+r ) (2) 

where gq, gu, and gd are the Z' coupling to left-handed quarks (assumed equal by SU{2)l 
invariance), right-handed up-type quarks, and right-handed down-type quarks, and BR{Z' — > 
is the branching ratio into £'^£~. 

2. Cu VERSUS Cd 

A measurement of the cross section becomes, through, Eq. (1), a determination of a combination 
of couplings times the branching ratio into leptons. A single measurement constrains only a 
combination weighted by Wu and Wd- As a first step toward separating out c„ and Cd individually, 
we consider the cross section differential in the rapidity distribution of the Z',da/dy. We expect 
that because there are more up than down quarks in the proton, that the contribution to forward 
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Z' rapdities will depend more strongly on the coupling to up quarks than to down quarks. This 
can be captured by introducing y-dependent Wu,d{y), 

da 

3 — = CuWuivz') + CdWdivz') (3) 
dyz' 

where the integral of Wu{d){y) over y reproduces the original Wu(d) in Eq. (1) for the inclusive 
Z' cross section. 

To illustrate how this works, in Figure 1, we present the ratio of Wu{y) /wd{y) as a function 
of y for a Z' of mass 2 TeV, computed at tree level. A more precise analysis would want to use 
higher order corrections, which are known to NNLO for this distribution [99], and are straight- 
forward to include. For this study which aims to examine a proof of principle, tree level is 
sufficient to estimate the utility of the measurement. (This quantity is equivalent to the ratio of 
cross sections da jdy for two Z's of equal masses and couplings, one coupling only to u-quarks 
and the other coupling only to d-quarks). We see that at high rapidities w^i^y) ^ Wd{y), and the 
sensitivity to the up-type quark couplings is enhanced. 




Figure 1: Left: The ratio of Wu{y) /wd{y) as a function of y. Right: The uncertainties on the combination of c„ 
and Cd obtained from inclusive (red) and differential (blue) Z' cross section measurements, for the Z' considered 
in the text. 



As an example to illustrate how this could work, we consider an example Z' with Mz' = 
2 TeV, gz' = 0.1, for which all quarks and leptons have the same charge (taken to be unity). 
This Z' has Cn = Cii — 3.2x 10 ^ and is compatible with LEP-II data. We imagine two 
measurements - one of the inclusive Z' — > e+e^,/i+/i^ and another with a cut of \yz'\ > 1, 
selecting only Z's in the forward region. No detector effects or efficiencies are included. We 
apply statistical errors assuming 100 fb~^ of collected data. These turn out to be on the order 
of 1% (3%) for the inclusive (cut) Z' cross sections we consider, and thus are probably similar 
in magnitude to the the residual theoretical error in an NNLO calculation, and coming from the 
PDFs. 

We unfold the la uncertainties in the combination of and for both measurements 
and plot the results in Figure 1 . Also shown on the plot is the correct value of Cu and for the 
example Z'. The plot illustrates that adding the measurement of the cut Z' cross section does 
help break the degeneracy in and Cd which is left after measuring the inclusive cross section. 
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3. CONCLUSIONS 

We examine the possibility that one can use the Z' rapidity distribution to help disentangle the 
coupling of the Z' to up- and down-type quarks. The result is encouraging, and would motivate 
further investigation, including realistic detector simulations, a better treatment of uncertainties, 
and a more sophisticated analysis such as comparison with the full Z' rapidity distribution (and 
not just a cut rate measurement). 

Between the time this study was performed at Les Houches 07 and the writing of this 
report, a more comprehensive theoretical treatment of Z' observables at the LHC has appeared 
[100], including discussion of this (and other) observables. The results of that study seem 
to continue to suggest that the Z' rapidity distribution is a useful quantity to unravel the Z' 
couplings to quarks. 
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Parts 

A search for top partners at the LHC 
using same-sign dilepton final states 

T. Bose, R. Contino, M. Narain and G. Servant 
Abstract 

A natural, non-supersymmetric solution to the hierarchy problem gener- 
ically requires fermionic partners of the top quark with masses not much 
heavier than 500 GeV. We study the pair production and detection at 
the LHC of the top partners with electric charge Qe = 5/3 (T5/3) and 
Qe = —1/3 (B), that are predicted in models where the Higgs is a 
pseudo-Goldstone boson. Both kinds of new fermions decay to Wt, 
leading to a ttWW final state. We focus on the golden channel with 
two same- sign leptons, that offers the best chances of discovery in the 
very early phase of LHC and permits a full mass reconstruction of the 
T5/3. Samples are processed with the CMS Fast Simulation. 

1. INTRODUCTION 

The most notorious example of symmetry protection for the light Higgs is Supersymmetry: 
according to its paradigm, the radiative correction of each SM field to the Higgs mass is fine 
tuned against that of a superpartner of opposite statistics. The top quark contribution, in partic- 
ular, is balanced by the contribution of its scalar partners, the stops. Another kind of symmetry 
protection, however, could be at work: the light Higgs could be the pseudo-Goldstone boson 
of a spontaneously broken global symmetry [101-103]. In this case the radiative correction of 
the top quark to the Higgs mass is balanced by the contribution of new partners of the same 
spin. The naturalness criterium suggests that these new heavy fermions should have masses 
below, or not much heavier than, 1 TeV. The production and detection of these top partners at 
the LHC was studied in [104] focussing on final states with two same-sign leptons. This note 
reports those findings and complements them with a preliminary analysis using the CMS Fast 
Simulation [105]. 

Particularly motivated is the possibility that the spontaneous breaking of the global sym- 
metry and the new states originate from a strongly-coupled dynamics. This would allow for a 
complete resolution of the Hierarchy Problem without the need of fundamental scalar fields, and 
would make it possible to generate a large enough quartic coupling for the Higgs via radiative 
effects. 

The LEP precision data are crucial in guiding our theoretical investigation, as they seem 
to be compatible only with a specific kind of strong dynamics: the new sector must possess a 
custodial symmetry Gc =SU(2)c to avoid large tree-level corrections to the p parameter [106]. 
This in turn implies an unbroken SU(2)ixSU(2)/jxU(l)x invariance of the strong dynamics 
before EWSB, meaning that its resonances, in particular the heavy partners of the top quark, will 
fill multiplets of such symmetry. It has been recently pointed out [107] that the LEP constraint 
on the ZhihL coupling is more easily satisfied if the custodial symmetry of the strong sector 
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includes a LR parity, Gc =SU(2)c x Plr, and couples linearly to a composite fermionic 
operator transforming as a (2, 2)2/3 under SU(2) l x SU(2)/j x U(1)x (hypercharge being defined 
as F = T]^ + X). In this case, the heavy partners of (t^, 6^) can themselves fill a (2, 2)2/3 rep- 
resentation. The latter consists of two SU(2)i doublets: the first, (T,B), has the quantum 
numbers of (t^, 6^); the second - its "custodian" - is made of one fermion with exotic electric 
charge Qe = +5/3, T5/3, and one with charge Qe = +2/3, T2/3. Since the Higgs trans- 
forms like a (2, 2)o, the partners of tji, if any, will form a (1, 1)2/3 or a [(1, 3) © (3, l)]2/3 of 
SU(2)ixSU(2)^xU(l)x [107]. 

As discussed in more details in Ref. [104], these heavy partners couple strongly to the 
third generation SM quarks plus one longitudinal W, Z gauge boson or the Higgs. As in [104], 
we focus on the pair production of the B and of its custodial partner T5/3. Once pair produced, 
both the heavy bottom B and the exotic T5/3 decay exclusively to W^W^W''W^bb, although 
with different spatial configurations as dictated by their different electric charges, see Fig. 1. 
In the case of the T5/3 the same-sign W's come from the decay of the same heavy fermion, 
while in the case of the heavy bottom they come from different B's. Selecting events with two 
same-sign leptons thus allows one to fully reconstruct the hadronically-decaying T5/3. Even 
though a full reconstruction of the B is not possible, it was found in [104] that the same-sign 
dilepton channel is probably the most promising one for its discovery (see [108, 109] for other 
recent studies of pair-production of 5-type quarks). 

In this work, we report some of the results derived in Ref. [104] and we improve on 
them by adding the CMS Fast Simulation [105] to the analysis. Section 2 presents our Monte 
Carlo simulations. In Section 3, we discuss our cuts and derive the discovery potential, while in 
Section 4 we reconstruct the W and t candidates and pair them to reconstruct the T5/3 invariant 
mass. 




Figure 1: Pair production of T5/3 and B to same-sign dilepton final states. 



2. FAST SIMULATION AND EVENT SELECTION 

We focus on the pair production at the LHC of B and T5/3 to two same-sign leptons : 

gg, qq^BB, T^/^n,^ ^ tiW^W- l^ul^ubbqq'qq' . (1) 



51 



The physical, observed final state is of the fomi 

pp — > l^l^ + n jets + , ^ = e, // , (2) 

where the number of jets depends on the adopted jet algorithm and on its parameters. We take 
M — 500 GeV for the mass of the heavy fermions (we refer the reader to [104] for the related 
study of the M = 1 TeV case) and set the coupling to tW to be Xt^^^ — Xb = 3. Such 
large values of the couplings are naturally expected if the heavy fermions are bound states of a 
strongly coupled sector, and is mainly composite. ^ 

The most important SM backgrounds to the process of eq.(2) are tiW+jets, tiWW+jets 
(including the tth + jets resonant contribution for > 2mw), WWW + jets (including the 
Wh+jets resonant contribution for nih > ^niw), W'^W'^ +jets and Wl~^l~ +jets (including 
the WZ + jets contribution) where one lepton is missed. To be conservative and consider 
the case in which the background is largest, we have set the Higgs mass to rrih = 180 GeV. 
This greatly enhances the tiWW and WWW backgrounds. The production cross sections for 
the signal and for the various backgrounds are reported in Table 1. No K-factors have been 
included, since those for the backgrounds are not all available (the K-factor for the signal is 
~ 1.8 for M = 500 GeV [110]). 





a[fb] 


a X BRilH^) [fb] 


T5/3T5/3/55 + jets (M = 500 GeV) 


2.5 X 10^ 


104 


T,/sT,/s/BB + jets (M = 1 TeV) 


37 


1.6 


tiW+W' + jets tth + jets) 


121 


5.1 


ttW^ + jets 


595 


18.4 


W+W-W^ + jets hW^ + jets) 


603 


18.7 


W^W^ + jets 


340 


15.5 



Table 1: Signal and background cross sections at leading order (left column). The right column reports 
the cross section times the branching ratio to two same-sign leptons final states (e or ji). 

Given its complexity, we were not able to fully simulate the Wl^l" + jets background, 
and for that reason we have not included it in our analysis. However, we estimated that after 
the main cuts presented below, the Wl'^l~ + jets background is expected to be smaller than 
~ 30% of the sum of the other backgrounds [104]. Even though this is not entirely negligible, 
the error due to its exclusion is within the uncertainty of our leading-order analysis. More- 
over, the Wl^l^ + jets cross section is expected to be strongly suppressed after requiring the 
reconstruction of one W and one top as done in the last section. 

Another potential source of background are ti+jets events where the charge of one of the 
two leptons from the top decays is misidentified. Given the large tt + jets cross section, even a 
charge misidentification probability s„ds ~ a few x lO^"^ would result into a same-sign dilepton 
background of the same order of ttW + jets. ^ The hardest lepton in the ti + jets events has 

'Notice, however, that our final results are largely independent of the specific values of AT5/3, A^, since the 
latter determine only the decay width of the heavy fermions. For our choice of couplings and M = 500 GeV, 
r = 31 GeV. 

^Requiring the reconstruction of one W and one top as in Section 3 is however expected to reduce significantly 
more the ti + jets events background than tiW + jets or tiWW + jets. 
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a pt distribution peaked at values smaller than 100 GeV (the second hardest lepton has instead 
a significantly softer pr distribution). From the latest ATLAS and CMS TDRs, probabilities as 
low as ~ 10~^ seem to be realistic in the case of 100 GeV muons, while slightly larger values 
are expected for electrons [8, 111]. In such case the tt + jets background would be safely 
negligible. In the absence of a realistic estimate of e„iis as a function of the lepton's px and 
pseudo-rapidity, we decided not to include the ti + jets background events in our analysis. It 
is however clear that a specific and accurate estimate of this background is required to validate 
our results. 

In these proceedings, we review part of the findings obtained with the simulation and 
analysis of Ref. [104], that we will refer to as generic analysis, and compare with the results 
obtained by adding the CMS Fast simulation of the detector. Both the signal and the SM back- 
ground events were generated at the partonic level with MadGraph/MadEvent [35,76, 112], ^. 

In the generic analysis, we have used Pythia [75] for showering and to include the initial 
and final-state radiation. Hadronization and underlying event have been switched off for sim- 
plicity. Jets have been reconstructed using the GetFet [113] cone algorithm with = 30 GeV 
and a cone size Ai? = 0.4. The parton-jet matching has been performed following the MLM 
prescription [1 14, 1 15]. No detector effects were taken into account, except for a simple gaus- 
sian smearing on the jets. Both the jet energy and momentum absolute value were smeared 
by ^E/E — 100%/>/E7GeV, and the jet momentum direction using an angle resolution 
A0 = 0.05 radians and A77 = 0.04. 

On the other hand, in the CMS Fast Simulation analysis [105], the samples are first pro- 
cessed via PYTHIA for showering, to include initial and final-state radiation and to fragment 
and hadronize quarks and gluons. In addition, the underlying event is switched on. A jet- 
matching algorithm, following the MLM prescription [116] is employed to ensure that there 
is no double counting due to the parton showering in PYTHIA. These samples are then pro- 
cessed via a dedicated fast simulation processor (FAMOS) of the CMS detector and the event 
reconstruction software. Jets with cone size Ai? = 0.5 are reconstructed with the iterative cone 
algorithm. Generic jet energy corrections are applied to the cone jets. The jet resolution used 
in the CMS Fast Simulation is AEt/Et = 1.25/^Er/GeV © 5.6/Er/GeV © 0.033 [8]. This 
jet energy resolution will be the main source of different results with the generic analysis. 

A first important information on the kinematics of signal and background events comes 
from the number of reconstructed jets. Fig. 2 compares the distributions for the number of jets 
in both analysis. The background and signal values are peaked at 2 and 5 (3 and 6) respectively 
in the generic (CMS Fast Simulation) analysis [105]. By signal here we mean either T^/sT^/s 
or BB events. The CMS Fast Simulation leads to a slightly larger number of reconstruted jets 
above pt > 30 GeV due to a combination of the underlying event settings, jet reconstruction 
algorithm, and energy corrections used. The total background distribution is peaked at smaller 
values. This is mainly due to the low jet multiplicity in the WWW + jets and VT^'^H^^'' -|- jets 
backgrounds. In the case of the signal, the hard scattering process produces 6 quarks, after 
the decay of the top and of the W. It turns out that for M = 500 GeV the 5-jet bin is mostly 
populated by events where the 6th jet is lost because it is too soft (i.e. it does not meet the 
minimum transverse energy requirement, Et > 30 GeV). 

^The factorization and renormalization scales have been chosen as follows: /j. = Mt,b for the signal; ji = 
2mt + mw for ttW + jets; fi = 2mt + mh for ttWW + jets; = mw + mh for WWW + jets; /x = 2mw 
for W^W^ + jets. 
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Figure 2: Fractions of signal and background events with a given number of jets for iS™" = 30 GeV. Left plot is 
from Ref [104] (using the GetJet cone algorithm with AR — 0.4) and right plot is from the CMS Fast Simulation 
(using iterative cone algorithm with AR = 0.5). 



3. DISCOVERY POTENTIAL 

In this section, we focus on the discovery of the top partners, proposing a simple strategy that 
does not rely on any sophisticated reconstruction, nor does it require 6-tagging. To isolate the 
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Figure 3: pr distributions of the first and fifth hardest jets (left) and of the first and second hardest leptons (right) 
before any cut (from the generic simulation). 
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signal, the following cuts are applied (both for the generic and the CMS Fast simulations): 

{PT(lst) > 50 GeV ( PT(lst) > 100 GeV 

PT(2nd) > 25 GeV jets : I pri^nd) > 80 GeV > 20 GeV , 

\r]i\ < 2.4 , ARi, > 0.4 [rijet > 5, < 5 

(3) 

where 1 st and 2nd refer respectively to the first and second hardest jet or lepton (electron or 
muon). These cuts were motivated by the pt distributions shown in Fig. 3. The relative effi- 
ciencies are reported in Table 2. 





signal 
(M = 500 GeV) 


ttw 


ttww 


WWW 


W^W^ 


Efficiencies (e^am) 

Cr [fb] X BRX Emain 


0.42 (0.55) 
44.2 (57.2) 


0.074 (0.15) 
1.4 (2.80) 


0.12 (0.06) 
0.62 (0.30) 


0.008 (0.02) 
0.15 (0.37) 


0.01 (0.03) 
0.16 (0.48) 



Table 2: Efficiencies of the main cuts of eq.(3) for the generic analysis. Here signal means either Tc^/^T^/^ 
or BB events. Numbers in parenthesis refer to the CMS Fast Simulation. [105] 



After the cuts, the signal is much larger than the background. A rough indication on the 
mass of the heavy fermions can be extracted from the distribution for the scalar sum of px of all 
jets in the event {Ht), see Fig. 4. This is a quantity that is not much affected by a poor jet energy 
resolution, and as such it is particularly appropriate in the first LHC phase. Further evidence 
on the production of a pair of heavy fermions comes from the distributions of the total invariant 
mass, the invariant mass of the hardest 5 jets, and transverse mass of the {llufj) system. We 
refer the interested reader to Ref. [104] for the relative plots. 



CMS Preliminary 1 fb"' 




Figure 4: Ht distribution (scalar sum of pT of all jets in the event). Comparison between the simulation of 
Ref. [104] (left) and the output of the CMS Fast Simulation for an integrated luminosity corresponding to 1 fb^^ 
(right). 

By counting the number of signal and background events that pass the main cuts of eq.(3), 
one can estimate the statistical significance of the signal over the background, as well as the 
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minimum integrated luminosity required for a discovery. In the generic analysis, we defined the 
latter to be the integrated luminosity for which a goodness-of-fit test of the SM-only hypothesis 
with Poisson distribution gives a p-value = 2.85 x 10"^, (see for instance [85]). ^ Results are 
reported in Table 3. The numbers in parenthesis are those from the CMS Fast simulation, were 
we used a different and more conservative estimate for the significance: S/ \/S + 2B. In the 







S 


B 


Ldisc 


M = 500 GeV 


B only 


864 (1124) 
424 (548) 


23 (42) 
23 (42) 


56 pb-' (240pb^i) 
147 pb-^ (530 pb-^) 



Table 3: Number of signal (S) and background (B) events that pass the main cuts of eq.(3) with L = 
10 fb^^ for M = 500 GeV. The last column reports the corresponding integrated luminosity needed 
for the discovery (Ldisc), as computed in the generic analysis by means of a goodness-of-fit test with 
Poisson distribution and p-value = 2.85 x 10~^. Values in parenthesis are from the CMS Fast Simulation 
and in the last column L^isc is computed using S/\/S + 2B as an estimate of the significance. 

most favorable case where both T5/3 and B partners exist and have mass M = 500 GeV, a 
discovery will need less than ~ 250 pb~^. Our estimates should be conservative, as we did not 
include any K-factor in our analysis, although it is known that next-to-leading order corrections 
enhance the signal cross section by ~ 80% for M — 500 GeV [110]. Even a common K-factor 
K, for both the signal and the background would imply a statistical significance larger by a factor 
~ a/k,, as well as a discovery luminosity smaller by the same factor. 

4. T5/3 MASS RECONSTRUCTION 

More direct evidence for the production of a pair of T5/3 or B comes from reconstructing the 
hadronically decayed top quark and W boson, as well as from the distribution of the invariant 
mass of their system. The mass reconstruction of T5/3 crucially depends on the jet energy 
resolution when using the full detector simulation. We find that the width of the reconstructed 
top quark mass as obtained in the generic analysis is 12.5 GeV, while that obtained from the 
CMS Fast Simulation is about 19.5 GeV. This latter value more realistically reproduces the 
expected resolution of in the early phase of the LHC. It will improve as the flavor dependent 
jet energy corrections and use of high purity fe-tagging for correct assignment of the 6-jets are 
incorporated in the analysis. 

We report here the results from the generic analysis using the energy resolution defined 
in Section 2, which is probably optimistic for the early LHC phase. We first select the events 
where two W's can be simultaneously reconstructed, each W candidate being formed by a pair 
of jets with invariant mass in the window \M{jj) — mw\ < 20 GeV. To avoid wrong pairings 
and reduce the fake ones from the background, we impose the following cuts: 

Ai?j-j- < 1.5 , \p{W) I > 100 GeV on the first W candidate ; (4) 

/\Rjj < 2.0 , \p{W) I > 30 GeV on the second W candidate . (5) 

The pt cuts, in particular, have been optimized using the pr distributions of the W and top 
from the signal events (see [104]). If more than one pair of W candidates exists which satisfies 



This p-value corresponds to a 5ct significance in the limit of a gaussian distribution. 
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the above cuts, we select that with the smallest = Ai?|^( 1st pair) + Ai?|^ (2nd pair). We 
then reconstruct the top by forming Wj pairs, made of one W and one of the remaining jets, 
with invariant mass in the window \M{Wj) — mt\ < 25 GeV. If more than one top candidate 
exists, we select that with invariant mass closest to m^. We discard events where no top can 
be reconstructed. The efficiencies of this reconstruction algorithm are reported in Table 4. The 





signal 
(M = 500 GeV) 


ttw 


ttww 


WWW 


WW 




0.62 


0.36 


0.49 


0.29 


0.15 


^top 


0.65 


0.56 


0.64 


0.35 


0.35 



Table 4: Efficiencies for the reconstruction of two W's {e2w) and one top (etop) using the algorithm and 
the cuts described in the text for the case M = 500 GeV. 

distribution of the Wt invariant mass is plotted in Fig. 5. As expected, in the scenario with 
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Figure 5: Invariant mass of the Wt system for M = 500 GeV with L = 10fb~^. The dotted curve refers to the 
case in which 6-tagging is performed in the reconstruction. It assumes two b tags, though no 6-tagging efficiency 
has been included. 

T5/3 partners there is a resonant peak centered at Mt^^^ = 500 GeV, while the distribution has 
a non-resonant, continuous shape if only Qe = —1/3 heavy fermions exist. The dotted curve 
refers to the case in which 6-tagging is performed in the reconstruction algorithm. More in 
detail, we have selected events with two b tags and we have reconstructed the top from Wb 
pairs, excluding at the same time the b jets when selecting the W jet pair candidates. No b- 
tagging efficiency has been included, in order not to commit to any specific value (hence, for a 
6- tagging efficiency the final distribution will be rescaled by a factor el). 

5. CONCLUSIONS 

Our results show that the analysis of final states with two same-sign leptons at the LHC is an 
extremely promising method to discover the top partners B and T5/3. By requiring two same- 
sign leptons one avoids the large ti background and selects a particularly clean channel where 
evidence for the existence of the heavy fermions could come in the early phase of the LHC. 
If both B and T5/3 exist, a discovery will require only ~ 60 - 250 pb^^, depending on how 



M = 500 GeV L=10fb-i 
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the statistical significance is estimated. Even without 6-tagging, and before reconstructing the 
hadronically decayed W and top, one can have a first crucial indication on the value of the mass 
of the heavy fermions from the distributions of the total invariant mass, the invariant mass of 
the hardest 5 jets or the Ht distribution (scalar sum of jet pr's). 

Ultimately, a crucial information to understand the origin and the role of the heavy fermions 
would come from the measurement of their decay width, which will in turn lead to a determina- 
tion of their couplings Xt^^^^b- A large value of Xt^/^,b will be strong circumstantial evidence 
for the compositeness of the heavy fermions. Extracting the decay width from the invariant 
mass distribution will be challenging, as one will have to cope with the issue of jet energy res- 
olution. Most likely, a measurement will be possible only with large statistics and will require 
sophisticated W and t reconstruction techniques. 

Given the strong theoretical motivations for a search of the heavy partners of the top, we 
think that our explorative study would also deserve to be followed by a more detailed investiga- 
tion. Our results suggest that the same-sign dilepton channel might be one of the golden modes 
to discover the top partners B and T5/3, but only a complete analysis with a full simulation of the 
detector effects, an exact calculation of the Wl^l~ +jets and ti+jets backgrounds, and the use 
of fully realistic reconstruction techniques will eventually establish its ultimate potentialities. 
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Part 9 

Top Quark Pairs as a Model-Independent 
Discriminator of New Physics at the LHC 

D.G.E. Walker 

Abstract 

We review top quark pair production as a way to probe new physics 
at the LHC. Our scheme requires identifying integer-spin resonances 
from tt semileptonic decays in order to favor/disfavor new models of 
electroweak physics. The spin of each resonance can be determined by 
the angular distribution of top quarks in their cm. frame. In addition, 
forward-backward asymmetry and CP-odd variables can be constructed 
to further distinguish the new physics. We parametrize the new reso- 
nances with a few generic parameters and show high invariant mass top 
pair production may provide a framework to distinguish models of new 
physics beyond the Standard Model. 

1. INTRODUCTION 

A considerable number of experimentally viable extensions to the Standard Model (SM) have 
been proposed to describe electroweak symmetry breaking. A sample of popular scenarios in- 
clude: the Minimal Super symmetric Standard Model (MSSM) [117], models with new strong 
dynamics [44,49, 118-122], composite Higgs models at the TeV scale [123], Little Higgs mod- 
els [124-126], and models with extra dimensions at the electroweak scale [127, 128]. In addi- 
tion, some string-inspired extensions [129, 130] can also lead to new signatures. In this review, 
we describe a scheme to sort out the particle content in a model independent fashion. 

The LHC will be a "top factory": About 80 million tt events will be produced via QCD 
production for integrated luminosity of 100 fb~^. With such a large number of events (and 
knowing natural models have a preferential coupling to top quarks), it seems worthwhile to 
study top pairs. If the new physics contributes to ti production as an s-channel resonance, we 
want to identify the signal as a bump [131] on the smoothly falling tt invariant mass distribu- 
tion. We want to then reconstruct the tt cm. frame so the integer spin of the resonance can 
be determined from the polar angular distribution of the top quark. Further, any asymmetry of 
this distribution could probe possible chiral couplings. In addition, any CP properties of the 
couplings could be elucidated with the help of CP-odd kinematical variables constructed from 
the final state particle momenta. 

We look at only the tt semileptonic decay mode: tt bjij2 b£~P + c.c. where £ = e or 
The purely hadronic decay mode not only suffers from a much larger QCD background, but 
also loses the identification of t from t. For the purely leptonic mode, with a small branching 
fraction of about 4/81, one cannot reconstruct the it invariant mass with two missing neutrinos. 
Thus, our signal will be an isolated charged lepton plus missing energy (fr), 2 6-jets plus 2 
light jets. Note: The branching ratio of the semileptonic to the hadronic channel is 2/3. 

In this article, we summarize the analysis of [132]. Since [132] appeared, many papers 
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have been published which look at the effectiveness of using top quarks to probe new physics. 
A brief sample of papers is [133-142]. Please note: [132] focused on early indications of new 
physics by restricting to only 10 fb~^ of data. When necessary, for emphasis, we include the 
relevant analysis for 100 fb~^. In all, the exposition of this review follows [132] closely. 

2. EVENT RECONSTRUCTION 

In order to search for new physics, we first have to reliably reconstruction ti semileptonic decays 
at high invariant mass on an event-by-event basis, ti semileptonic decays contain a neutrino in 
the final state - thereby complicating the reconstruction. The transverse momentum of the 
neutrino is identified with the observed The neutrino longitudinal momentum is subject to 
a two-fold ambiguity from solving the kinematic quadratic equation. 

We choose a reconstruction method useful for high invariant mass ti events [132]. Several 
top reconstruction methods have been used at the Tevatron [143, 144]. There, however, the top 
quarks are produced near threshold and the kinematics of the subsequent decay products are 
very complicated. High invariant mass ti events tremendously simplify the kinematics, espe- 
cially by distinguishing the b quark from the b. Throughout this review, we use a 2 — > 6 partonic 
level monte-carlo simulation that incorporates full spin correlations from production through 
decay [145, 146]. [132] made a Pythia simulation, including gluon radiation and hadronization, 
that confirmed the results presented here. 

To ensure only high invariant mass ti events, we impose a cluster transverse mass cut on 
the ti system 

Mt = ^JiPb+Pb+Ph +Ph +Pt? + + It > 600 GeV. 

Following ATLAS and CMS, we adopt the following kinematical cuts: 

E{, > 20 GeV \r]j\ < 2.5 
> 20 GeV \r]t\ < 2.5. 

In addition, we require a lepton-jet and jet-jet isolation cut and a ^ cut of 

Ai? > 0.4 $r> 20 GeV. 

The hadronic energy is smeared according to a Gaussian error given by ^Ej/Ej = 0.5/ y^Ej/GeV © 
0.03. Additionally, the lepton momentum is smeared by Ap|./p^ = 0.36(]9f./TeV)©0.013/Vsin^. 
Here 9 is the polar angle of the lepton with respect to the beam direction in the lab frame. 

We choose a reconstruction scheme that takes Mw and as inputs for their on-shell 
production and decays. We break the scheme into the following steps: 

Step I: Demand mf^ = M^. The longitudinal momentum of the neutrino is formally expressed 
as 

where A = + 2peT- $r- If A? - Ap^j, ^ > 0, the value of p^L that best yields the known 
top mass via mf^^ = is selected. This ideal situation may not always hold when taking into 
account the detector resolutions. For cases with no real solutions, we then proceed to the next 
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Figure 1: The W and top mass reconstructions from our reconstruction scheme. Step I (Step 11) is solid (dashed), 
respectively. 



Step. 

Step II: To better recover the correct kinematics, we instead first reconstruct the top quark 
directly by demanding mf^,j — m^. The longitudinal momentum of the neutrino is expressed as 

P.L = A'pul/2{EIi-pIil) ± 



X (pLA'' + {El - pI,) {X^ - AEim"\ 

where A = — M^i + 2phiT ■ The two-fold ambiguity is broken by choosing the value 
that best reconstructs = mf^. A plot of the top and W mass distributions is shown in Fig. 1. 
The solid histogram is from the procedure Step I, and the dashed histogram from Step II. With 
these two steps, there could still be some events that do not lead to a real solution. Because 
we want accurate reconstruction of ti at high invariant mass, we thus discard them in our event 
collection. The discard rate is about 16%. 



3. BACKGROUNDS 

The major backgrounds to our ti events include the processes W+ jets, Z+ jets, WW, W Z and 
ZZ. Both the ATLAS and CMS Technical Design Reports [37, 147] performed detailed studies 
of the selection efficiencies for these background processes in comparison to a reconstructed ti 
semileptonic signal. The ATLAS (CMS) group found for an integrated luminosity of 10 fb~^ 
(1 fb^^) a signal to background ratio of S/ B = 65 (S/B = 26). This ratio has been obtained 
using the kinematical cuts listed in the previous section. Because of the expected high S/B 
ratio, our analysis is concentrated solely on the events without including the small background 
contamination. Our analysis does not include misidentification of faked leptons from jets in ti 
total hadronic decays. 

Although the t (i) is primarily identified by the charged lepton, (f-^), a concern is the 
matching of the b-jet associated with this top quark decay. Both ATLAS and CMS studies [37, 
147] show a combination of kinematic fits, designed to properly reconstruct the W boson and 
the hadronically decaying top significantly reduces misidentification. The cut on Mt helps 
significantly in this regard. 
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4. SEARCHING FOR NEW PHYSICS 




Figure 2: Reconstructed tt invariant mass distributions. The plot features a 1 TeV resonance with a total width of 
2% (soUd), 5% (dashed), and 20% (dotted) of the resonance's mass. 

We want to search for new resonances that couple to tt in a model independent manner. We 
consider ti production via 

gg^ (/)^tt, qq^ V ^ tt, qq, gg ^ h ^ tt, 

where 0, V and h are the spin-0, spin-1, and spin-2 resonances. We characterize the effects on 
the invariant mass spectrum with three parameters: mass, total width, and the signal cross sec- 
tion normalization (o;^). The normahzation uj — 1 defines our benchmark for the spin 0, 1 and 
2 resonances. They correspond to the SM-like Higgs boson, a Z' with electroweak coupling 
strength and left (L) or right (R) chiral couplings to SM fermions, and the Randall- Sundrum 
graviton h with the couplings scaled as for hqq, and (Aln(Mp^/A))~^ for hgg, respec- 
tively.^ Numerically, we take A = 2 TeV. 

In Fig. 2 we show the reconstructed tt invariant mass distribution for our reconstruction 
scheme. The SM tt total cross section is theoretically known beyond the leading order in QCD 
[150-152]. We thus expect to have a good control of this distribution even at high invariant 
masses. As for new physics, we include the contribution of a 1 TeV vector resonance for 
illustration, for uj^ = 1, with total widths specified in the caption of Fig 2. 

We maximize the signal observability by isolating the resonance within an invariant mass 
window of ±100 GeV, ±30 GeV and ±25 GeV for the scalar, vector and graviton resonance, 
respectively. Given a resonance mass and total width, we can quantify how large u needs to be 

for a 5a discovery. With the number of events for a signal (S) and background (B), we require 
S/^/B + S > 5. This translates to a bound cj^ > (25 + 5^25 + AB)/2Si where is the 
benchmark signal rate for a; = 1. This is illustrated by Fig. 3 versus the mass for a scalar. 



^More precisely, we use the Feynman rules given in [148] and include the additional warp correction factors 
from [149]. 
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Figure 3: Normalization factor versus the resonance mass for the scalar (dashed) with a width-mass ratio of 20%, 
vector (dot-dashed) with 5%, and graviton (solid) 2%, respectively. The region above each curve represents values 
of w that give 5cr or greater statistical significance with 10 fb~^ (left panel) and 100 fb~^ (right panel) integrated 
luminosity. 

vector and graviton resonance for total widths of 20%, 5%, and 2% of its mass, respectively, for 
an integrated luminosity of 10 fb~^ (Fig. 3, left panel) and 100 fb~^ (Fig 3, right panel). 

It is of critical importance to reconstruct the cm. frame of the resonant particle, where 
the fundamental properties of the particle can be best studied. In Fig. 4, we show the top 
quark angular distribution, cos 6'*, with d* defined as the angle in the tt cm. frame between 
the top-quark momentum and the incident quark momentum, with the latter determined by the 
longitudinal boost direction of the cm. system. Although events in the forward and backward 
regions are suppressed due to the stringent kinematical cuts, we still see the impressive features 
of the rf-function distributions^ in Fig. 4: a flat distribution for a scalar resonance (dashed), rfj^^ 
distribution for the left/right chiral couplings of a vector (dotted), and from qq (solid) and 
c?2±i from gg (dot-dashed) for a spin-2 resonance. To illustrate the statistical sensitivity for 
observing the characteristic distributions, we show in Fig. 5 the expected SM ti events (solid) 
with Icr statistical error bars in each bin for a 10 fb^^ and 100 fb~^ integrated luminosity, along 
with a 5(j signal of a chirally coupled vector summed with the tt background in the resonant 
region (dashed). Due to the large event sample, the statistical significance is evident in the 
central and forward region. The forward-backward asymmetry in cos Q* can thus be constructed 
to probe the chiral couplings of the particle to the top quark. With the identification of the 
charged leptons, one may even form kinematical triple products to test the CP properties of the 
couplings [153, 154]. 

5. CONCLUSIONS 

In summary, we reviewed the use for top quarks in discovering new physics in the form of 
integer-spin resonances. We showed the use of angular distributions of the top in the recon- 
structed CM frame to reveal the spin of the resonance, and the relative contribution from the 



^For the definition and convention of the rf-functions, we follow the PDG. [85] 
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Figure 4: Polar angular distributions for the top quark in the cm. frame, (a) Signal only by the for a scalar (dashed), 
a vector (dots), and a graviton from qq (sohd) or from gg (dot-dashed). 
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Figure 5: Polar angular distributions for the top quark in the cm. frame. Number of events for the SM tt back- 
ground (solid) with 1(7 statistical error bars, and the backgroimd plus a vector resonance (dashed) for (a) 10 fb~^ 
and(b) 100 fb"! 



initial states qq or gg. The forward-backward asymmetry and CP-odd variables can be con- 
structed to further differentiate models. Because SM top quark physics is well predicted, high 
invariant mass top pair production may provide an early indicator for new physics beyond the 
Standard Model at the LHC. 
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Part 10 

Production of Kaluza-Klein excitations of 
gauge bosons at the LHC 

A. Djouadi, G. Moreau andR.K. Singh 
Abstract 

We consider the Randall-Sundrum model with fields propagating in the 

bulk and study the production of the strongly and weakly interacting 
gauge boson Kaluza-Klein excitations at the LHC. These states have 
masses of order of a few TeV and can dominantly decay into top quark 
pairs. We perform a Monte Carlo study of the production process pp — > 
tt in which the Kaluza-Klein excitations are exchanged and find that 
the latter can lead to a significant excess of events with respect to the 
Standard Model prediction. 

1. INTRODUCTION 

If the extra-dimensional model suggested by Randall and Sundrum (RS) [128, 155] is to solve 
the gauge hierarchy problem of the Standard Model (SM), the masses of the first KK exci- 
tations of gauge bosons, M^k, must be in the vicinity of the TeV scale. Direct experimen- 
tal searches for KK gluon excitations at the Tevatron lead to Mkk ^ 800 GeV [156, 157], 
while high-precision measurements impose stronger bounds as the exchanges of the KK ex- 
citations lead to unacceptably large contributions to the electroweak observables [158-163]. 
Nevertheless, it was shown [16] that if the SM symmetry is enhanced to the left-right structure 
SU(2)l X SU(2)r X U(1)b-l, with B and L the baryon and lepton numbers, electroweak preci- 
sion data can be fitted while keeping the KK masses down to the acceptable value of Mkk — 3 
TeV. 

In the RS model with SM fields in the bulk, one can generate through a simple geometrical 
mechanism the large mass hierarchies prevailing among SM fermions [164] by placing them 
differently along the extra dimension: their different wave functions overlap with the Higgs 
boson, which remains confined on the so-called TeV-brane for its mass to be protected, generate 
hierarchical patterns among the effective four-dimensional Yukawa couplings. In this case, the 
KK gauge bosons dominantly couple to heavy SM fermions as they are localized toward the 
TeV-brane (this typical feature can only be avoided in some particular situations [165]). In this 
case, the processes involving the third generation b and t quarks are those which are expected 
to be significantly affected by the presence of the new vector states. 

In the extended gauge symmetry originally proposed in Ref. [16], aU(l)B-L group was 
included; other U(l)x groups with different fermion charges can be considered [107, 166]. An 
important motivation would be that with specific charges of the new Abelian group and for 
specific fermion localizations, the three standard discrepancy between the forward-backward 
asymmetry Ap^ in Z — > 66 decays measured at LEP and the SM prediction [167] is naturally 
resolved, while keeping all the other observables in agreement with data. There are also con- 
straints on the t, b couplings from flavour changing neutral current processes, but those can be 
satisfied for Mkk values around the TeV scale [168-175]. 
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At the LHC, one can produce directly the KK excitations of the gluons and the elec- 
troweak gauge bosons with masses in the multi-TeV range [134, 137, 142, 149, 176]. In this 
note, we study the main production mechanism, the Drell-Yan channel pp — > Vkk = 
Z^^\ Z'^^\ and perform a Monte Carlo simulation of the process qq ti. We base our analysis 
on the framework which resolves the A^p^ anomaly [166], but the results that we obtain can be 
easily generalized to other scenarios. Including only the dominant QCD backgrounds, we show 
that for a set of characteristic points of the parameter space, the exchange of KK gauge bosons 
can lead to visible deviations with respect to the SM production rates. 

2. PHYSICAL FRAMEWORK 

We consider the RS model in which SM fields propagate along the extra spatial dimension, like 
gravity, but the Higgs boson remains confined on the TeV-brane. In the RS scenario, the warped 
extra dimension is compactified over a 5'^/Z2 orbifold. While the gravity scale on the Planck- 
brane is Mp = 2.44 x 10^^ GeV, the effective scale on the TeV-brane, = e-'^^^'^Mp, is 
suppressed by a warp factor which depends on the curvature radius of the anti-de Sitter space 
1/k and the compactification radius Re- The product kR^ ~ 11 leads to = 0{1) TeV, thus 
addressing the gauge hierarchy problem. The values for the fermion masses are dictated by their 
wave function localization. In order to control these localizations, the five-dimensional fermion 
fields ^i, with i = 1,2,3 being the generation index, are usually coupled to distinct masses rrii 
in the fundamental theory. If = sign(y)QA;, where y parameterizes the fifth dimension and Cj 
are dimensionless parameters, the fields decompose as ^i(a;'^, y) — Y^^=o V'i"^ i^'^)fniy)^ where 
n labels the tower of KK excitations and /o(|/) = e^^"'^*^'^'^' /N^ with A^^q being a normalization 
factor. Hence, as q increases, the wave function /o(y) tends to approach the Planck-brane at 
y = 0. 

We consider the scenario 'RSb' developed in Ref. [166] where, in order to protect the 
electroweak observables against large deviations and, at the same time, resolve the anomaly in 
Ap^, the electroweak gauge symmetry is enhanced to SU(2)l x SU(2)r x U(1)x with some 
specific fermion representations/charges under the group gauge. The usual symmetry of the SM 
is recovered after the breaking of both SU(2)r and U(l)x on the Planck-brane, with possibly 
a small breaking of the SU(2)r group in the bulk. Note the appearance of a new Z' boson (but 
without a zero-mode) which is a superposition of the state associated to the SU(2)r group 
and B associated to the U(l)x factor; the orthogonal state is the SM hypercharge B boson. 
Solving the gauge hierarchy problem forces the masses of the first KK excitations of the SM 
gauge bosons, Mkk = M^w = Mg(i) ~ M^(i) ~ Mz>{i), to be (9(TeV) and we will fix the 
common mass value to Mkk = 3 TeV in the present study. 

The light SM fermions [leptons and first/second generation quarks] are characterized by 
ciight > 0.5, c being the parameter which determines the fermion localization to cope with 
high-precision data. The large value of the top quark mass requires ct^ < 0.5 and Cgs < 0.5, 
with cq3 = ctj^ — cj,^ [as the states bi and belong to the same SU(2)l multiplet] so that 
the top and bottom quarks have to be treated separately. In the framework of Ref. [166], the 
precision data in the b sector, that is Ap^ and R^, are correctly reproduced with Mkk = 3 
TeV and e.g. gz' = 0.3\/47r for the coupling of the new Z' boson. Then, the best fit of Rb and 
^fb{V^) (which also lead to the correct range for the top and bottom quark masses) is obtained 
for cq3 ~ 0.36 and c^^ ~ 0.135. One gets the values (5(c6^) — 3.04 and Q'{cbj^) ~ 3.19, where 
Q{c) (Q'{c)) is the ratio of the four-dimensional effective coupling between the g'^^y^'^^yz^^^ 
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boson and the SM fermions, over the coupling of the gluon/photon/Z (would be Z') 
boson zero-mode. A small Q charge holds for the light fermions, (^(ciight) — —0.2, as well 
as zero Z' charge, Q'(ciight) = 0. This means that the couplings of the KK excitations of the 
gluon, photon and Z boson are an order of magnitude smaller to light fermions compared to the 
couplings to top and bottom quarks, while the KK excitations of the Z' boson do not couple to 
these light fermions at all. 

From the previous discussion, one concludes that in the chosen scenario, the sum of the 
branching ratios BR(VkK ^ ti) and BR(Vkk ^ hi), with Vkx = 7^^\ Z^'^\ Z'^^^ and g^^\ 
is close to unity which means that the KK excitations decay almost exclusively into the heavy 
t, b quarks and that little room is left for decays into light quarks and leptons. In the case of 
(yf(^) for instance, one has BR(5f*^^^ ^ tt) = 0.69 and BR((?'^^) bb) = 30. Because of the 
large couplings to fermions, the total decay widths of the KK excitations Vkk [which grow 
proportionally to the mass Mvj^j^ > 3 TeV] are very large. For instance, the decay width of g^^^ 
is of the order of a few hundred GeV and is between 10% and 20% of its mass; the KK state 
can be thus considered as a relatively narrow resonance. 

3. PRODUCTION CROSS SECTIONS AT THE LHC 

The most straightforward way to produce the KK excitations of the gauge bosons Vkk at the 
LHC is via the Drell-Yan process, pp ^ qq ^ Vkk QQ , Q = t,b with Vkk subse- 
quently decaying into top and bottom quarks. The relatively small couplings of the initial 
quarks q = it, s, c to Vkk lead to smaller production rates compared to, for instance, the 
production of Z' bosons from GUT's. Since the KK gauge bosons have different couplings to 
left- and right-handed fermions, one expects the produced t/b quarks to be polarized and to 
have a forward-backward asymmetry. Because Vkk have substantial total decay widths, the 
narrow width approximation in which the production and decay processes are factorized is not 
sufficient and one needs to consider the virtual exchange of Vkk in which the total width is 
included in a Breit-Wigner form, together with the exchange of the zero modes; the full inter- 
ference should be taken into account. For bb production, the subprocesses are also initiated by 
bottom partons and one should also consider the channel in which Vkk are exchanged in the 
t-channel. 

The signal qq — > Vkk QQ and the main SM background qq QQ and gg — > QQ 
[which gives a much more substantial contribution] have to be considered simultaneously. For 
the signal reaction, we have calculated the matrix element squared of the process pp QQ 
with polarized final state quarks and incorporated the exchange [including the t-channel contri- 
butions] of all the SM gauge bosons as well as their KK excitations and those of the Z' boson; 
we use the CTEQ5M1 set of parton distributions with the factorization and a renormalization 
scales set to the invariant mass of the QQ system, np — n^^ — mgQ .The significance S of the 
signal in the RS model can be then defined as Sc = (ct^^+sm _ ^sm) /(^sm)i/2 ^ £i/2 ^ith C 
the total LHC luminosity. 

In order to enhance the signal, which is peaked at rnqQ and to suppress the continuum 
background, one needs to select events near the KK resonance. We thus impose the cut [rngg — 
Mkk\ < Tyj^j^. To further reduce the backgrounds, we also impose the following cuts on the 
transverse momenta of the two final jets and their rapidity p^'*^ > 200 GeV and |?7q,q| < 2, 
as in the signal, the pt of the jets is peaked close to \Mv^^ and the production is central, 
while in the background, the jets are peaked in the forward and backward directions and the 
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bulk of the cross section is for low pr jets. These cuts can certainly be optimized but in this 
preliminary and simple parton-level investigation, we will simply compare the signal and the 
main corresponding physical background to determine if, grossly, one can have a detectable 
signal. More efficient and realistic cuts and detection efficiencies will not be discussed. 

The invariant mass distribution da/dmu of the process pp tt is shown in the left-hand 
side of Fig. 1 for the chosen scenario with Mkk = 3 TeV, including the cuts mentioned above. 
As can be seen, there is a substantial contribution of the KK excitations to the invariant mass 
distribution, in particular around the peak niu ~ 3 TeV. At higher m^f, the KK contribution 
becomes small, while at lower mu, it is significant even for rritt ~ 2 TeV; only for nitt ^ 1 TeV 
the KK contribution becomes negligible. Outside the KK mass peak, the RS effect is mostly due 
to the interference between the excited state and SM contributions; this interference is positive 
below and negative above the peak. The dominant contribution compared to the SM case is 
by far due to the exchange of the excitation of g'^^^ which has the largest (QCD versus EW) 
couplings to the initial state partons; the contributions of 7'-^^ and Z^^"* increase the peak only 
slightly. In turn, Z'^^'' has a negligible impact as it does not couple to the initial light quarks and 
the parton density of the heavier bottom quark in the proton is small. 

The significance of the excess of events in the RS scenarios when all KK excitations are 
included is large, S^q ~ 30 for a moderate luminosity and Sf§^ ~ 95 for a high luminosity. 
Since the excess over the SM background is mainly due to g'^^^ exchange, the significance 

iSfoo' ~ 90, is almost the same as when the full signal is considered, S^qq. In the case where 
only the first KK excitation of the photon or the Z boson is considered (assuming that the peaks 
can be disentangled, which is not obvious), the significance is much smaller SJ^^ ~ 5 and 
SiqI ~ 7. This is a mere consequence of the fact that the EW 7*^^^ Z'^^^ couplings are much 
smaller than the g'^^'> QCD couplings, leading to limited production cross sections. The smaller 
rates are, however, partly compensated by the smaller total decay widths . 




Figure 1 : The invariant mass distribution of the cross sections for the pp — > tt (left) and bb (right) processes for 
the scenario E2 discussed in the text including the RS signals and the SM backgrounds with the relevant cuts 

The discussion for the pp bb process is quite similar to the one of pp — > ti except 
for the fact that the small and not peaking t-channel bb bb contribution. The invariant 
mass distribution dcr/dm^^ for bottom quark pair production is shown in the right-hand side 
of Fig. 1 for ruf^i between 2 and 4 TeV with the relevant cuts. Here, we simply show the SM 
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background and the signal excess in the case where the contributions of all KK excitations are 
simultaneously included; again, this excess is largely dominated by the exchange of g'^^\ The 
signals are less striking than in the pp — > tt case, the main reason being that BR(g'(^) — > bb) 
is smaller than BR(5f(^) tt). However, the significances are large enough, <S^^ ~ 15 and 
iSfoo ^ 50, to allow for detection if no other background or experimental problem is included. 
Note that when combining the pp — > ti and pp — > bb processes, one would in principle be able 
to have access to the couplings of the KK states (at least g'^^^ to t, b quarks. 

Note that we have also calculated QCD higher order processes to the heavy quark pro- 
duction involving exclusively KK gauge couplings to these heavy quarks, which are the fa- 
vored couplings in the present framework. First, we have studied the one loop level reaction 
99 Vkk tt where the anomalous 99VKK vertex was regulated via the Stiickelberg mixing 
term. Secondly, we have analyzed the four-body reactions pp — > tibb, titi, bbbb and three-body 
reactions gb — > btt, bbb. Our result there is that the RS effects will be difficult to test at the 
LHC, due to the small amplitudes involved relatively to the SM background. 

Each of the KK excitations of the gauge bosons, g^^\^^^^ and Z^^\ has a different cou- 
pling to the right- and left-handed top quarks [which are themselves different from the SM 
ones]. These couplings appear in the forward-backward asymmetry as well as in the polariza- 
tion of the produced top quarks. While the enhancement in the production rate due to a single 
KK gauge boson is proportional to the sum of squared couplings, the polarization and forward- 
backward asymmetries are proportional to their difference. Thus, a combined measurement of 
the cross section together with asymmetries would determine the couplings of the vector boson 
Vkk- However, since the process is mediated by the exchange of several KK gauge bosons with 
differing right- and left-handed couplings, it will be challenging to measure precisely these 
couplings for each of the KK excitations. This is particularly true as the major contribution to 
the total rate for a{pp QQ) is coming from g'^^\ as the contribution from the electroweak 
excitations 7^^^, Z^^^ and Z'^^^ is relatively small. Nevertheless, the measurement of the polar- 
ization and forward-backward asymmetries for top quarks can be instrumental in establishing 
the presence of parity violating KK gauge bosons. 

4. CONCLUSION 

We have considered the version of the RS model with SM fields in the bulk, concentrating on 
quark geometrical localizations and gauge quantum numbers that allow to solve the the LEP 
anomaly on the forward-backward asymmetry Ap^. We have studied the main potential effects 
from KK excitations of gauge bosons at LHC, which come from their exchange in the pair 
production of third generation quarks constituting new contributions with respect to the pure 
SM ones. Based on the computation of cross sections and estimations of the SM backgrounds, 
it has been shown that simple kinematical cuts permit to detect the excesses of events in tt 
production due to the KK resonances (essentially the KK gluon) for some characteristic points 
of parameter space and A^kk ft 3 TeV (to satisfy EW constraints). For the case of bb production, 
large significances are also obtained. Furthermore, relevant top quark polarization and angular 
asymmetries were computed and turn out to provide a good test of the chiral structure for top 
quark couplings to KK gauge modes. In this simple parton-level analysis, we did not take 
into account the non-leading and non-physical backgrounds as well as as detector efficiencies, 
etc. An implementation of this process in a Monte Carlo generator, interfaced with a detector 
response simulator, have started at this Les Houches workshop [177]. 
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Part 11 

LHC studies of the left-right twin Higgs 
model 

X. Miao, S. Su, K. Black, L. March, S. Gonzalez de la Hoz, E. Ros and M. Vos 
Abstract 

The twin Higgs mechanism has recently been proposed to solve the 
little hierarchy problem. We study the LHC collider phenomenology 
of the left-right twin Higgs model. We focus on the cascade decay of 
the heavy top partner, with a signature of multiple b jets + lepton + 
missing energy. We also present the results for the decays of heavy 
gauge bosons: Wh tb, Wh — > 0^0*^ and Zh — > e+e". 

1. INTRODUCTION 

Naturalness requires the stabilization of the Higgs mass against large radiative corrections. The 
scale of new physics needs to be around the electroweak scale to avoid the fine-tuning of the 
Higgs potential. On the other hand, electroweak precision measurements push the cutoff scale 
for new physics to be likely above 5—10 TeV. This conflict in the energy scale of new physics 
is the so-called 'little hierarchy' problem. 

Recently, the twin Higgs mechanism has been proposed as a solution to the little hierarchy 
problem [178-183]. Higgses emerge as pseudo-Goldstone bosons once a global symmetry is 
spontaneously broken. Gauge and Yukawa interactions that break the global symmetry give 
masses to the Higgses, with the leading order being quadratically divergent. Once an additional 
discrete symmetry (twin symmetry) is imposed, the leading one-loop quadratically divergent 
terms respect the global symmetry. Thus they do not contribute to the Higgs masses. The Higgs 
masses do obtain one-loop logarithmically divergent contributions, resulting in masses around 
the electroweak scale when the cutoff is around 5—10 TeV. 

The twin Higgs mechanism can be implemented in left-right models with the discrete 
symmetry being identified with left-right symmetry [181]. Many new particles which have 
order of one interaction strength with the Standard Model (SM) sector are predicted and rich 
phenomenology is expected at the Large Hadron Collider (LHC). 

This paper is organized as follows. In Sec. 2, we describe the left-right twin Higgs (LRTH) 
model briefly. In Sec. 3, we present the results on the LHC studies of the cascade decay of the 
heavy top partner. In Sec. 4, we present the results on two channels for the heavy counterpart of 
the charged gauge boson: Wh tb and Wh 0^0° and the decay Zh — > e+e" of the neutral 
heavy gauge boson. Finally, in Sec.5, our conclusions are presented. 

2. THE LEFT-RIGHT TWIN HIGGS MODEL 

The LRTH model was first proposed in Ref. [181] and the details of the model as well as the 
Feynman rules, particle spectrum, and collider phenomenology have been studied in Ref. [184, 
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185]. Here we briefly introduce the model and focus our attention on the heavy top partner and 
heavy gauge bosons. 

In the LRTH model proposed in Ref. [181], the global symmetry is U(4)xU(4), with 
the diagonal subgroup of SU(2)l x SU(2)r x U(1)b_l gauged. The twin symmetry which is 
required to control the quadratic divergences of the Higgs mass is identified with the left-right 
symmetry which interchanges L and R. For the gauge couplings g2L and g2R of SU(2) ^ and 
SU(2)/i;, the left-right symmetry implies that g2L = g2R = 92- 

Two Higgs fields, H = {Hi, Hr) and H = {Hi, HrY, are introduced and each trans- 
forms as (4, 1) and (1, 4) respectively under the global symmetry. Hl^r (Hl,r ) are two com- 
ponent objects which are charged under SU(2) ^ and SU(2)r, respectively. Each Higgs obtains 
a vacuum expectation value (vev): (H) = (0, 0, 0, /), {H) = (0, 0, 0, /) with f ^ f, breaking 
one of the U(4) to U(3), respectively. The Higgs vevs also break SU(2)i^ x U(1)b_l down to 
the SMU(l)y. 

Below the cutoff scale A, the effective theory can be described by a nonlinear sigma 
model of the 14 Goldstone bosons. After spontaneous global symmetry breaking by / and /, 
three Goldstone bosons are eaten by the massive gauge bosons and Zh, and become their 
longitudinal components. The masses of the heavy gauge bosons are given approximately by 

V 2 sm 6w V cos 20^ 

After the SM electroweak symmetry breaking, three additional Goldstone bosons are eaten by 

the SM gauge bosons and Z. With certain re-parametrizations of the fields, we are left with 
four Higgses that couple to both the fermion sector and the gauge boson sector: one neutral 
pseudoscalar 0°, a pair of charged scalars 0=*", and the SM physical Higgs hsM- In addition, 
there is an SU(2)i doublet h that couples to the gauge boson sector only. It could contain a 
good dark matter candidate [186]. 

The fermion sector of the LRTH model is similar to that of the SM, with the right handed 
quarks (ur, (Ir) and leptons (/r, vr) form fundamental representations of SU(2)/j. In order to 
give the top quark a mass of the order of the electroweak scale, a pair of vector-like quarks 
and qR are introduced. The mass eigenstates, which contain one SM-like top t and a heavy top 
T, are mixtures of the gauge eigenstates. Their masses are given by 

mt ~ yv/ v^, niT ^ yf. (2) 

The mixing angle q.l and aR are controlled by the mass mixing term MqiqR. The collider 
phenomenology differs significantly for a very small value of M < 1 GeV or for not so small 
values of M. In our analysis below, we assume that M = 150 GeV. 

The new particles in the LRTH model are: heavy gauge bosons Zh, Wh, heavy top quark 
T, neutral Higgs (fP, a pair of charged Higgses 0=*=, and an SU(2)i complex Higgs doublet h. 
The free parameters in the model that are relevant for the collider studies are (/, A, M, ji^), 
where / is the vev for Higgs boson^ H, A is the cutoff scale, M is the top quark vector singlet 
mass mixing parameter, and /x^ is the mass parameter for (fp. Table. 1 shows the masses for the 
new particles in the model for several benchmark values of /. 

'The introduction of H and the requirements that it couples to gauge boson sector only and has a vev f ^ f 
are due to the electroweak precision constraints. 

^The vev / can be determined by minimizing the Coleman- Weinberg potential for the SM Higgs and requiring 
that the SM Higgs obtains an electroweak symmetry breaking vev of 246 GeV. 
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Table 1 : Mass spectram and the LHC production cross sections for Tj, Wh and Zh in the LRTH model for several 
benchmark values of /. The other parameters in the model are chosen as A = 47r/, M = 150 GeV, and ji^ =50 
GeV. 



3. CASCADE DECAY OF HEAVY TOP PARTNER 

The dominant production mode for the heavy top T at the LHC is a single heavy top quark 
produced in associated with a jet (most likely a 6 jet). For a heavy top mass of 500—1500 GeV, 
the cross section is in the range of 7 x 10^ fb — 10 fb, (predominantly via the on-shell decay 
of Wh)- The LHC cross section for Tj associated production for several benchmark points are 
given in Table. 1 . The cross section for QCD pair production is about a factor of five smaller 
due to the heavy top quark mass. 

The dominant decay channel for the heavy top is T — > Considering the subsequent 

decay of 

(f)+^tb, t^W+b^ l+ub, (3) 

and taking into account the additional energetic jet (most likely a 6-jet) that accompanies T 
from single heavy top production, the signal is typically four 6-jets + one charged lepton (e or 

fx) + missing Et. The SM backgrounds dominantly come from tt, Wjjjj, Wcjjj, Wccjj, 
and Wbbjj. Two independent studies have been performed to identify this process at the LHC. 
Their procedures, cuts and results are summarized below. 

The study by Miao and Su used MadGraph [35] and BRIDGE [84] to generate the sig- 
nal processes. Background tt events are generated using Madgraph while W+ jets events 
are generated using Alpgen [77]. Both the signal and background events are passed through 
PYTHIA [75] and PGS4 [71] for hadronization and detector simulations. For /=600 GeV, they 
have adopted the following cuts: 

• At least three jets with px > 30 GeV, with leading jet has pr > 400 GeV, and the second 
leading jet has px > 250 GeV. 

• One energetic lepton (e or /x) with px > 30 GeV and \r]\ <2.5. 

• Missing £'t > 15 GeV. 

• Reconstructed transverse top mass (from blu ) within rrit ± 20 GeV. 

• Reconstructed transverse 0^ mass within m^±± 30 GeV. 

• Reconstructed transverse heavy top mass within myi 50 GeV. 

• At least one jet is tagged as b jet. 

Similar cuts are imposed for / = 1000 GeV and / = 1500 GeV. The results for the signal 
and background cross sections before and after the cuts are summarized in Table. 2. For LHC 
integrated luminosity of 30 fb~^, the significance for heavy top discovery is more than 10 cr for 
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/ = 600 and 1000 GeV. For /=1500 GeV, more luminosity is required to reach a significant 
discovery. The significance level could be further increased with higher luminosity or when 
both detectors at the LHC are taken into account. 
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Table 2: Results for the cascade decay of heavy top T from two studies. The LHC integrated luminosity is taken 
to be L = 30 fb~^. In the analysis of Miao/Su all contributions to the T cross-sections are included, while in the 
analysis of Vos et al. only the pp Wh Tb is considered. 



In the study by Vos et al. the signal production processes and some key decay modes 

are implemented in Pythia [75]. The response of the ATLAS detector is simulated using the 
ATLAS fast simulation package ATLFAST [73]. The decay is reconstructed step by step, start- 
ing from the W decay into lepton neutrino. At each step kinematical constraints - on the mass 
and transverse momentum of the reconstructed particles - are applied. The values of the cuts 
employed in the selection vary with model parameter / that governs the masses of the involved 
particles. In the following, the selection criteria for / = 555 GeV (i.e. a reconstruction aimed 
at a 1.25 TeV Wn-hoson) are given. 

• A lepton (e"^, fi^) with transverse momentum greater than 25 GeV. The presence of the 
lepton ensures that the events can be triggered efficiently. 

• A minimum missing transverse energy of 25 GeV. A W candidate is reconstructed from 
the missing transverse energy and the lepton momentum using the coUinear approxima- 
tion (i.e. assuming = ). 

• The W candidate is combined with all jets with 25 < pr(j ) < 200 GeV. The combination 
that gives the best match with the top mass is selected. If none of the combinations yields 
a mass < 250 GeV, the event is discarded. 

• A second jet with 25 < pr < 100 GeV is added to reconstruct the charged Higgs boson 
0^. Again, events with a reconstructed 0^ mass greater than 250 GeV are discarded. 

• A third jet with pr(j) > 100 GeV is required to reconstruct the heavy top quark T. 
The T-candidate is required to satisfy the following constraints: mr < 700 GeV and 
Pt{T) > 150 GeV. This latter cut, that takes advantage of the Jacobean peak in the signal, 
is particularly useful to reduce the dominant ti background. 

• Finally, a fourth jet with prij) > 150 GeV is used to form the Wh candidate. As the 
dominant T production process is through Wh Tb, the explicit reconstruction of the 
Wh is instrumental in reducing the background. 

The width of the reconstructed mass peaks is dominated by the experimental resolution for jet 
energy and missing transverse energy. For the lightest Wh boson of 1.25 TeV, the Gaussian 
width of the peak is approximately 100 GeV, large compared to the natural width of 30 GeV. 
The total efficiency for the kinematical reconstruction is 9 %. 

The tt and W+ jets backgrounds are generated using Pythia. The former is found to be 
the dominant contribution to the background after the kinematical reconstruction is performed. 
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At this stage, a significance (when estimated as S/ \fB ) well above 5 is found for all mass 
points considered. The shape of signal and background mass distributions are, however, quite 
similar. Experimental methods to normalize the background, like the sideband estimate, cannot 
readily be employed. Taking into account the large uncertainty on the background production 
cross section and selection efficiency, the significance of the signal is greatly reduced. 

To further reduce the background, the lifetime signature of the multiple 6-jets is used. The 
signal topology contains very highp^ &-jets. To correctly describe the ATLAS performance, jets 
with transverse momenta up to 1 TeV have been studied using a detailed GEANT4 simulation 
of the detector response [187]. Several detector effects are found to lead to a quite significant 
degradation of the performance for the highest pt bins. In this study, a parameterization of the 
full simulation results are used. 

To discriminate the signal topology, with four 6-jets, against the dominant ti background 
a four jet likelihood is constructed by summing the tag likelihood of the four leading jets. The 
four-jet likelihood allows for a significant reduction of the background. For the studied mass 
points, the signficance (S/y/B) improves slightly or remains unaltered, but the S/B is greatly 
improved, thus rendering the analysis much more robust against uncertainties in the number of 
background events. The results - listed in Table 2 - indicate that the discovery potential reaches 
a Wh mass of 3 TeV, even for a relatively small integrated luminosity (30 fb~^). For larger 
Wh masses, the discovery potential is rapidly degraded by the small absolute number of signal 
events. 

A closer look into the dominant tt background for different versions of the generator 

yields significantly different results for the high-p^ tail of the top quark spectrum, and therefore 
of the number of background events that pass the kinematical reconstruction and selection cuts. 
Therefore, the ATLAS study is being repeated using the MC@NLO [188, 189] generator. These 
results will be published at a later time [190]. 

The heavy top could also decays into bW, thgu, tZ, and t(jP. The branching ratios for 
those channels quick drop for larger mr and smaller M. Therefore, we will not discuss those 
channels further here. 

4. HEAVY GAUGE BOSONS 

The dominant production channels for heavy gauge bosons at the LHC are the Drell-Yan pro- 
cesses: pp — > WhX and pp — > ZhX. The production cross sections for several benchmark 
points are given in Table. 1. 

The dominant decay modes for Wh are into two jets, with a branching ratio of about 60%. 
Such modes suffer from the overwhelming QCD di-jets background for large pt jets. Wh could 

also decay into a heavy top plus a &-jet, with a branching ratio of about 20%— 30%. This is the 
main channel for Tj associated production as discussed earlier in Sec. 3. 

The decay Wh ^ (p^cj)'^ has a branching fraction of 3%. The subsequent decay of the 
charged and neutral Higgs bosons 0=*" — > th and 0° — > hh have large branching ratios. Thus, 
the same final state with four 6-jets and a lepton and neutrino as for the decay Wh — > Th 
is obtained. A kinematical reconstruction of the decay chain along the lines of the previous 
analysis described in Sec. 3 is quite successful for small values of /. For larger Wh mass the 
boost of the relatively light Higgs bosons increasingly leads to difficulties in the reconstruction 
of the jets. The reduction in reconstruction efficiency, in combination with the sharply dropping 
production cross section, limits the discovery range of this signature to Wh masses below 1.5 
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- 2 TeV. 

The branching ratio for Wh tb is of the order of 4%. Due to the much reduced 
branching ratio this channel is a priori less promising than in the Littlest Higgs model, studied 
by ATLAS [191]. Recent work [190] investigates the reduction of the dominant tt background 
by exploiting the presence of additional jets or leptons in the background sample. Preliminary 
results indicate that isolation of this signal for Wh masses up to 1.5 — 2 TeV may well be 
possible. 

Although the dominant decay modes of Zh are into dijets, the discovery modes for Zh 
would be Zff — > (with a branching ratio of 2.5% for e+e^, and t+t~ individually). 

The natural width of the heavy Zh ranges from 25 to 75 GeV for a Zh boson from 1.2 to 3.6 
TeV. The di-lepton modes e+e" and therefore provide clean signatures, which can be 

separated from the SM background by studying the invariant dilepton mass distribution. 

A study of the di-lepton signature has been performed by ATLAS. The excellent momen- 
tum resolution for high pt electrons yields an error in the invariant mass of the parent boson 
that is inferior to the natural width throughout the studied mass range. For the di-muon final 
state the invariant mass resolution is limited by the resolution of the combined measurement 
of inner tracker and muon spectrometer. Therefore, this first exploration only considers the 
di-electron signature. The discovery potential of this channel is evaluated using a classical anal- 
ysis counting the (small) number of signal and background events in a narrow mass window. 
This signature may give rise to very early discovery, with only a few inverse fb of data, pro- 
vided the Zh mass is less than 2.5 TeV. With an integrated luminosity of 75 fb~^ the discovery 
reach is extended up to 3.5 TeV. Potentially, the LHC experiments are sensitive to much larger 
masses through the interference of the heavy neutral gauge boson with the Standard Model Z 
and photon. 

Zh could also decay into ti final states, with a branching ratio of about 2 — 5%. Searches 
of tt resonance have been studied in Refs. [37, 192]. The reach for Zh — > tt at the LHC is very 
Umited, due to the small decay branching ratio into the tt final states. 

5. CONCLUSIONS 

The twin Higgs mechanism provides an alternative method to solve the little hierarchy problem. 
The LRTH model has rich collider phenomenology that could be studied at the upcoming LHC. 
In this paper, we presented LHC studies on the searches of the heavy top partner and the heavy 
gauge bosons in the LRTH models. For the heavy top quark partner, a significance level of over 
3 a could be reached for almost the entire interesting parameter regions of the LRTH model 
with an integrated luminosity of 30 fb~^. An independent study in the ATLAS fast simulation 
framework has used a parameterization of the full simulation results for high pt 6-tagging. This 
study finds a significant signal up to / = 1100 GeV with an integrated luminosity of 30 fb~^. 
Several Wh decays have been studied. The decays Wh 0^0° and Wh — > tb give rise to a 
final state with a lepton and a neutrino and four, respectively two, &-jets. The search for both 
channels is expected to yield a significant signal only for small values (of the order of 1 TeV) 
of the Wh mass. Zh e+e" is likely to be the discovery channel for the LRTH model. For 
relatively light Zh (up to 2.5 TeV) this signature could be observed with only a few inverse fb 
of data. The discovery potential for the LRTH model at the LHC is very promising and further 
studies towards the identification of the twin Higgs mechanism and the distinction between 
various electroweak models are currently under investigation. 
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Part 12 

WiWi Scattering 

A. Delgado, C. Grojean, E. Maina and R. Rosenfeld 
Abstract 

In this report we intend first to review the main models where strong dy- 
namics are responsible for EWSB. An overview of tests of new models 
through the production of new resonant states at the LHC is presented. 
We illustrate how different models can be related by looking at two gen- 
eral models with resonances. 

1. INTRODUCTION 

One of the main goals of the LHC is to find the mechanism responsible for electroweak sym- 
metry breaking (EWSB) at the TeV scale. In the SM of electroweak interactions, this is accom- 
plished by postulating the existence of a complex scalar Higgs field with a potential crafted in 
such a way as to result in the breaking of SU{2)l x U{1)y into the residual electromagnetic 
U{1) symmetry. The couplings of the Higgs field with gauge bosons and fermions generate the 
masses we observe. This rosy picture has its thorns: a scalar sector is unstable under radiative 
corrections and the Higgs sector has to be understood as an effective theory valid up to some 
energy scale A. 

One of the ideas that has been intensely studied in the past is that there are new strong 
interactions responsible for EWSB and the Higgs sector effective description will break down 
at a scale A ~ TeV. Some of these ideas have resurfaced recently in the form of various models. 
The longitudinal components of the gauge bosons have their origin in the EWSB mechanism 
and hence provide a window to study these models. 

The charge of our working group is to review the large amount of work that has been done 
in the past and to identify improvements that can be made in the light of these novel models 
of strongly interacting EWSB sector and of the new tools available to study possible signatures 
at the LHC. In particular, we want to concentrate on the more model-independent features of 
these models. 

There are three main themes we would like to address in this report: 

• WW scattering and unitarity; 

• Drell-Yan versus WW fusion as discovery processes at the LHC; 

• constraints from electroweak precision measurements. 

2. MODELS 

Heavy, relatively wide resonances are the hallmark of strong interactions. They may or may 
not be at the LHC reach for direct detection, depending on their masses. Hence it is convenient 
to classify different models of strong interactions in terms of having light or heavy resonances, 
where light means within LHC reach. It is also convenient to separate out the scalar from 
the vector resonances, since they are expected to have a very different phenomenology. An 
incomplete list of models is classified in Table 1. This classification is of course arbitrary. In 
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fact, many of these models can actually move among different classes for a different set of 
parameters. For an alternative classification of models, see e.g. Cheng [193]. 





No (light) vector resonances 


Light vector resonances 


No (light) Higgs 


Chiral lagrangians [194, 195] 


LSTC [53], Higgsless [15, 17] 
(D)BESS [33] 


Light Higgs 


SM, SILH [196] 


Warped/Composite [197] 

Holographic [198], Little [124, 125], 
Gaugephobic [199], Twin [178, 181] 
LDBESS [200], Gauge-Higgs [201-204] 



Table 1: Attempt to classify different models of strong dynamics 



Many of these non-SUSY models share the same low-energy phenomenology. Models 
based on extra dimensions, such as Gauge-Higgs unification in either flat or warped extra di- 
mension, may resemble conventional 4D models with new particles such as Little Higgs models 
since only the lightest KK modes will be accessible at the LHC. In fact, they can both be de- 
scribed by a so-called three-site moose model [205]. Also, warped Higgs models and Higgsless 
models can be smoothly interpolated in the so-called gaugephobic Higgs model [199]. In the 
latter Sections we will compare 2 general models with new scalar and vector particles, namely 
LDBESS and CompositeAVarped Higgs models. 

We will focus primarily on the optimistic view that new resonances are within the LHC 
reach and construct a general model that captures some general features of the phenomenol- 
ogy of these resonances. In particular, we will not consider here the more model dependent 
possibility of composite fermions, since we will be interested mostly in the interaction of these 
resonances to light fermions, which is important for Drell-Yan production, for example. If the 
resonances are heavy and outside the LHC reach one will have to test their indirect effects in 
anomalous gauge boson couplings coming from higher dimensional operators. 

We will denote generically by V these new vector resonances. They will in general mix 
with the SM gauge bosons. They may or may not have direct couplings to light fermions. In 
the latter case, the couplings will arise solely from the mixings with gauge bosons. The models 
will then be mainly characterized by: 

• the mass eigenstates and My of the gauge bosons and resonances; 

• the couplings gvww and ghww between gauge bosons and resonances and the Higgs boson 
arising from the mixings. 

In the SM the WlWl scattering ampUtudes (here W generically denotes massive gauge 
boson, the as well as the Z) are unitarized by the Higgs boson contribution. If the Higgs 
coupling ghww is modified, unitarity is not exactly restored by the Higgs boson alone and 
the contribution from other resonances may be relevant. In models with no Higgs at all, the 
contribution arising from a tower of resonances can help to either unitarize or delay unitarity 
violation to higher energies. 

2.1 Electroweak constraints 

We will study the conditions imposed on these models by the electroweak observables. This 
goes as far as the early nineties when the oblique parameters S, T and U [206] and the e's [207] 
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where introduced, it was shown that minimal technicolor had difficulties accommodating the 
results from the different accelerators running at that time, specially SLC and LEP. 

The effects of new physics can be classified as "oblique" or universal with respect to the 
fermionic currents and those which affect differently to different flavors. The second ones tend 
to be more restrictive because they imply FCNC and they are also more model dependent so we 
will only describe here in some detail the first ones because, in general any model predicting 
FCNC will have great difficulties accommodating actual data . 

These universal effects can only affect the self energies of the electroweak bosons since 
anything that affect equally all fermionic currents can be move into the kinetic term for the 
gauge bosons via field redefinition. Those self energies can be expanded in powers of the 
momentum, following the notation of [208]: 

liABip') = n^B(o) +p'iIab{o) + ^(/)'n:;^(o) + . . . (i) 

where AB = W~^, W~ , W^W^,, BB, BW^. There are twelve coefficients, three of them are 
reabsorbed into g, g' and v and other two are fixed to ensure that the photon is massless and 
couples to the correct current. We are left with seven parameters, three of them correspond to 
the usual e's or 5", T and U and four new ones only relevant for higher order corrections. 

Let as focus to the most important ones that are S or €3 and T or ei, their definition are as 
follows: 

5,63 ~ nV3B(o) 

T,ei ~ ^WzwX^) -^w+w-{^) (2) 

The latter has a direct connection to the p-parameter which measures the breaking of the cus- 
todial symmetry. In the SM that breaking is very weak leading to very small deviations from 
1 in the p-parameter. The way to avoid large contributions to the T parameter is to ensure that 
the strong interaction describing the EW breaking sector has some kind of custodial symmetry 
embedded on it [16]. This, in turn, implies that the resonances fall into representation of the 
custodial group, and also that there are vectorial right-handed resonances. 

The second big issue is that technicolor or Higgsless theories tend to predict very large 
contributions to the S parameter. In that sense having a light scalar resonance, i.e. Higgs, is 
preferred from the electroweak constraints. We can rewrite the S parameter in the following 
suggesting way: 

,5 = 47r 



(3) 



where Fi and Mi are the decay constants and masses for vectorial (V) and axial (A) resonances, 
in minimal technicolor inspired by QCD the value is 0.3, too large compare to the experimental 
data, so a way to suppress big contribution to S is to have vector and axial resonances almost 
degenerate which in turn will mean that experimentally these resonances should be discovered 
with similar masses. 

For theories where such resonances do not follow that rule, like in Ref. [27], another 
way to cancel those large contribution to the S parameter is to suppose that light fermions do 
have suppressed couplings to heavy resonances, this will make Vector Boson Fusion (VBF) the 
favorite channel to produced these resonances. 

To summarize, the electroweak data impose two main conditions into theories with strongly 
coupled sectors: 
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1. Some implementation of custodial symmetry is needed which means extra right-handed 
states in the spectrum. 

2. A particular choice of the spectrum of the resonances or cancelation with fermionic con- 
tributions is granted to reduce dangerous contributions to S. This will also have particular 
implications in the spectrum. 

Before finishing this discussion let us mention that the third family may potentially deviate 
from this picture since getting a big enough top mass can be in conflict with the requirement 
that the bottom quark couples to the Z the same way the rest of the quarks do but since we 
are just dealing with WlWl scattering those potential problems will not affect this scattering 
process. 

2.2 Strong interactions from mixings 

In this subsection we show the origin of a possible strong interaction between the vector reso- 
nance composite fields and the SM gauge fields. Let us consider the decay of a vector resonance 
to two gauge bosons. Due to the gauge structure of the vertices before mixing one would have: 

^XMX. ^ W^'{qi)W;:^iq2)) = gvww [{p - qi),gau+ (4) 

+{qi - q2)ag^. + {q2 - pUa^^] e^''^{q2) 

Using the rest frame of the decaying vector resonance and the usual representation for the 
polarization vector one finds for the decay amplitude into longitudinal polarizations of the gauge 
bosons: 

A^'^^ = gvww^iqi - q2) ■ e\ (5) 

Comparing this result with the one following from a simple QCD inspired effective la- 
grangian describing the couplings of p's to pions: 



>C ~ S'p^ttPm 9" (6) 
one immediately obtains the correspondence 

QpTTTv — gvww (7) 
and strong interactions can arise from mixing when 

My 

9VWWWV72- Z 1 (8) 



2M^ 



2.3 Drell-Yan versus Vector Boson Fusion 

In general there are two competing processes for the production of resonances at hadron collid- 
ers: Drell-Yan (DY) and vector boson fusion (VBF) processes. Let us for the moment focus on 
the production of a generic neutral resonance R at the LHC which couples to first generation 
fermions and to longitudinally polarized Wl's. The DY cross section can be written as: 

a^y^pp ^ WlWlX) -Y. f (dL/dr)^/^^^^ a{q,qj ^R^ W^W^) , (9) 
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where the sum is over quarks flavors that can produce the pair of Wl's, t = My^-^/s, Tmin — 
AM^/s, a is the partonic cross section at a center-of-mass energy of ^/TS and the partonic 
luminosity function is given in terms of the parton distribution function by: 

(dL/dr)^/,^^^ = / - [x.(«)(x)^zW(r/x) + d^'^\x)S'^\r/x) + (a) - (b)] (10) 
Jt ^ 

In the narrow-width approximation of the Breit-Wigner form for the partonic cross section one 
has 

which results in 

^DY^^ ^ W^WlX) = ^ m,)BR{R WlWl) {rdL/dr)^,^^^^ (12) 



Analogously, the VBF cross section can be estimated as 

^^T{R^WlWl)BR{R- . 



(13) 

where {dL/dT)^^^^^^-^^^ is the luminosity of a pair of Wl's inside the proton and is approxi- 
mately given by: 

/I Wr' dr 
— - [u(^\x)S'\r' M + {a) ^ {b)] {dL/dO^,/^.^-, 

(14) 

where ^ = t/t' and 

idL/dO^,/^.^- = (^)'^[(l + C)ln(l/e) + 2(e-l)] (15) 

The relative DY and VBF contributions can be easily estimated in the context of these 
approximations: 

a^BP{pp ^ WlWlX) T{R ^ WlWl) {dL/ dr)^^,^^^^ 

- (16) 



a^y^pp^WLW^X) T{R^q,q,) {dL/dr)^/^^ 



Estimating the ratio of luminosities one finds that in order for the VBF process to be 
competitive, the coupling of the resonance to light quarks must be suppressed so that the ratio 
of the partial widths can compensate for the (9(10^®) smaller luminosity. This is usually the 
case for scalar resonances but it can also be the case for some models with vector resonances. 



3. A BRIEF HISTORY OF VECTOR BOSON FUSION 

Vector boson fusion has come a long way since first proposed by Cahn and Dawson [209]. 
The development of the effective W approximation (EWA) [210], used in conjunction with the 

equivalence theorem (FT) [1,211,212], provided an important tool for estimations of the signals 
for these processes. A first realistic study of signatures of strong interactions in WW scattering 
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at the LHC was perfomied for different models in the context of EWA/ET in [38], see Fig. 1. 
For a model with vector resonances, it was recognized that the Drell-Yan production is more 
important than VBF for resonances lighter than C(TeV). 

Similar conclusions were reached in models of chiral electroweak lagrangians with reso- 
nances generated via unitarization procedures, again in the EWA/ET context [213,214]. How- 
ever, a detailed simulation for the ATLAS detector of this type of models with a 1.2 TeV vector 
resonance using VBF with EWA/ET implemented in Pythia has been reported [215,216], where 
it is claimed that 100 (300) fb~^ is needed to detect this resonance in qq ^ q'q'WZ using the 
WZ liyjj (WZ Ml) mode. 

Realistic simulations using the CMS detector have been recently performed for the case 
of a vector pre in the context of LSTC implemented in Pythia, using the dominant Drell-Yan 
process [70]. A 5a discovery is reported for a integrated luminosity of 30 fb~^ for Mp^^ up to 
700 GeV. 

The signatures for a strongly interacting sector becomes even more challenging when 
there are no resonances at the LHC reach. They manifest at low energies in anomalous gauge 
boson interactions generated by integrating out the heavy resonances. An exact tree-level 
study of anomalous quartic gauge boson couplings in VBF was performed, where bounds on 
some coefficients of the chiral lagrangian were determined [217] and refined recently in the 
Ref. [218]. The VBF production of a Higgs boson with non-standard couplings to electroweak 
gauge bosons was analysed in [40]. Significant differences between a complete calculation of 
VBF and the EWA/ET approximation were found in the large Mww invariant mass distribution 
region, see Fig. 2, where a full calculation of 2 — > 6 amplitudes was also performed [219,220]. 




Figure 1 : Signal and background for a vector resonance in WW (left panel) and Drell-Yan (right panel) (from [38]). 



Turning now to the production of resonances in recent developments, the production at 
the LHC of a neutral KK excitation of the electroweak gauge bosons, generically called Z' , 
in the context of a warped extra dimension model was studied by Agashe et al. [137]. They 
obtained that the Drell-Yan process is dominant over the VBF. The best discovery mode is 
Z' — > HZ — > bbl'^l~, see Fig. 3. They conclude that a 2 (3) TeV Z' can be discovered with 
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approximately 100 fb ^ (1 ab ^) integrated luminosity. Production of KK excitations of gauge 
bosons and its contribution to top quark pair production at the LHC was discussed in [176]. 

In the case of Higgsless models, an opposite situation occurs for the production of W [36] . 
Due to the vanishing of fermionic couplings to the new vector resonance states, which guaran- 
tees a null correction to the electroweak precision parameters, the VBF is the dominant produc- 
tion process. A detailed study, including a complete leading order computation of the signal and 
relevant background, shows that with a 100 fb^^ integrated luminosity the LHC will completely 
cover the parameter space of this model, with a5 a discovery of W up to a mass of 1.2 TeV, 
see Fig. 4. 

It is also interesting to point out that it could be possible to use gluon-gluon fusion Higgs 
production rate (and the Higgs decay rate to 2 photons) to discriminate among different models, 
such as UED and Gauge-Higgs unification models [222-224]. 

The next frontier in VBF is certainly the next-to-leading order (NLO) computations. A 
parton level Monte Carlo implementation of NLO QCD corrections to vector boson pair pro- 
duction via VBF was recently reported [225]. The implementation of new models in the same 
framework may be available in the near future. 

See also [92,226] for recent studies of strong WW scattering. 

4. RESONANCES IN A STRONGLY COUPLED EWSB SECTOR 
4.1 Inspiration from low energy QCD 

Low energy QCD (by low energy we mean ^/s -C m^) with two flavors is a theory of strongly 
coupled pions described by a chiral lagrangian (non-linear a model): 

C = ^Tr [d^Ud^U^] (17) 
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Figure 3: Cross section for Z' production at the LHC for different Z' masses, with the Drell-Yan and VBF pro- 
cesses shown separately (left panel). Z' mass reconstruction from pp ^ Z' ^ HZ — > bbl~^l~ (right panel) 
(from [137]). 



with U = e*^'^"^"/-^'^ transforming as a (2, 2) under a global chiral symmetry G = SU{2)l x 
SU{2)r: 

U^gnUgl; gReSU{2)R, Ql e SU{2)l. (18) 

This lagrangian is the lowest order term of an infinite expansion in an increasing number of 
derivatives of the field U. The pions are the Nambu-Goldstone bosons resulting from the spon- 
taneous symmetry breaking SU{2)l x SU{2)r SU{2)v generated by a non-zero vaccum 
expectation value (U) = 1. 

It follows from the lagrangian (17) that tttt scattering amplitudes grow as s/f^, as pre- 
dicted by current algebra arguments. They violate s— wave perturbative unitarity when ^/s ^ 
4:7: /t^. Higher derivative terms in the chiral lagrangian can delay the energy scale for unitarity 
violation. 

There are several ad-hoc methods to unitarize tttt scattering but we will not dwell on 
them. In QCD, violation of perturbative unitarity points to the existence of vrvr resonances, like 
the p and ai. Hence, one must extend the formalism of chiral perturbation theory to include 
these resonances in order to describe physics at energy scales ^/s nip. Possibly one of 
the most successful attempts to introduce resonances consists in describing them as dynamical 
gauge bosons of a local symmetry, as first proposed by Sakurai [227]. This idea led to the 
so-called hidden symmetry approach [228], where the masses of the resonances are generated 
via a Higgs mechanism. The hidden symmetry model can nicely accomodate features such as 
vector meson dominance and gauge invariance. It can also be used to describe resonances in a 
strongly coupled EWSB sector. 
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Figure 4: Cross section for W' production at the LHC for different W' masses, with the associated W'Z and 
VBF processes shown separately (left panel). Significance for W' detection for the different modes (right panel) 
(from [221]). 



4.2 Strongly coupled EWSB sector 

The EWSB sector in the SM is a gauged linear cr-model which can be nicely written as [194, 
195]: 

CsB = \Tr [{D^U)^D^U] - h (^-Tr [U^U] - v^^ (19) 

where 

D^U = d^U + \igW;nU - (20) 

In the limit (?' — >^ this lagrangian has an additional global symmetry G = SU (2) l x SU{2)r, 
with U transforming as a (2, 2): 

U e^''^^'/2f/e-^'«^'/2. (21) 

Electroweak symmetry is spontaneously broken by a non-zero vacuum expectation value (f/) = 
V, which generates masses for the electroweak gauge bosons. A so-called custodial global 
SU{2)v symmetry (e^ = en) survives, which is responsible for keeping the parameter p = 
M'^/Ml cos 6'vi/ = 1 at tree level. 

One way to introduce resonances from a strongly coupled sector in an effective lagrangian 
is to extend the SM global symmetry to a hnearly realized [SU{2)l x SU (2)^]''' x [SU{2)l x 
SU{2)^^°'^P symmetry. Here we are using the superscripts {el) and (comp) to denote the 
elementary and composite sectors respectively. The subgroup [SU (2)^ x U il)xY^ x [SU {2)i x 
SU{2)rY°"^'p will be gauged. The masses of all gauge bosons are generated via an extended 
Higgs mechanism. The SM Higgs can be part of the elementary sector (LDBESS model by 
Casalbuoni et al. [33, 229,230]) or part of the composite sector (2-site composite Higgs model 
by Contino et al. [197]). 
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Figure 5: Moose diagrams of the LDBESS model (left panel) and a 2-site composite Higgs model (right panel). 



4.3 LDBESS model 

The scalar sector of the LDBESS model is given, as depicted on the moose diagram of Fig.5, 
by 3 sets of fields transforming under [SU{2)l x SU{2)b]^^ x [SU{2)l x SU{2)rY"^p as 

Lg (2,0,2,0) =^ L^QLLhl (22) 
i? €(0,2,0,2) =^ R^gRRh^j, (23) 
f/G (2,2,0,0) =^ U^gLUgl (24) 

where g^ R G SU (2)f ^ and Hl^r g SU (2)^°7. Notice that the Higgs field U transforms under 
the elementary group. 

The covariant derivatives acting on the scalar fields are (suppressing Lorentz indices): 

DL = dL + tg^^W'L-tglL'^Vl (25) 
DR = dR + igJ^YR-ig*^R^-^Vk (26) 
DU = dU + tg2'^W'U-tg,U^Y, (27) 

where W, Y, Vl and Vr are the gauge bosons of the local \SU{2)l x U{1)yY x [SU{2)l x 
SU {2)r\^°'^p symmetry. The couplings of the elementary sector are gi and g2, whereas g*j^ and 
g\ are the couplings of the composite sector. 

The scalar potential is chosen in such a way that (L) = {R) = /I2 and (U) = VI2, with 
f ^ V, resulting in the symmetry breaking pattern 

[2Lxlxf X [2Lx2Rr"^P^''ti''^ [SU{2)lxU{1)y]''' ^ U{lU (28) 
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As usual, the scalar kinetic terms generate masses and mass-mixing terms for the gauge 
bosons: 



4^^ 



(29) 



with the result: 



Even in absence of EWSB (v = 0), the L link field induces a mixing between the elemen- 
tary and the composite fields. The mass eigenstates are a mixture of elementary and composite 
fields. Diagonalizing the mass matrices after EWSB, we obtain the following spectrum 

• Electrically charged sector 

- a composite massive W^^ — with a mass m^^± = lg*jif 



'hi 



- a mostly composite massive which is a combinaison of and V^- 



W^^ = sin ^(1 + r cos' e)W^ - cos ^(1 - r sin' e)V^, (30) 
with a mass given by 

Mk. - \9ff (^ + -tan'^ + 0(.Vr)) . (31) 

- a mostly elementary light which is the linear combination orthogonal to 

^ cose {l-rsm'^e)W^ + sine {l + rcos'^e)V^ , (32) 
with a mass proportional the SM Higgs vev 

M^, = -/v^ (1 - r sin' e + 0{v'/f')) . (33) 

• Electrically neutral sector 

- a massless photon 7 = F. 

- a mostly elementary light Z which is a linear combination of and Vl 

Zi = cos e{l - r sin' e)W^ + sin e{l + r cos' ^)yt^ (34) 

with a mass proportional to the SM Higgs vev and equal to Mwi 

- a mostly composite heavy Zfi2 corresponding to the linear combination orthogonal 
to Zi. Its mass is equal to Mwh2 

- a composite heavy Z^i = with a mass m^± = ^g*^/. 

The mixing angle e corresponds to the ratio of the gauge couplings of the elementary and 
composite sectors 

V 9^+91^ 

It measures the amount of compositeness in the massless gauge bosons before EWSB. g is the 
gauge coupling of the low-energy SM SU{2) gauge coupling 
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after the breaking SU{2Y^ x 5^7(2)'=°'"^ to SU{2)1^'. Note that c/ ~ ^2 if the composite sector 
is strongly coupled {g\ 512)- Finally, the parameter r measures the mass gap between the 
Ught and the heavy states 

We have neglected here the coupling gi of the elementary U{1) gauge group. Turning it on will 
induce a mixing of V^, V|, and Y . 

Trilinear couplings between the light and the heavy states arise from the gauge kinetic 
terms, 

J2 ^TrK.]. (38) 

W,Y,Vl,Vr 

For instance, we obtain the coupling 

9zhiW+W[- ^ gr sine cos 9. (39) 

Notice that the coupling is suppressed by a factor r and that, in the absence of symmetry break- 
ing, this coupling vanishes. 

This model can also be extended to include more resonances by simply introducing more 
copies of [SU{2)l x 5[/(2)i?]^»™P or [SU{3)c x SU{2)l x SU{2)bY"'"'p if color octet vector 
bosons (gluon resonances) are present in the spectrum. However, the basic idea is the one shown 
above. Notice that in the case of gluon resonances, their coupling to gluons would vanish since 
the color symmetry is unbroken, as pointed out in [231]. 



4.4 A 2-site composite Higgs model 

The basic difference between the LDBESS model and a 2-site composite Higgs is that the Higgs 
is now part of the strong sector and interacts directly with the composite gauge bosons. So 

U e (0, 0, 2, 2) =^ hiUh^'ji (40) 

and 

DU = dU + tgl^VlU - tg*nU^V^, (41) 

in the notation adopted in the previous section. One still keeps (L) — (R) — f and (U) — v. 
Consequently we now have the mass-mixing lagrangian given by: 

J^mass = I [v' {gl'Vl, + g*j,'Vl, - 2glg},VlVk) + f [glW^ + g^Vl, - 2g^glW,Vl) 

As before, the mass eigenstates are linear combinations of elementary and composite 
states. Diagonalizing the mass mixing terms, we easily obtain the spectrum after EWSB. For 
example, in the electrically charged sector, we get two massive VF^^_2 with masses given by 

= -/Lf' + ^ '^'^"^ + 0{vVn) . (43) 
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and one light Wi whose mass is proportional to the Higgs vev and that can be identified with 
the SM 

Ml.^ = \gW + 0{v'/f). (44) 
As for the LDBESS model, g is the gauge coupling of the low-energy SM SU (2) gauge coupling 

/ = 4f^2^ (45) 

gi + gl 

after the breaking SU{2Y^ x SU{2Y""''p to SU{2)1^'^ at the scale / and the mixing angle is 
the ratio of the gauge couplings of the elementary and composite sectors 

sing^ r#--f f (46) 

^gl + g*L 

Finally, r is the mass gap between the light and the heavy states 

The spectrum is the same as in the LDBESS model up to differences of 0{v'^ / p) in the heavy 
sector. Still, these differences have big effects in the couplings of the Higgs boson to the heavy 
resonances. These couplings can be simply obtained by substituting vby v-\-h'va the lagrangian: 
therefore, in the LDBESS model the Higgs coupling to VF^ are reduced two powers of the 
compositeness angle 6, while in the composite Higgs model, the coupling is enhanced by the 
same factor. 
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Part 13 

Vector-like quarks: a toolkit for 
experimenters 

J. Santiago 

Abstract 

We review the motivation and main features of vector-like quarks with 
special emphasis on the techniques used in the calculation of the fea- 
tures relevant for their collider implications. 

In four space-time dimensions, Dirac fermions can be decomposed into left- and right-handed 
chiralities, ip = ipL + ipR, with ipL,R = (1 T 7^)/2^/'. Chiral fermions, for which the two 
components have different charges under the electroweak SU (2) ^ x U{1)y gauge group, have 
masses proportional to the electroweak symmetry breaking (EWSB) scale {v ~ 174 GeV) 
and are therefore expected to be relatively light, mchirai = < v. They can only be made 
heavy at the expense of introducing a large dimensionless coupling A ^ 1, thereby inducing 
large one loop corrections to electroweak observables (they do not decouple). Thus, new chiral 
fermions are quite constrained experimentally from electroweak precision tests (EWPT). If the 
fermions are vector-like (the two chiraUties have the same quantum numbers) we can write 
down a gauge invariant (Dirac) mass that is unrelated to the EWSB scale. They can therefore be 
naturally heavier than v without introducing large dimensionless couplings. As the Dirac mass 
becomes large, their low energy effects become negligible {decoupling) and EWPT impose no 
constraints. Furthermore, new vector-like fermions are a very common prediction of theories 
beyond the Standard Model (SM), from Kaluza-Klein modes in models with extra dimensions 
to partners of the SM fermions in Uttle Higgs models, or members of extended multiplets in 
GUT theories. 

For their relevance to LHC physics, we will restrict ourselves to vector-like quarks, trans- 
forming as triplets under the SU (3) gauge group. We will also consider only the case of new 
quarks that can sizably mix (through mass terms) with SM quarks, generating flavour changing 
electroweak processes that can lead to interesting signatures at the LHC. Finally, we will assume 
that a (relatively) light SM-like Higgs boson (doublet under SU (2) exists in the spectrum and 
is the main source of EWSB. 

1. HOW TO COMPUTE THE RELEVANT FEATURES: MASSES AND COUPLINGS 

New vector-like quarks with sizable mixings with the SM quarks cannot have arbitrary quan- 
tum numbers. They can only couple in a gauge invariant way to a SM quark and the Higgs, 
which gives only the seven distinct possibilities shown in Table 1. Given these possibilities, the 
question is, how do we get the relevant features of the model in order to study it at colliders. Or 
said otherwise, what are the masses and couplings of these new quarks? The aim of this short 
review is to present the general algorithm to compute these relevant features and to exemplify it 
in a particular case that will be further explored below. The general Lagrangian can be written, 
in the current eigenstate basis {i.e. in terms of fermions with well defined SU{2)l x U{1)y 
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Table 1 : Vector-like quark multiplets Q mixing with the SM quarks through Yukawa couphngs. The electric charge 

is Q = T3 + y. 



Q 



isospm 

hypercharge 



U 



D 





1/2 1/2 

2/3 -1/3 1/6 7/6 



1/2 1 
-5/6 2/3 



1 

-1/3 



quantum numbers), as 

A + A + A^, (1) 

where Ci is the SM Lagrangian, Ch, contains the kinetic, Dirac mass terms (that can always 
be taken diagonal) and Yukawa couplings involving only the (heavy) vector-like quarks, and 
finally Cih contains the (linear) mass (and Yukawa) mixing between SM and vector-like quarks. 
The general expressions for these Lagrangians can be found in Ref. [232]. In order to obtain 
the properties of the physical particles {i.e. those with a well defined mass) we have to diag- 
onalize the corresponding mass matrices that are contained in Ch + Cih. Due to the unbroken 
electromagnetic U{1)q gauge invariance, the mass matrix will be block diagonal according to 
the charges of the different fields. 



fM^^T + h.c., (2) 



where ipf^^ is a vector that contains the uq quarks with charge Q and M'^ is a ng x uq matrix. 
The overscript (0) denotes current eigenstates. We can then diagonalize each block with two 
unitary matrices U^r, 

{U^yMgUS = Vq, (3) 

with Vq a diagonal matrix containing the masses of the physical particles. C/^^ can be com- 
puted as the unitary matrices that diagonalize the mass matrix squared, {1/^^ M.qMqU^ — 

{U^y MqMqU^ — Vq. Once we have computed the rotation matrices U^ji, we only have to 
replace in the gauge and Yukawa couplings the mass eigenstates in terms of the physical quarks, 

= (4) 

where ^/^^^ are now the physical states. This final step will have two main implications. First, it 
will modify the SM quark couplings to the electroweak gauge bosons and the Higgs, and second 
it will introduce off-diagonal (electroweak) gauge and Yukawa couplings between heavy and 
SM quarks. In general, the first effect strongly contraints large mixing of vector-like quarks with 
first and second generation quarks (as their couplings have been measured to agree very well 
with the SM prediction). ^ Only mixing with the top is at present poorly constrained and could 
be large. If the new vector-like quarks play a role in the resolution of the hierarchy problem, it 



'One exception is the possibility of cancellations between the contributions of several vector-hke quarks, some- 
thing that can happen quite naturally in some models with extra dimensions. 
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is indeed natural for them to be relatively light and to mix sizably with the top. Note that in this 
case they might induce sizable corrections at loop level to precision observables, most notably 
the T parameter and the Zbb coupling, that should be checked in particular models [233-235]. 

Let us make this procedure more explicit with one relevant example whose phenomenol- 
ogy will be further studied in this report. In particular, we will consider an extension of the SM 
with two vector-like quark doublets, ql^r = {<f-,<fTLR and Xl,r = {x^^X'^Tlr^ with hyper- 
charges = 7/6 and Yq = 1/6, respectively. The electric charge of the different components 
are given by Q = + Y (Q^u = 5/3, Qqu = Q^d = 2/3 and Qqd — —1/3). Note in par- 
ticular the exotic charge of x". For simplicity, we will also assume that these two new quark 
doublets can only mix with the top but not with the bottom or any other light quark. This ex- 
ample is motivated by recent ideas using a subgroup of the custodial symmetry to protect large 
corrections to the Zbibi coupling [107] for which these simplifying assumptions are naturally 
realized. This mechanism has been successfully implemented in composite Higgs [234-238] 
and Higgsless models [239] in warped extra dimensions. 

The mass Lagrangian for the charge 2/3 quarks (the other ones do not have any mass 
mixings and are therefore mass eigenstates) can be written as. 



-m(0) 



-d(0) 
Xl 



rrit 



m. 



q,t 





Mq 





/ 4°' \ 

fa™ 



h.c. 



(5) 



EWSB masses, which are proportional to the Higgs vev, are denoted with a lower case m 
whereas Dirac masses, that can be arbitrarily larger than v are denoted by a capital M. The 
hierarchy between these two types of masses, allows us to diagonalize the mass matrix pertur- 
batively, in a power expansion of m/M <^ 1 in order to obtain simple analytic expressions for 
the masses and couplings of the physical quarks. For instance, to leading order in m/M, we 
have 



t 



(0) 
R 



tR 



m 



q,t 



m 



x,t d 
Xr-, 



«(0) 

1r 



m 



Qr 



q,t 



tR, 



x'k 



m, 



" M^^"' Mq M^ 



tR, (6) 



where fields without the (0) superscript are physical (mass eigenstate) fields. All other fields 
are not modified (i.e. they are already mass eigenstates) at this order. Now we only have to 
introduce these rotations in the gauge and Yukawa couplings to obtain the couplings among the 
physical fields. The couplings to the Z are, in the current and mass eigenstate basis, respectively 



9 



2cw 



v«(0)^,At m(0) 



xfVxT 



"^2/3] 



2cw 



m 



tRi'xi - 



+ - 



(7) 



where J^g 



r7<0)^„,(0) 



(2/3)(tl,-7% 



+ xTrxT 



{2/mRl^tR + rRl^tR + 



XrI^Xr) is the electromagnetic current for the RH charge 2/3 quarks and £f is equal to the 
original Lagrangian with the replacement ^ ^. Similarly we obtain, for the charge currents 



AC 



w 



£W _ £W 



m 



x,t 



xirt 



R 



+ h.c. + 



m 



\M^ 



and Yukawa couplings 



H 



C" - 



II 



[mq,tqltR + m^^txitR + h-c] + O 



m' 



(8) 



(9) 
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Diagonal couplings and masses are only modified at 0{m? /M'^) [240] which, for the purpose 
of their collider implications can usually be considered a negligible correction. (Of course if 
the vector-Uke quark masses are close to the electroweak scale, M ~ v, this expansion breaks 
down and a full numerical diagonalization would be required). Using these couplings we can 
compute the decay width of the different heavy quarks. The result is extremely simple in the 
large M limit, in which we obtain 

r(x" ^ wt) ^ 2r(g« ^ zt) ^ 2r(g" ^ m) ^ f — g^, (lO) 

b47r My^ 

Tiq" ^ Wt) ^ 2r(x" ^ Zt) ^ 2T{x'' - Ht) ^ (H) 

Note the 2 : 1 : 1 pattern of decays into W, Z, H as predicted by the Equivalence Theorem. This 
means that, for heavy vector- like quarks, the decay is mainly into longitudinal gauge bosons (or 
equivalently in the Goldstone bosons in a gauge other than the unitary one). 

Vector-like quarks can be pair or singly produced. Pair production is governed by QCD 
and is therefore model-independent. There is a very strong suppression at large masses which 
makes the reach at the LHC ~ 1 — 2 TeV [192, 241]. However, pair production can give rise 
to spectacular signatures. For instance, pair production of or q"^ results in a final state with 
four W plus two h quarks. An analysis of the reach in that case, using like-sign leptons will be 
discussed below. Single production occurs through the off-diagonal couplings between SM and 
heavy quarks. Although it is model-dependent, the suppression at large masses is milder than in 
pair production. Thus, if the new quarks mix sizably with light quarks (including the bottom), 
single production can give a better reach than pair production [242]. 



2. CONCLUSIONS 

Vector-like quarks are a common prediction of models of new physics. They are in general 
compatible with EWPT and can give rise to spectacular signals at the LHC. We have reviewed 
the tools needed to compute their collider implications (masses and couplings) and discussed a 
particular example relevant for realistic composite Higgs and Higgsless models. Other possibil- 
ities have been discussed elsewhere, with special emphasis on the case of a vector-like singlet 
(for a recent review see [243], page 72, and references there in). 
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Part 14 

Top-partner mass reconstruction by using 
jets 

M.M. Nojiri and M. Takeuchi 
Abstract 

At the LHC, the top-partner (T ) and its antiparticle is produced in pairs 
in the Littlest Higgs model with T-parity. Each top-partner decays into 
top quark (t) and the lightest T-odd gauge partner Ah, and t decays into 
three jets. We demonstrate the reconstruction of t decaying hadroni- 
cally, and measure the top-partner mass from the mr2 distribution. We 
also discuss the dependency on four jet reconstruction algorithms (sim- 
ple cone, kt, Cambridge, SISCone). 

1. TOP PARTNERS IN THE LITTLEST HIGGS MODEL WITH T PARITY 
1.1 Productions and decays 

We consider the Littlest Higgs model with T parity (LHT). In the following, we assume the top 
partner is the lightest in the SM fermions' partners. The top partner may be produced in pairs 
and decays at LHC as, 

pp T_T_ UAhAh bW-^bWAHAn ^ 6j + (1) 

This process is similar to the top squark (i) production process in the MSSM. However, the top 
partner production cross section is larger than that of top squarks when the masses are same, 
because the top partner is a fermion. At the LHC, the T_T_ production cross section is 0.171 

pb for the top partner mass mT_ set to 800 GeV. Once the top partner is produced, it decays 
into top and the heavy photon A^ (We set to 150GeV). The branching ratio is 100% when the 
top partner is the next lightest T-odd particle. To identify this process, it is important to identify 
two tops and missing transverse momentum ^j,. 

To simulate the top partner reconstruction we generated 8,550 it* events with HER- 
WIG6.5 (expected number of events for 50fb~^). We use a similar MSSM model point instead 

of the LHT for the event generation. The distribution of production and decay is almost correct 
except for the spin correlations of T_ decays, so we do not simulate top partner polarization ef- 
fects but top polarization is taken into account. The main source of SM background is ti events, 
whose production cross section is 463 pb at tree level. We also generated 4,630,000 ti events 
(for 10fb~^). When we compare signal and background distributions, we simply multiply ti 
distributions by a factor of 5. We ignore the tiZ background because this can be calibrated with 
tiZ{-^ events easily. 

The process was already analysed in [244], where AcerDETl.O [245] was used for the 
detector simulation and jet reconstruction. We also use AcerDETl.O for the detector simula- 
tion. In addtion to that, we interfaced the calorimeter information to FastJet2.2beta [246] which 
allows us to compare different jet reconstruction algorithms (kt, Cambridge, SISCone). To 
compare the four jet algorithms under the same conditions (in Sec 2.), we switch off jet energy 
smearing of AcerDET. In this section, we only show the results with the Cambridge algorithm. 
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generated 


Mcft- cut 


"Icp = 


PT,H cut 








both JTt 


relaxed m 




8,550 


6,590 


4,384 


2276 


1433 


437 


380 


118 


708 


tt 


23,150,000 


199,640 


88,540 


9475 


6835 


2105 


765 


235 


1835 


s/Vn 


1.777 


14.75 


14.73 


23.38 


17.33 


9.525 


13.73 


7.70 


16.53 



Table 1: Summary of the cuts. The number of events are for JdtC = 50 fb^^. pt,h denotes the cuts that requires 
both hemispheres' momenta Pt,Hi,Pt,H2 > 200GeV. mp^. denotes the cuts that the hemisphere i's momentum 
is consistent with top mass (ISGGeV < mp„ < 190GeV). The relaxed mp^ cut denotes that both hemispheres' 
momenta satisfy 50GeV <mp^ < IQOGeV. 



1.2 Event selection and Top mass reconstruction 

We apply successive cuts to select T_T_ events, the summary of the cuts is shown in Table 1. 
First, we imposed the standard cut to collect events related with new physics, 

^j. > 200GeV and > 0.2Mes, niepton = 0. (2) 

The lepton cut is for dropping ti production events, in which large I^t is dominantly caused 
by neutrinos from the leptonic decay of top. We do not assume any b tagging, because for the 
given kinematics jets are very coUinear with each other, and the b tagging efficiency has not 
been studied in that case. 

We applied the hemisphere analysis to find top candidates [247]. Each of the jets was 
assigned to hemispheres which were defined as follows; 

Vi e ifi, j e H2 d{pH^,Pi) < d{pH2,Pi) and d{pH2,Pj) < d{pH,,Pj)- 

d{pi,p2) = / J^^^ ^^f^^^ is the angle between pi and P2). (3) 
+ 1^2) 

Here, pn^ = Y^i^H-Pi^ required both hemispheres' transverse momenta to be larger 

than a threshold, 

Pt,h„Pt,h, > 200GeV. (4) 

Top quarks from T_T_ production are highly boosted. Therefore decay products from 
two top quarks are correctly grouped into two hemispheres with high probability. We analysed 
the events in which each hemisphere has up to 3 jets to drop contribution from other QCD jets. 

The distributions of invariant masses of hemispheres (mp^ = a/p|^) for the T_T_ and ti 
events are shown in the Fig 1. We can see the top mass peak both for the T_T_ and ti events 
in the Hi. On the other hand, the top mass is not well reconstructed in the H2 for ti events. 
This is because at least one of tops must decay leptonically to make large ^t- Two dimensional 
scattering plots of m vs. m are also shown in Fig 2. We regard a hemisphere's momentum 
as the top momentum in this paper if its mass is consistent with the top mass (ISOGeV < mp^ < 
190GeV). The matching of a hemisphere's momentum with the parton momentum is discussed 
in Section 2.1. 

1.3 Measurement of the end point of mra distribution. 

To extract the mass of T_, we used the Cambridge rriTi variable [248]. This variable is consid- 
ered in events like CC' ~^ {aa){bp), where we assume ( and (' have the same mass m^, a and 
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Figure 1: The distributions of invariant masses of a) iJi for the T_T_ Figure 2: nip^^ vs. mp^^ for 
events, b) H2 for the T_T_ events, c) Hi for the ti events, d) H2 for T_T_ (upper figure) and for ti 
the if events. (lower figure). Cambridge algo- 

rithm is used for Jet reconstruction. 




Figure 3: The niT2 distribution for the Cambridge algorithms, where 
proper m,^o = 150 GeV is used. The dashed line shows for events 
with 150GeV < mp^. < 190GeV The solid fine shows for events 
with 50GeV < mp„^ < 190GeV. The endpoints are 790.95 GeV for 
nominal niAH- (tot_ = 800 GeV). 



h are visible objects, and a and (5 are invisible particles and have the same mass m^o. In the 
events, the mT2 variable is defined as follows: 

^r2(PT>PT>I^Ti"^x°) = -"^^^ max < m^(p^, pJ^; m^o), m^(p^, pJ^; m^o) > . (5) 

Here, the transverse mass m-r is defined as m|,(py, pi^; m^") = m^+m?(, + 2 [E^]^'^ — Pxl^^] ■ 
By the definition, the mT2 variable is a function of m^o, and satisfies the condition mT2{'m^o) < 
m^. Then we can measure by measuring the upper endpoint of the 171x2 distribution. 

For our purpose, visible particles are two top quarks and invisible particles are two Ah. 
We may regard hemispheres' momenta as top momenta for the hemisphere whose masses are 
150GeV < mpjj^ < 190GcV. However, there are not enough events left under this cut. There- 
fore we use a relaxed criterion 50GeV < mp^. < 190GcV. Indeed, when nip^ is less than 
rrit, all of the top decay products do not contribute to the hemisphere reconstruction due to the 
minimum jet energy cuts. The endpoint of 771^2 does not change because uit is always under- 
estimated compared to the true m^. This is easy to understand if you consider the system of 
a neutralino and other sources of missing momenta ( for example neutrino ) as one invisible 
particle. The invisible particle mass (minvisibie) is always larger than the neutralino mass which 
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results mmT2{m^o) < ("^invisible) < "Ne- 
llie mT2 distribution is shown in Fig 3, where the nominal value m^o = 150 GeV is 
used for the calculation of mT2- We show only the region of mT2 > 350GeV. The tt events 
have low mT2 values, so we can neglect the BG to fit the endpoint. The dashed lines show the 
mT2 distributions for events with ISOGcV < mp^. < 190GeV and have the same endpoint as 
those with the relaxed criterion of 50GeV < mp^^ < 190GeV (solid line). This fact supports 
the validity to relax the criterion to determine the endpoint. We fitted the distribution near the 
endpoint by a linear function and obtained mj^f^ = 790.95 GeV. This value is consistent with 
the proper value mT_ — 800 GeV. 

It is not possible to determine the top partner mass itself unless the A h mass is determined. 
If we assume the Ah thermal relic density is consistent with the dark matter density in our 
universe, tuah is related to the Higgs mass (niH) so that it is determined with two fold ambiguity 
[244]. 

2. COMPARISON AMONG JET RECONSTRUCTION ALGORITHMS 

We now make a comment on the dependence on the jet reconstructing algorithm. Four algo- 
rithms for jet reconstruction (Simple cone, kt, Cambridge, SISCone) are used in the following 
analysis. For this study, we interfaced FastJet2.2beta [246] to AcerDETl.O. The energy deposits 
in calorimeter cells are regarded as "particle momentum". Cone sizes (or the counterparts for 
clustering algorithms) R are chosen as 0.4 and overlap parameter / as 0.5 for SISCone. Simple 
cone denotes the algorithm used in AcerDETl.O. 

2.1 Invariant mass distributions 

The results are summarized in Table 2. Acerdet has an option to calibrate jet momenta, therefore 
both calibrated and non-calibrated numbers are given for the simple cone algorithm in the table. 
The calibrated result is shown in the figures. 

The distributions of mp^^ are shown in Fig 4. The shaded regions denote that the hemi- 
sphere's invariant mass is consistent with mt (150GeV < mp^ < 190GeV). The peak for 
the Simple cone algorithm is dull and has a broad tail and the position of the peak is located 
below rrtt. This is because AcerDET takes massless jets. Therefore the probability to find the 
jets consistent with the top mass is relatively low among the four algorithms. The acceptance 
becomes worse as the lower bound of an invariant mass cut for nip^ is increased. If the other jet 
definitions are feasible in the LHC environment, the reconstruction efficiency can be increased 
significantly. On the other hand, the endpoints of the mT2 distributions are 795.69 (779.37) GeV 
(Simple cone), 780.62 GeV (kt), 790.95 GeV (Cambridge), 803.67 GeV (Siscone) respectively. 
The dependence on jet algorithms is not severe for this analysis. 

Fig 5 shows the deviation of the hemisphere momentum from the true top quark momen- 



For r_T_ events 


^, Meff cut 


niet,H < 3 






both nipjj 


relaxed nip^ 


mT2 endpoint 


Simple Cone (calibrated) 


6673 


945 


306 


326 


110 


439 


795.69 GeV 


(not calibrated) 


6673 


1213 


283 


323 


82 


545 


779.37 GeV 


kt 


6673 


1436 


444 


384 


105 


604 


780.62 GeV 


Cambridge 


6673 


1433 


437 


380 


118 


621 


790.95 GeV 


SISCone 


6673 


1656 


512 


437 


150 


608 


803.67 GeV 



Table 2: Sunmiary of the cuts for various jet algorithms, cut is the same as Table 1. 
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turn {Apt)/pt- We selected signal events with 150 GeV< rrip^^ < 190 GeV for the Simple 
cone and Cambridge algorithms. They are mostly distributed in a ±5% region. We can see the 
(Apr) /pt are larger than by about 2% for the Simple cone. The positive contribution is due 
to jet pt calibration of AcerDET. 




Figure 4: The distributions of mp^^ for the T^T^ 
events for a) Simple Cone (Acerdet), b) kt, c) Cam- 
bridge and d) SISCone. Acerdet has dull peak and broad 
tail and the position located below rrit 



Figure 5: {Apt)/pt dis- Figure 6: The comparison 

tributions of the hemisphere between Cambridge (Solid) 

Hi selected by 150 GeV< and kt (dashed) algorithms 

mpjj^ < 190 GeV for with underlying events. 

AcerDET with jet calibra- Cone sizes of jets are 

tion (upper) and Cambridge AR = 0.4 (upper) and 0.6 

(lower). (lower). 



2.2 Effect of Underlying Events 

So far we have ignored the effect of underlying events when generating events. We now inves- 
tigate the reconstruction with underlying events by comparing the kt algorithm and Cambridge 
algorithm. We generated the underlying events using JIMMY with HERWIG6.5. Fig 6 shows 
the distributions of invariant masses of P^i for the kt and Cambridge algorithms with underly- 
ing events. The event selection cuts are the same as in section 1.2. We can see that the position 
of the peak for the kt algorithm is larger than that for the Cambridge algorithm and the distri- 
bution is smeared for increasing R. This is because the kt algorithm over-collects contributions 
which are far from the jet direction (large Rij) due to the factor min(A;^-, k'fj) in the definition 
of distance (splash-in effects) . The Cambridge algorithm does not have the factor min(A;|, kfj) 
therefore it is not too sensitive to the existence of underlying events. 

CONCLUSIONS 

We have reconstructed top quarks from T T production and the decay into top and stable 
gauge partner in the LHT for the top quark decaying hadronically. The main background from 
SM processes is ti production events. They can be reduced by applying the hemisphere analysis 
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and imposing a cut for hemisphere momenta. The top partner mass can be measured by using 
mT2 variable. 

We have also investigated the dependence on the jet reconstructing algorithms. AcerDET 
takes massless jets, so invariant masses are significantly underestimated for the boosted top 
quark. We therefore present our result for the Cambridge jet reconstruction algorithm. We 
find that the kt algorithm over collects far activities and the invariant masses are overestimated 
(Splash-in effect) compared with the Cambridge algorithm. 
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Part 15 

Searches for a paraphoton in associated 
production with tt in e+e~collisions 

E. Boos, V. Bunichev and H.J. Schreiber 
Abstract 

We discuss prospects to search for a new massless neutral gauge boson, 
the paraphoton, in e+e~ collisions at center-of-mass energies of 0.5 and 
1 TeV. The paraphoton naturally appearing in models with abelian ki- 
netic mixing has interactions with the Standard Model fermion fields 
being proportional to the fermion mass and growing with energy. At 
the ILC, potentially the best process to search for the paraphoton is 
its radiation off top quarks. The event topology of interest is a pair 
of acoplanar top quark decaying to jets and missing energy. Applying 
a multivariate method for signal selection expected limits for the top- 
paraphoton coupling are derived. Arguments in favor of the missing 
energy as the paraphoton with spin 1 are shortly discussed. 

1. INTRODUCTION 

Although the Standard Model does not require any additional gauge bosons it is possible to 
introduce gauge invariant operators in the Lagrangian which involve new gauge fields. An 
example is given in [249] (see in addition [250-252]) by the abelian kinetic mixing of the SM 
Uy{1) field with a new Up{l) field in a gauge invariant manner. The mixing term of the two 
U{1) fields can be diagonalized and canonically normalized by an SL(2,R) transformation 
in a way that one linear combination of the fields corresponds to the ordinary photon which 
couples in the usual manner to all electrically charged particles within the SM. The other linear 
combination appears as a massless spin-1 neutral particle, referred to as the "paraphoton" in 
[253-255] and denoted by 7' in this note. The paraphoton couples only indirectly to the SM 
fields via higher dimensional operators as was worked out in [256]. The effective interactions 
of the paraphoton with the SM fermions following from the higher dimensional operators have 
the chirality flip structure, proportional to the SM fermion masses and inversely proportional to 
new physics scale M squared: 

^F;: (qLa^'^C^HuR + QLa^^CdHdR + ha'^'^CeHeR + h.cj , (1) 

where F^l is the paraphoton field strenth, Ql^Il are the quark and lepton doublets, ur, (Ir the up 
and down-type SU (2) singlet quarks, cr the electrically-charged SU (2)-singlet leptons, and H 
is the Higgs doublet. The dimensionless coupling parameters C/ = C„, C^, after the Higgs 
field gets its vev Vh are re-expressed in the form Cf = Cf*mf/{vh/ a/2) . The coefficients c / are 
unknown, but various phenomenological constraints exist. Discussions on possible lower limits 
on 7' interactions with fermions can be found in ref. [256]. 

As follows from the Lagrangian (1) and existing bounds on the couplings Cf for light 
fermions the paraphoton couples most strongly to the top quark and very weakly to light 



103 



M/^t [TeV] 


^/s = Q.bTeV 


= 1 TeV 


0.2 


5700 


42500 


0.3 


1100 


8500 


0.5 


40 


1100 


1 


10 


70 



Table 1 : tt"f' event rates for several values of M/ ^Jct at \fs = 0.5 and 1 TeV and an integrated luminosity of 0.5 
respectively 1 ab~^. 

fermions. Therefore, one expects the most interesting process to search for the paraphoton 
is 7' radiation off the top. Since so far no constraint on ct exists, access to Mj ..Jct seems pos- 
sible or corresponding limits might be set for the first time. It seems a priori very difficult to 
perform 7' searches at hadron colliders because of very large tt + multi-jet background. The 
next generation e^e~ linear collider (ILC) provides potentially a possibility to search for the 
paraphoton via the channel 

e+e~ ^tti . (2) 

2. THE SIGNAL REACTION 

The parameter scan was done for various couplings and detailed simulations of ti'^' signal events 
were performed for a 'reasonable' value of the effective coupling parameter Mj y/ct {Mj y/ct = 
0.2 TeV) when the signal is large enough to be clearly distinguishable from the SM background. 
All the computations and simulations were performed for the ILC collision energies ^/s — 0.5 
and 1.0 TeV and an integrated luminosity of 0.5, respectively, 1 ab~^. 

The characteristics of the signal reaction were computed and partonic events were gen- 
erated by means of the program package CompHEP [257, 258]. The Feynman rules for the 
fermion-fermion-7' vertices following from the effective Lagrangian (1) 

§^-mfp2'{Yr-rr) (3) 

have been implemented into CompHEP. An interface with PYTHIA 6.202 [259] allows one to 
simulate initial and final state radiation and jet hadronization. Also, beamstrahlung effects [260] 
are taken into account. The signal event rate is given in the Table 1. In order to estabUsh a search 
strategy for the paraphoton in tt events it is advantageous to know whether an off-shell or on- 
shell top quark radiates the 7'. Fig. 1 (left) shows the invariant mass of the 'y'Wb system of that 
top which radiates the paraphoton. Clearly, in most cases the paraphoton is radiated off a top 
being off-shell, and 7' search strategies should be based on on-shell top with t Wb decays in 
association with the 7'. The energy of the 7' shown in Fig. 1 (right) reflects that the paraphoton- 
top coupUng is proportional to the paraphoton momentum, so that large missing energy, ^, will 
be a tag of signal events. 

3. SIGNAL EVENT SELECTION 

After event generation using CompHEP, PYTHIA and the CompHEP-PYTHlA interface pack- 
ages with the Les Houches Accord implemented [261] an approximate response of an ILC 
detector was simulated by means of SlMDET_v4 [262]. 
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Table 2: Background events at \/s = 0.5 and 1 TeV for an integrated luminosity of 0.5, respectively, 1 ab ^. 

The most important background consists of ti + (7) events, where photons from initial 
state radiation (ISR) are not detected. The number of events expected for both energies are given 
in Table 2. They exceed substantially the number of signal events (see Table 1) for the chosen 
value for the parameter M/ y^. The next significant background to consider is e'^e~ ti+uu, 
with the same signature as for the signal. The corresponding event numbers also given in Table 
2 are comparable to the signal event rates for not too small M j values. An invariant mass 
cut of e.g. Myjy < 80 GeV, i.e. a cut on the event missing mass, removes most of these events. 

To discriminate signal and a very large background we use a multivariate technique based 
on likelihood method. 18 kinematics variables, such as missing energy, invariant masses of 
various jet combinations etc. were combined into a global discriminant variable Pp, designed 
to give a measure of the 'Paraphoton-likeness' of any particular event. As it should be the 
background events are preferentially distributed at low Pp values while for signal events Pp is 
concentrated close to unity. By choosing optimal values of cuts for the discriminant Pp one gets 
the signal selection efficiency of 49% (76%) at y/s = 0.5 (1) TeV, while only 9% of background 
events survive. At = 0.5 TeV, S/Vb =11.96 for M/^ = 0.2 TeV, while S/v^ = 162.6 
at 1 TeV, i.e. the probability of measuring the total event rates as a result of a background 
fluctuation is 0.5 ■ 10^^^ and < 10^^^ at 0.5, respectively, 1 TeV, using Gaussian sampling of 
uncertainties. In this way, an almost background-free signal event sample can be extracted for 
further measurements at 1 TeV. The situation is much less convenient at 0.5 TeV. 

4. DISCUSSION OF THE RESULTS 

If an excess of signal events over the SM background is established, limits on the inverse cou- 
pling parameter M/y/ct accessible for a significance of S/\fB = 5 can be derived being suf- 
ficient for the paraphoton discovery. The numbers of surviving 7' events for 5a discovery at 
0.5 and 1 TeV energies and an integrated luminosity of 0.5 and 1.0 ab^^ can be converted into 
limits for M/ ^fct. These limits are found to be of 0.33 and 0.61 TeV for 0.5 and 1 TeV cases 
respectively. The value Mj ^Jct = 0.61 TeV, is expected to be the most stringent limit accessible 
at the ILC. 
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The signal-to-background ratio, S/ B, is about 1.79 at 1 TeV, sufficiently large to under- 
stand a spin assignment of the radiated massless particle. In order to demonstrate the spin-1 
nature of the 7', we follow studies performed to establish the vector nature of the gluon in 3-jet 
e^e' annihilation events at PETRA [263-266] and LEP [267-269] energies, based on predic- 
tions that a spin-^ quark radiates the spin-1 gluon. The analysis of the EUis-Karliner angle [270] 
distribution and the polar angle distribution of the normal to the reaction plane at 1 TeV clearly 
shows that spin-1 assignment for the paraphoton is highly favored over spin 0. 

5. CONCLUSIONS 

Some realistic extensions of the Standard Model suggest the existence of a new massless neu- 
tral gauge boson, denoted as the paraphoton 7' in this study. This particle is similar to the 
ordinary photon, but the interactions of the 7' are very distinct: couplings to SM fermions are 
proportional to fermion masses and therefore strongest to the top quark, and grow with the 
7' momentum. Hence, the paraphoton radiation off the top at the ILC is studied as the most 
promising process for the paraphoton observation. Only the all-hadronic top decay mode was 
selected to ensure a high signal-to-background ratio and to avoid complications due to final state 
neutrinos in leptonic W decays. A multivariate search strategy was used to better separate the 
signal from backgrounds. Allowing for a 5a paraphoton discovery significance, e+e" collisions 
at 1 TeV allow to bound the 7' -top quark inverse coupling to M/ ^Jct ^ 0.61 TeV, which is ex- 
pected to be the most stringent limit accessible at the next generation colliders. For the sake of 
demonstration two angular variables, the ElUs-Karhner angle and the polar angle of the normal 
to the t i 7' plane as a function of a thrust cut-off, were studied to establish the vector nature 
of the 7'. Both angular distributions are in accord with the spin-1 assignment of the paraphoton 
and inconsistent with e.g. a scalar hypothesis. 
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Part 16 

High Hadronic Top Quark 
Identification 

G. Brooijmans 

Abstract 

At the LHC objects with masses at the electroweak scale will for the 
first time be produced with very large transverse momenta. In many 
cases, these objects decay hadronically, producing a set of collimated 
jets. This interesting new experimental phenomenology requires the 
development and tuning of new tools, since the usual reconstruction 
methods would simply reconstruct a single jet. This study describes the 
application of the YSplitter algorithm in conjunction with the jet mass to 
identify high transverse momentum top quarks decaying hadronically. 



1. INTRODUCTION 

At the LHC, top quarks, W and Z bosons are relatively light and can be produced with very 
high transverse momenta with respect to their masses. In the case of hadronic decays, the 
quarks can be so close together in the detector that they are in principle reconstructed as a 
single jet. This is illustrated in Figure 1, which shows generator-level distributions of angular 
distance between the b quark and W boson, and quarks from W boson decays in tt events. Here 
dR = a/ (A0)2 + (A?])^ with (j) the azimuthal angle and r] the rapidity. For top quark transverse 
momenta larger than about 200 GeV, the distance between the decay products is often smaller 
than twice the typical jet radius. These events were generated using PYTHIA [271] and no top 
polarization effects are included. Such effects are very model dependent and therefore beyond 
the scope of this generic study. 




1000 1200 
Top pT (GeV) 




800 1000 1200 
Top pT (GeV) 



(a) 



(b) 



Figure 1: Angular distances between decay products in top quark decays as a function of top quark transverse 
momentum: (a) between the b quark and W boson, and (b) between quarks from W boson decays. 
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Identifying top quarks at high transverse momentum with high efficiency is of particular 
interest in searches for new physics. In addition to a number of recent theoretical models specif- 
ically proposing the existence of high mass resonances decaying dominantly to top quarks (see 
for example [104, 133, 142, 176, 244]), the large top quark mass suggests it might be closely 
linked to forms of new physics that would manifest themselves at very high energies. It is 
therefore quite probable that new heavy objects decay to top quarks at least some fraction of the 
time, if not exclusively. 

2. DATASETS AND TOOLS 

The datasets used in this study are a) sequential standard model Z'-bosons of mass 2 and 3 
TeV decaying to top-antitop pairs for signal, and b) multijet events with transverse momenta 
ranging from 300 to 2200 GeV for the background. For the signal events, one of the top quarks 
is forced to decay hadronically, and the other semileptonically. The events were all generated 
using PYTHIA, passed through the full ATLAS detector simulation and reconstructed using the 
ATLAS reconstruction program. 

The jet mass, which is the invariant mass of all the jet's constituents (typically calorimeter 
cells or towers), and "YSplitter" [214], which determines the scales at which jets can be resolved 
into two or more subjets are used as discriminating variables. Top "monojets" are selected in 
the signal samples by selecting events in which only one jet with pt > 20 GeV has dR < 1.0 
from the closest top quark and dR < 2.0 from the hadronically decaying W boson. For these 
jets. Figure 2 shows the distribution of jet mass as a function of jet transverse momentum, and 
the scales at which the jet splits up into two and three jets. The distributions have only limited 




Jet pT (GeV) YScale 1-2 (GeV) YScale 2-3 (GeV) 



(a) (b) (c) 

Figure 2: (a) Jet mass as a function of transverse momentum for jets passing the top monojet selection: jets from 
the M = 2 (3) TeV Z' sample in blue (red), (b) and (c) Solid (dashed): scales at which the top monojet splits into 
two and three jets respectively for events in the AIz' — 2(3) TeV samples. 

dependence on the jet transverse momentum distribution and are therefore well suited to the 
identification of top monojets over a wide spectrum. The splitting scales into two and three jets 
cluster around half the top quark mass and half the 11/ -boson mass respectively as expected. 
The same distributions are shown for the background samples in Figure 3. 

3. QUANTITATIVE ANALYSIS 

In principle, with this number of variables a multivariate tool like an artificial neural network 
might yield optimal results in terms of high/low efficiency for signal/background. For the sake 
of clarity however this study is based on simple two-dimensional cuts, keeping in mind that the 
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Figure 3: (a) Jet mass as a function of transverse momentum, (b) and (c) scales at which the jet splits into two and 
three jets respectively for events in the background samples. 



use of multivariate tools typically leads to a factor of approximately 1.5 improvement in the 
signal over background ratio. 

Cuts are applied on 

1. Jet mass (> 170 GeV), 

2. jet mass as a function of jet transverse momentum, 

3. YScale 1-2 (split from one to two jets) as a function of YScale 2-3 (split from two to three 
jets), 

4. YScale 2-3 as a function of YScale 3-4, 

5. YScale 1-2 as a function of jet mass, 

6. YScale 2-3 as a function of jet mass, and 

7. YScale 3-4 as a function of jet mass (two cuts). 

As an illustration, three of these (2, 4 and 7) are shown in Figure 4. The resulting efficiencies 
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Figure 4: Cuts applied on (a) jet mass as a function of jet transverse momentum (events are required to he below 
the line), (b) YScale23 as a function of ^5*00/634 (events are required to lie above the line), and (c) YScale34 as 
a function of jet mass (events are required to lie above the line). In each plot, the background is in blue and the 
signal in red. 



for signal and background events are shown in Figure 5. For the background the efficiency 
plateaus at about 10% starting at jet transverse momenta of 1300 GeV, where the signal effi- 
ciency reaches 60%. The ratio of signal over background efficiency increases for smaller jet 
transverse momenta. A more detailed description of this analysis is available in [272]. 
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Figure 5: Selection efficiency as a function of transverse momentum for top monojets (solid, black) and jets in the 
background samples (blue, dashed). 



4. CONCLUSIONS 

In this analysis, the measured jet mass and jet splitting scales have been used to distinguish 
high transverse momentum "top monojets" from jets originating from light quarks. The com- 
bination of algorithms allows for good separation of signal and background, with the ratio of 
selection efficiencies for signal and background evolving from approximately 30 for jets with 
Pt = 600 GeV to 10 for 1000 GeV and 7 for 1500 GeV. Further work using subjets and tracking 
information is underway. 
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Part 17 

A Les Houches Interface for BSM 
Generators 

7. Alwall, E. Boos, L. Dudko, M. Gigg, M. Herquet, A. Pukhov, P. Richardson, A. Sherstnev and 
P. Stands 

Abstract 

We propose to combine and slightly extend two existing "Les Houches 
Accords" to provide a simple generic interface between beyond-the- 
standard-model parton-level and event-level generators. All relevant in- 
formation — particle content, quantum numbers of new states, masses, 
cross sections, parton-level events, etc — is collected in one single file, 
which adheres to the Les Houches Event File (LHEF) standard. 

1. INTRODUCTION 

The simulation of interactions at the LHC is characterized by the use of many different pro- 
grams specializing in different stages of the calculation, such as matrix-element-level event 
generation, decay of resonances, parton showering, hadronization, and underlying event simu- 
lation. The communication of simulation parameters between those stages can be complicated 
and program-specific. For supersymmetric models, this situation has been greatly improved 
by the introduction of the SUSY Les Houches Accord [273] (SLHA) and its upcoming exten- 
sion [274]. For general models however, there is still no corresponding standard. In this note, 
we suggest an addition to the SLHA to allow for the specification of the quantum numbers, 
masses, and decays of arbitrary new states, thus generalizing the accord beyond its original 
supersymmetry-specific scope. We also make a proposal for how to include these model pa- 
rameter files into Les Houches Accord event files [275] (LHEF) in a standardized way. This 
both reduces the number of files that need to be passed around and minimizes the possibility for 
error by keeping all relevant model information together with the actual events. 

2. DEFINITION OF THE INTERFACE 

The concrete proposal consists of the following three points: 

1. Introduce new SLHA-like blocks QNUMBERS (for "quantum numbers") with the format: 

BLOCK QNUMBERS 7654321 # balleron 

1 # 3 times electric charge 

2 1 # number of spin states (2S+1) 

3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 

4 # Particle/Antiparticle distinction (O=own anti) 
where this example pertains to a fictitious neutral spin-0 color-singlet self-conjugate par- 
ticle to which we assign "PDG" code 7654321 and the name "balleron". That is, the 
BLOCK declaration should define a PDG code and, optionally, a human readable name 
after the # character (if no name is given, the PDG code may be used). We advise to 
choose PDG numbers in excess of 3 million for new states, to minimize the possibility 
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of conflict with already agreed-upon numbers [85]. The entries so far defined are: 1: the 
electric charge times 3 (so that most particles will have integer values, but real numbers 
should also be accepted); 2: the particle's number of spin states: 25" + 1; 3: the colour 
representation of the particle, e.g., 1 for a singlet, 3 (-3) for a triplet (antitriplet), 8 for an 
octet, etc.; 4: particle/antiparticle distincition, should be (zero) if the particle is its own 
antiparticle, or 1 otherwise. 

2. Use the existing SLHA blocks MASS and DECAY [273] to define particle masses and 
decay tables. If the model in question is a SUSY model, a full SLHA spectrum [273] can 
also be included. We propose that the reader should "turn on" SUSY whenever the SLHA 
SUSY model definition block MODSEL is present. 

3. Include the information from points 1 and 2 enclosed within the subtags < s lha> < / s lha> 
in the <header> part of Les Houches event files [275]. 

3. IMPLEMENTATIONS 

For the purpose of this contribution, the above proposal was tested explicitly by interfacing 
MadGraph/MadEvent with Pythia. Below we summarize the main aspects of these im- 
plementations. 

3.1 MadGraph/MadEvent implementation 

Starting from version 4 [35], the multi-purpose MadGraph/MadEvent parton-level event 
generator by default includes a detailed summary of all simulation parameters in the output 
LHEF [275] parton-level event file. From version 4.1.47, this information is stored in the XML 
<header> section. For the interface considered here, the relevant part of this section is a copy 
of the so-called param.card . dat MG/ME input file. 

The MG/ME param_card.dat uses an extension of the SUSY Les Houches Accord 
[273,274] for model parameters in all implemented models. In particular, it always includes the 
SMINPUTS, MASS, and DECAY blocks. This file is used by MadGraph/MadEvent as an 
input for cross section computations and event generation but is not modified by the program. 
The file is instead assumed to be created by an external "Model Calculator". Such calculators 
are currently available on the web for the SM, MSSM and 2HDM models. Starting from the 
parameters in the Lagrangian (primary parameters), they calculate all needed secondary param- 
eters (such as masses, decay widths, and auxiliary parameters). Note that widths and branching 
ratios can also be evaluated in an intermediate step by MG/ME itself or by external tools like 

DECAY or BRIDGE [84]. 

In previous versions of MadGraph/MadEvent, the param_card. dat file did not 
contain information regarding the particle content of the physical model considered. This infor- 
mation is stored in the particle . dat file filled by model writers during the model creation. 
Starting from version 4. 1 .43, the template for inclusion of user defined models (called USRMOD) 
in MadGraph/MadEvent automatically generates the QNUMBERS blocks described above 
from the information contained in the particle . dat file. These blocks are then included in 
the default param_card . dat for the new model (and from there are copied into the LHEF 
output), such that no extra intervention is required to pass them to parton shower programs af- 
ter parton-level event production. The script only outputs information for particles which have 
PDG numbers not identified as standard SM or MSSM particles, since those are assumed to be 
defined in the parton shower generators. 
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Note that in the current version, the spin, color and particle/antiparticle information is 
automatically extracted, but not the electric charge, which is set to zero by default. This is due 
to the fact that, in MadGraph/MadEvent, the electric charge does not appear in the list of 
particle properties and is only defined through the value of the coupling to the photon. This 
issue will be addressed in future versions of USRMOD, but can currently be circumvented by 
fixing the electric charge information by hand at the end of the model implementation process. 



3.2 Pythia implementation 

The following capabilities are implemented in PYTHIA 6.414 [75] and subsequent versions. 

Already for some time it has been possible to use the QNUMBERS blocks described above 
to define new particles in Pythia via its SLHA interface [276]. What is new is that, when 
reading an LHEF event file, Pythia now automatically searches for QNUMBERS blocks in the 
header part of the LHEF file, updating its internal particle data tables accordingly. It then pro- 
ceeds to search for MASS and DECAY tables, and finally looks for other SLHA blocks contained 
in the header. If the SUSY model definition block MODSEL is found, SUSY is automatically 
switched on and the remaining SLHA blocks are read, without the user having to intervene. 
The read-in of LHEF files containing general BSM states, masses, and decay tables, should 
therefore now be relatively " plug-and-play". 

A note on decay tables: only 2- and 3-body decays can currently be handled consistently. 
They are then generated with flat phase space, according to the branching ratios input via the 
DECAY tables. The colour flow algorithms have been substantially generalized, but if too many 
coloured particles are involved (e.g., an octet decaying to three octets) Pythia will still not 
be able to guess which colour flow to use, leading to errors. Please also read the warnings 
in the section on decay tables in the SLHA report [273] concerning the dangers of double 
counting partial widths and obliterating resonance shapes. To get around the restriction to 
flat phase space, either 1) use Pythia's internal resonance decays whenever possible (e.g., 
do not read in decay tables for particles for which Pythia's internal treatment is not desired 
modified), 2) perform the decays externally, before the event is handed to Pythia (e.g., with 
MadGraph/Bridge [84] or CalcHEP [277]), or 3) do a post facto re-weighting of the 
generated events, based on the kinematics of the particle decays stored in the event record. 

The interfaces can of course still also be used stand-alone, independently of LHEF. The 
user must then manually open a spectrum file containing QNUMBERS and MASS information 
and give Pythia the logical unit number in IMSS (21). New states can then be read in via 
either of the calls 

CALL PYSLHA (0, KF, IFAIL) ! look for QNUMBERS for PDG = KF 

CALL PYSLHA (0, 0, IFAIL) ! read in all QNUMBERS 

and MASS information can be read by 

CALL PYSLHA (5, KF, IFAIL) ! look for MASS entry for PDG = KF 

CALL PYSLHA (5, 0, IFAIL) ! read in all MASS entries 

where IFAIL is a standard return code, which is zero if everything went fine. (For read-in 
of a complete SLHA SUSY spectrum file, these direct calls should not be used, instead set 
IMSS ( 1 ) =11 before the call to PYINIT.) For stand-alone decay table read-in, the unit num- 
ber of the SLHA decay table file should be given in IMSS ( 2 2 ) , and the corresponding read-in 
calls are 

CALL PYSLHA (2, KF, IFAIL) ! look for DECAY table for PDG = KF 

CALL PYSLHA (2, 0, IFAIL) ! read in all DECAY tables 
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4. CONCLUSIONS AND OUTLOOK 

We have proposed a simple file-based interface between parton- and event-level generators fo- 
cusing on the particular problems encountered in the simulation of beyond-the-standard-model 
collider physics. To deal with general BSM models, we add a new block QNUMBERS to the 
SLHA structure, which defines the SM quantum numbers of new states for use in subsequent 
resonance decay, parton showering, and hadronization programs. We also integrate the SLHA 
file into the existing LHEF format to minimize the number of separate files needed. The 
proposal has been tested explicitly by implementations in the MadGraph/MadEvent and 
Pythia6 Monte Carlo event generators. 

In the near future, also the HERWIG++ [11] and PythiaS [278] generators will be ex- 
tended to automatically read in SLHA spectra from LHEF headers. Likewise, forthcoming ver- 
sions of the CalcHEP [34] and CompHEP [258] parton-level generators will include write-out 
of this information in their LHEF output, including also the QNUMBERS extension. 

In the longer term, with the XML format emerging as the de facto standard for file-based 
interfaces, we note that it could be worth investigating the merits of formulating an XML-SLHA 
scheme, that is, transforming the current ASCII SLHA format conventions into a native XML 
form that could be parsed with standard XML packages. A concrete first realization of such a 
strategy is HepML [279] which aims to unify the description of generator information in the 
form of standard XML schemes, in which an XML-SLHA scheme would form a natural part. 
The first release of the public HepML library has been implemented into CompHEP version 
4.5, including also HepML headers in the LHEF output. 
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